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PREFACE 


This  report  documents  work  conducted  under  Task  Number 
3  of  Contract  Numbur  D0T"CG-72O7^,*A  from  January  I/,  1978  to 
August  21,  1979.  The  work  was  performed  at  Clemson  University 
under  the  auspices  of  the  U.  S.  Coast  Guard,  with  LTjg  Steven 
F.  Wlker  and  Ens.  John  A.  Budde  serving  as  program  technical 
monitors.  The  contract  principal  Investigator  was  Dr.  R.  Michael 
Harnett.  The  leader  of  the  part  of  the  work  addressed  In  this 
report  was  Dr.  E.  R.  Baker,  !V.  The  graduate  assistant  In  this 
work  was  Jeffrey  L.  Rlnguest. 
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1.0  INTRODUCTION 


The  use  of  human  subjects  for  testing  hypothermia  protection  devices 
Is  not  without  risk.  Additionally,  such  In  vivo  tests  can  be  used  only 
to  mild  state  of  hypothermia  (rectal  temperatures  not  less  than  35°C).  A 
mathematical  model  capable  of  accurately  simulating  the  thermal  responses 
of  a  protected  man  In  a  cold  environment  would  be  an  attractive  alternative 
to  human  experimentation.  In  addition,  the  time  required  to  perform  the 
human  experiments  makes  this  alternative  method  of  device  evaluation 
extremely  attractive. 


I . I  Background 

A  variety  of  models  have  been  proposed  for  the  human  thermal  system. 

Some  of  the  models  have  been  established  from  experimentation  while  others 
have  beer,  constructed  from  the  theories  of  thermodynamics  and  fluid  mechanics. 
Some  models  predict  the  behavior  of  the  whole  body  while  others  are  specialized 


to  a  particular  body  segment.  The  Interest  In  this  research  lies  In 
models  of  the  entire  human  body.  Excellent  reviews  of  the  work  In  this 
field  are  available  by  Hardy  (1972),  Mitchell,  et  al .  (1972),  Shltzer 
(1972),  Fan,  et  al.  (1971),  and  Hwang  and  Konz  (1977). 

Mathematical  models  of  the  whole  human  body  may  be  generally  classified 
as  single  cylinder  or  multi-segment  models.  Further  classification  Is 
accomplished  by  determining  whether  thermoregulation  ts  Internal  or  external 
to  the  model.  Models  with  internal  thermoregulation  functions  may  be 
viewed  as  a  composite  of  two  submodels,  one  for  the  passive  system  (physical 
system)  and  one  for  the  controlling  system. 

The  complexity  of  models  of  the  passive  systems  depends  upon  the 
number  of  body  segments  modeled,  their  geometries,  the  number  of  nodes  and 
shells  (layers  of  segment  composition)  attributed  to  each  segment  and  the 
sophistication  of  the  model  circulatory  system.  The  complexity  of  con¬ 
troller  models  range  from  simple  function-evaluation  types  to  those  which 
compute  error  signals  based  on  variables  such  as  average  skin  temperature, 
core  temperature  and  skin  heat  flux.  All  controllers  determine  metabolic 
ate  and  In  some  models,  the  control  system  also  determines  the  sudometer 
(sweating)  and  vasometer  (variable  blood  flow)  responses. 

The  models  are  most  generally  expressed  as  a  set  of  differential  equa¬ 
tions.  Early  models  were  solved  using  analog  computers.  The  advent  of 
larger,  faster  digital  machines  has  resulted  In  most  models  now  being 
programmed  for  digital  computers.  The  solution  methodology  most  often 
employed  Is  that  of  finite  differences. 

1 . 2  Object  1 ves 

The  objective  of  this  research  was  to  determine  the  feasibility  of 
using  an  existing  human  thermal  model  In  the  evaluation  of  Immersion 
hypothermia  protection  devices.  To  accomplish  this  task  a  review  of  the 
literature  was  conducted  leading  to  the  selection  of  candidate  models 
representing  the  general  types  available.  Computer  codes  for  these  models 
were  obtained  and  modified  as  necessary  to  Implement  them  on  the  computer 
system  at  Clemson  University.  These  modifications  of  the  computer-based 
models  were  not  Intended  to  accomplish  basic  structural  changes  to  the 
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models  or  to  improve  their  intrinsic  capabilities. 

The  selected  candidate  models  were  to  be  vei (dated  against  the  data 
collected  from  the  human  Immersions  discussed  In  Part  I  of  this  report, 
see  Harnett,  at  al .  (1979).  Based  on  the  results  of  these  validations, 
recommendations  were  formulated  regarding  the  potential  usefulness  of  each 
selected  model.  Also,  the  most  Important  areas  for  model  Improvement  were 
determined  and  recommendations  regarding  their  priorities  were  developed 
for  consideration  prior  to  undertaking  any  efforts  aimed  at  improving  the 
capabilities  of  a  model. 

1.3  Scope 

The  scope  of  this  effort  included  only  an  evaluation  of  existing  models. 
Ho  new  model  development  or  modification  of  existing  models  was  required. 
However,  in  the  course  of  implementing  the  computer-based  models  and 
experimenting  with  them,  it  was  also  possible  within  project  cost  and 
schedule  constraints  to  develop  and  Implement  a  number  of  substantive  model 
modifications.  These  modifications  Improved  the  performance  of  the  affected 
moduls  and  in  some  cases  were  necessary  In  order  to  give  the  models  any 
chance  of  performing  as  required  for  the  evaluation  of  cold  water  protection 
equipment. 
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2.0  MODEL  SELECTION 


After  a  review  of  the  literature  It  was  decided  that  this  Investigation 
would  be  limited  to  non-steady  state  models  containing  controllers.  This 
decision  was  predicated  on  two  facts:  (l)  the  human  response  to  cold  water 
Immersion  was  unlikely  to  be  steady  state  and  (2)  we  would  be  unable  to 
externally  control  the  model  since  no  experimental  Immersion  data  was 
available  which  correlated  metabolic  rate,  vaso-constrlctlon  and  surface 
blood  flow  to  the  physical  parameters  of  the  modal  such  as  core,  skin 
and  water  temperatures. 

References  were  found  In  the  literature  to  but  a  few  models  fitting 
these  criteria.  Conversations  and  correspondence  with  the  authors  of  these 
models  revealed  several  other  unpublished  models.  However,  these  models 
were  all  reported  to  be  very  similar  to  those  found  through  the  literature 
search.  The  following  five  models  were  obtained:  Stolwljk,  (1972); 

Montgomery,  (1972);  Gordon,  (1972);  Kuznet ,  (197*0;  and  Wlnton  and  Llnebarger, 
(1971).  An  additional  model,  without  a  controller,  Wlssler  ( I S66)  was  re¬ 
quested  from  Its  author';  but  this  request  was  not  granted. 

These  models  represented  several  different  philosophies  In  the  modeling 
of  both  the  passive  and  control  systems.  Stolwljk's  and  Kuznet's  models 
had  been  used  primarily  for  Investigation  of  hyperthermia.  The  applications 
were  in  support  of  the  NASA  manned  space  flight  program.  Both  models  were 
implemented  on  Clemson's  IBM  370/165  and  exercised  with  test  data  as  specified 
by  their  authors.  Neither  model,  however,  was  modified  for  further  evaluation, 
Stolwljk's  model  was  dropped  from  further  Investigation  because  Montgomery's 
model  was  determined  to  be  an  extension  of  It  which  had  been  used  for  diving 
studies.  Kuznet's  model  was  dropped  because  of  !*s  similarity  to  both  the 
Gordon  and  Montgomery  models.  The  three  remaining  models,  Wlnton's, 
Montgomery's  and  Gordon's  were  chosen  because  of  the  contrasts  among  them 
in  terms  of  their  general  approaches  ar.d  complexity. 

The  Wlnton  and  Llnebarger  Model 

The  Wlnton  and  Llnebarger  (1971)  model  Is  the  simplest,  of  the  models 
studied.  This  model  was  Intended  for  study  of  both  the  steady-state  and 
transient  response  of  the  human  thermoregulation  system  to  various  degrees 
of  internal  and  external  thermal  stress.  Emphasis  In  this  model  was  placed 


on  the  feedback  structure  and  controller  mechanisms  Involved  In  thermo¬ 
regulation.  This  model  has  been  exercised  using  both  analog  and  digital 
simulation. 

Wlnton  and  LInebarger  represent  the  shape  of  the  body  as  a  cylinder 

having  three  concentric  layers.  The  Inner-most  layer  represents  the  core 

of  the  body,  which  Is  composed  of  the  deep  tissues  and  Internal  organs. 

Surrounding  the  Inner  core  Is  a  middle  ’ayer  made  up  of  muscle  and  fat. 

♦ 

The  skin  comprises  the  outer  layer  of  the  cylinder.  For  purposes  of 
analysis  the  thermal  properties  of  the  three-layer  model  are  represented 
by  an  analogous  electrical  circuit.  The  thermal  system  and  Its  electrical 
analog  are  governed  by  Identical  differential  equations  which  form  the 
bas Is  of  this  model . 

Three  primary  control  mechanisms  are  Included  In  this  model:  sudometer, 
vasometer  and  metabolic.  These  control  mechanisms  are  Incorporated  Into 
the  model  by  varying  the  related  parameters  of  the  circuit  analog  In  an 
appropriate  manner.  Physiological  studies  have  shown  the  Importance  of 
both  core  and  skin  thumperatures  In  thermoregulation.  Based  on  these 
studies  feedback  signals  from  the  core  and  skin  have  been  Included  In  this 
model . 

The  Montgomery  Model 

The  Montgomery  model  Is  an  extension  of  the  Stolwljk  (1970)  model 
Intended  to  allow  Investigation  of  heat  ioss/gain  during  underwater  diving 
work.  Its  thermoregulatory  system  Is  divided  Into  two  distinct  subsystems: 
the  phys ical -control  1 ed  subsystem  and  the  dynamic-controlling  subsystem. 

The  controlled  subsystem  contests  of  the  physical  portions  of  the  body. 

The  controlling  subsystem  contains  the  central  hypothalamic  thermo¬ 
integrator,  the  central  set  point  temperature  and  associated  afferent 
and  efferent  signal  pathways.  The  controlling  subsystem  receives  afferent 
signals  from  all  portions  of  the  body,  Integrates  the  signals,  compares  the 
results  to  the  central  set  point  and  distributes  the  appropriate  effector 
command  signals  to  all  portions  of  the  body. 

The  controlled  subsystem,  Figure  I  I  I  -  I  ,cons I st s  of  the  head  which  Is 
considered  a  sphere  and  cylinders  representing  the  trunk,  arms,  hands,  legs  and 
feet.  Both  arms,  hands, legs  and  feet  are  represented .because  of  symmetry, by  one 
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FIGURE  I  I  1-1 
SCHEMATIC  OF  THE 
MONTGOMERY  CONTROLLED  SYSTEM 


ARMS  HANDS 


FIGURE  111-2 
SCHEMATIC  OF  THE 
GORDON  CONTROLLED  SYSTEM 


SECTION  A-A 


Body  elements  that  comprise  the  passive  system 
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cylinder  each.  Each  body  segment  Is  composed  of  eleven  concentric 
layers:  four  representing  the  body  core,  four  representing  the  muscle 
layer,  and  one  each  representing  the  body  layers  of  fat  and  skin,  and 
one  representing  the  outer  wet  suit  layer.  A  central  blood  compartment 
simulates  the  body's  central  blood  pool  which  exchanges  heat  with  all  other 
body  compartments,  via  convective  heat  transfer,  through  simulated  blood 
flow  to  each  compartment.  Each  of  slxty-one  body  compartments  Is  represented 
by  a  heat  balance  equatlc/n  whtch  accounts  for  Internal  heat  generation, 
conductive  heat  transfer  between  adjacent  compartments,  and  convective 
heat  exchange  with  the  central  blood  compartment.  Where  applicable, 
respiratory  heat  generation  and  heat  loss  are  Included  In  the  body  compart¬ 
ment  heat  equations.  Additional  heat  balance  equations  represent  a  wet 
suit  on  the  six  body  segments  and  Include  the  effects  of  conductive  heat 
transfer  with  the  skin  and  conduction-convection  with  the  ambient  water. 

Each  of  the  sixty-seven  heat  balance  equations  Includes  the  thermal  capaci¬ 
tance  of  the  compartment  enabling  the  transient  response  of  the  compartment 
to  be  simulated. 

The  Gordon  Model 

The  Gordon  model  Is  the  most  complex  of  the  three  models  examined. 

Gordon  also  divides  the  temperature  regulatory  system  Into  two  major  subsystems: 
the  passive  subsystem  and  the  control  subsystem.  The  passive  subsystem  1$ 
divided  Into  ten  body  segments.  Cylinders  are  used  to  represent  the  neck, 
thorax,  abdomen,  arms,  hands,  legs  and  feet.  A  cylindrical  segment  represents 
the  face  and  spherical  segments  represent  the  head  and  forehead.  All  seg¬ 
ments  consist  of  four  concentric  layers.  Both  the  number  of  segments  and 
the  number  of  Integration  nodes  per  layer  Is  variable.  The  model  was 
Implemented  at  Clemson  with  the  segmentation  described  above  and  eleven 
nodes  per  segment.  Both  of  these  choices  were  used  by  Gordon  In  his  Initial 
validation  of  the  model.  They  represent  a  logical  balance  between  computational 
burd*n  and  numerical  precision. 

Since  the  body  is  modeled  as  concentric  spherical  or  cylindrical  shells, 
(Figure  I  I  1-2)  the  governing  equations  were  obtained  by  considering  a  s he  1 1  of 
uniform  properties.  An  energy  balance  equation  Is  used  for  each  shell.  To 
complete  the  passive  subsystem  model,  heat  transfer  Involving  the  blood  pool  was 
modeled,  Including  counter-current  heat  exchange  between  certain  body  elements 
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and  the  blood  pool.  The  partial  differential  equations  (which  from  the 
energy  balancing)  are  nonlinear  and  are  solved  by  finite  difference 
techn I ques. 

The  control  subsystem  In  this  model  generates  several  "error  signals" 
which  are  deviations  of  average  skin  heat  flux,  average  skin  temperature, 
and  hypothalmus  temperature  from  set-point  values  obtained  passively  during 
basal  conditions  (resting  with  no  food  In  stomach).  Combinations  of  these 

I 

error  signals  are  weighted  to  produce  the  control  subsystem  signal.  The 
equations  which  describe  this  control  signal  make  up  the  control  system 
model.  Although  the  original  emphasis  of  this  model  was  to  simulate 
exposure  to  cold  air,  provisions  were  made  for  the  future  addition  of  a 
warm-environment  controller.  The  computer  code  for  the  Gordon  model  was 
written  in  FORTRAN. 
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3.0  THERMAL  TESTING  PROCEDURE  AND  RESULTS 

Any  mode!  to  be  used  for  evaluations  of  a  cold-prctect ton  device 
must  be  provided  a  description  of  the  thermal  properties  of  the  device. 

The  models  considered  In  this  stud/  were  modified  to  accept,,  as  Input 
data*  the  thermal  conductivity  (resistivity)  of  the  devices.  Since  this 
data  was  not  available  for  many  of  the  devices  Included  In  this  study,  It 
was  necessary  to  determine  them  through  experimentation.  This  chapter 
presents  a  brief  description  of  the  testing  methodology  and  the  results 
obtained  with  It,  A  detailed  description  of  the  methodology  was  given 
by  Baker  (1979). 

Briefly,  the  procedure  was  to  monitor  the  temperature  In  a  bag  r f 
water  (Initially  warm)  which  had  been  surrounded  by  the  test  material, 
while  It  was  vigorously  agitated  In  a  large  tank  of  colder  water.  From 
this  data,  the  BTU  loss  could  be  calculated.  Assuming  linearity  over  a 
brlaf  span  (3~5  minutes)  an  estimate  of  the  thermal  resistance  (R)  of  the 
material  was  calculated.  Values  obtained  In  this  way  are  only  approxi¬ 
mations  but  were  seen  to  compare  well  with  published  values  obtained  using 
tne  methodology  of  ASTM-C15I8-70. 

Each  device  Included  In  the  cold- immers Ion  test  (except  the  PFD)  was 
tested  to  determine  Its  thermal  resistance,  A  number  of  tests  were  per¬ 
formed  on  each  device.  Any  test  whose  result  was  suspect  due  to  obvious 
nonlinearity  was  discarded.  This  resulted  In  sample  sizes  ranging  from 
three  to  five  for  Individual  protective  devices.  Table  ISI-I  presents 
therm* I  resistance  values  (R  values)  determined  In  this  testing  and  co¬ 
efficient  of  heat  transfer  values  (U  values). 

The  suits  which  relied  upon  Aramld  underwear  to  provide  thermal  protection 
were  tasted  with  a  swatch  of  It  placed  between  the  bag  and  the  device  being 
tested.  These  suits  are  Identified  In  Table  MI-1  with  the  comment  (W/A) . 

St  w®s  noted  that  In  the  case  of  the  CWU-2I/AP,  the  Artmld  swatch  was  wet 
at  the  completion  of  each  test.  A*  was  pointed  out  In  Part  I  of  this  report, 
Harnett,  et  al.  (1979),  this  suit  was  also  observed  to  leak  during  Its  In 
vivo  cold  Immersion  testing. 

The  1 LC  Industries  prototype  end  Dr.  Rentsch's  prototype  rely  mainly 
on  an  InfUtatla  elr  bladder  to  provide  thermal  protection.  The  conductivity 
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TABLE  ! 1 1- I 

RESULTS  OF  THERMAL  TESTING 


st  number  is  the  vaiue  determined  experimentally  for  the  deflated  suit, 
second  number  is  the  calculated  vaiue  for  thesuit  inflated, 
standard  deviation  is  shown  for  the  calculated  value. 


of  these  two  devices  could  noc  be  determined  by  direct  tests  as  previously 
described.  Instead,  we  measured  the  Intrinsic  conductivity  of  the  suits 
deflated.  These  values  along  with  the  average  thickness  of  the  air  space 
with  the  suit  Inflated  were  used  to  calculate  an  overall  conductivity. 

The  calculated  values  shown  In  Tabla  lll-l  are  actually  the  average  of 
two  calculated  values.  One  value  was  calculated  assuming  that  the  air 
space  was  devoid  of  convection;  the  second  assumed  that  the  convection- 
conduction  was  equivalent  to  that  between  two  smooth  parallel  plates. 

The  first  assumption  leads  to  an  optimistic  value,  the  second  to  a  pessi¬ 
mistic  value.  The  average  was  taken  with  uniform  weighting.  This  approx! 
matlon  Is  necessitated  by  the  geometry  of  the  problem  which  does  not  lend 
Itself  to  solution  by  known  techniques. 

The  accuracy  of  this  method  Is  unknown.  It  Is  suggested  that  In  the 
future,  conductivity  tests  be  run  using  a  heat  balance  system  similar  to 
that  recommended  In  ASTM-C5I8-70  but  adapted  to  testing  In  "the  wet". 

For  example,  a  brass  sphere  containing  water  and  a  heating  coll  could  ce 
covered  with  the  material  to  be  tested.  By  monitoring  the  temperature  In 
the  vessel  during  Immersion,  the  power  through  the  coll  can  be  adjusted 
until  a  steady-state  temperature  Is  achieved  within  the  vessel.  The  power 
being  supplied  to  the  coll  Is  then  directly  related  to  the  conductivity 
of  the  protective  material. 


4,0  EVALUATION  OF  THE  SELECTED  MODELS 

The  models  selected  were  modified  to  accept  the  Inputs  from  the 
thermal  testing  described  In  the  last  chapter.  This  simplified  the 
process  of  exercising  the  models  for  the  variety  of  test  articles  con¬ 
sidered  In  the  study.  It  was  also  necessary  to  modify  by  adding  a  layer 
to  represent  the  protection  equipment.  Further,  It  was  necessary  to 
modify  the  modeling  of  convective-conductive  heat  transfer  av  the  Interface 
with  the  environment  to  reflect  the  differences  resulting  from  water  as 
opposed  to  air  Immersion. 

4.1  The  Wlnton  and  Linebarger  Model 

The  Wlnton  and  Linebarger  model  Is  the  simplest  of  the  three  selected  for 
evaluation.  The  three-layered  cylinder  used  to  model  the  body  was  represented 
by  the  electrlcal  circuit  analog  shown  In  Figure  I  1 1  -  3 • 

FIGURE  1 i 1-3 

WINTON-l I NE BARGER  3" LAVER  MODEL 


The  model,  originally  Implemented  on  an  analog  device  was  later  written  In 
CSMP.  A  fSMP  version  of  the  model  was  Implemented  and  evaluated  In  this 
study. 

In  Figure  I  I  I - 3 *  H£  Is  a  current  source  which  represents  basal  metabolic 
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thermogenesls,  shivering  thermogenes  Is  and  respiratory  heat  loss.  G  Is 
a  resistor  representing  the  net  heat  conduction  of  the  body  tissues  and  G 
represents  the  net  conduction  from  body  to  environment.  C  Is  a  capacitor 
representing  the  heat  capacitance  of  the  body  tissues  except  for  the  skin 
which  Is  represented  by  C#.  H,  Is  a  current  source  representing  evaporative 
heat  loss,  T#  I*  a  voltage  source  representing  ambient  environmental 
temperature  and  T.  and  T#  represent  core  and  skin  temperatures,  respectively 

To  accommodate  modeling  of  a  thermal  protective  device  It  was  necessary 

to  modify  the  circuit  In  Figure  lll-3.  Figure  111-4  shows  the  modified 
model . 

FIGURE  Ml 

MODIFIED  Wl NTON-U NE BARGER  MODEL 


Htre  Gw  represents  the  conductance  of  the  protective  device,  and  C  Its 
thermal  capacitance.  The  current  source,  «w,  was  Included  to  allow  the  modelfnc 
o  evlces  which,  themselves,  produce  heat.  Is  the  temperature  of  the  pro-  ’ 
tectfve  device.  For  Immersion  simulation  He  Is  set  for  zero.  AM  of  the 
protective  devices  considered  In  this  study  are  passive,  so  H  was  set  to  zero. 

t,n9  °f  the  CSHP  Program  used  to  Implement  the  model  In  Figure  IN-4  Is 
shown  In  Appendix  A. 

The  rode  I  performed  well  considering  Its  simplicity.  Core  t.mper.tor. 
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losses  were  of  the  same  magnitude  as  those  observed  during  a  human 
Immersion.  The  time  trace  Itself  however  was  not  the  same  as  Is  seen  In 
Figure  S 1 1  -  5*  The  model  loses  a  significant  amount  of  heat  before  the 
metabolic  process  can  begin  to  compensate.  The  metabolic  rate  Is  driven 
as  high  as  necessary  to  overcome  the  rate  of  heat  loss  and  the  model 
finally  exhibits  a  plateaulng. 

FIGURE  III -5 
TIME-TEMPERATURE  PROFILE 
FROM  WINTON'S  MODEL 


TIME  FROM  IMMERSION  (min) 


The  Initial  large  loss  of  heat,  as  reflected  by  a  rapidly  declining 
core  temperature,  Is  attributed  In  the  ma*n  to  the  controllers  Inability 
to  react  to  the  high  gradient  experienced  upon  Immersion.  As  with  many 
models,  Wlnton's  was  developed  for  use  In  estimating  steady-state  conditions. 
As  a  consequence,  large  continuing  changes  are  required  In  core  and  skin 
temperatures  before  the  metabolic  rate  becomes  sufficiently  elevated  to 
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overcome  the  rate  of  heat  loss  due  to  the  Immersion. 

Making  the  mode!  more  useful  for  evaluating  protective  devices  vould 
require  modification  of  the  controlling  subsystem.  Because  the  mode 
represents  the  surface  of  the  subject  as  a  single  cylinder,  It  Is  applicable 
oniy  to  the  study  of  full  body  protective  devices.  Some  general  I zat  ons 
to  this  model  are  possible.  With  modification  of  the  controller  goo<  estimates 
of  cooling  rate  may  be  possible  for  full  body  suits.  Since  the  model  requires 
small  amounts  of  both  computer  time  and  storage,  there  may  be  sufficient 
Justification  to  warrant  this  modification  effort. 

4.2  The  Gordon  Model  * 

As  described  In  Chapter  2,  Gordon  represented  the  body  by  a  set  of 
spherical  and  cylindrical  segments  consisting  of  concentric  layers.  Each 
layer  was  subdivided  to  finally  represent  each  segment  as  a  set  of  partial 
differential  equations  across  eleven  nodes.  Each  node  was  centered  In  a 
shell  of  uniform  material.  The  resulting  equations  are  solved  using  finite 
difference  techniques. 

The  addition  of  a  protective  device  to  this  model  required  only  two 
small  changes  to  the  model.  An  additional  layer  was  modeled  for  each  body 
segment  covered  by  the  device.  Each  segment  originally  had  four  layers. 

Because  of  the  rapid  transient  temperatures  to  be  experienced  In  this  layer, 

It  should  be  modeled  with  not  less  than  two  Integration  nodes, Wi ssler  (1971). 
Thus,  at  a  minimum,  the  model  would  consist  of  thirteen  differential  equations 
for  each  protected  segment. 

The  second  change  required  was  In  the  control  subsystem.  In  the 
original  model,  the  controller  Is  hooked  to  the  last  Integration  node,  that 
of  the  skin.  Since  the  model  was  built  to  accept  a  varying  number  of  nodes 
per  segment,  the  addition  of  the  protective  device  would  cause  the  controller 
to  use  as  Input  temperatures  and  fluxes  calculated  for  the  furthest  node 
from  the  segment  core  which  would  be  protective  device  surface  temperatures 
and  fluxes.  This  Is  easily  remedied  In  either  of  two  manners.  First,  one 
may  fix  the  number  of  nodes  representing  body  In  each  segment  thus  affixing 
to  the  skin  a  constant  node  number.  Alternatively  the  node  number  representing 
the  skin  may  be  calculated  by  Including  as  Input  the  number  of  nodes  used 
to  represent  the  protective  device  on  each  segment. 
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As  mentioned  above,  the  three  models  to  be  modified  were  first  examined 
under  conditions  for  which  they  had  been  designed.  In  this  testing  Gordon's 
model  was  found  to  be  the  best  overall  analog.  Before  modifications  were 
made  to  model  the  protective  device,  the  model  was  set  up  to  simulate  a 
nude  cold  Immersion.  As  a  result  of  these  runs,  a  major  problem  was 
observed  In  the  Gordon  model.  One  would  expect  that,  after  an  Initial 
rise  in  core  temperature,  a  decrease  would  occur  which  may  or  may  not  find 
an  equilibrium  level.  Th‘e  period  of  Initial  rise  and  the  magnitude  of  the 
rise  can  vary.  However,  from  experimental  evidence  one  would  expect  the 
period  to  be  relatively  short  with  a  magnitude  of  at  most  a  few  tenths  of 
a  degree  celslus  for  a  nude  man  In  cold  water.  The  Gordon  model  exhibited 
the  initial  Increase  upon  Immersion.  However,  the  core  temperature  continued 
to  climb  as  the  simulated  Immersion  continued.  At  the  end  of  a  5-hour  simulated 
Immersion,  the  core  temperature  had  reached  approximately  42°C. 

A  re-examinat Ion  of  the  model  under  the  conditions  for  which  It  was 
validated  by  Its  author  showed  that  the  problem  could  have  existed  at  that 
time.  It  was  validated  during  relatively  short  (2  hour) s fmul at  Ions  of  ex¬ 
posure  to  cold  air  (5°C).  Increasing  core  temperatures  were  expected  during 
this  period.  Their  appearance  In  the  simulation  results  seemed  to  Indicate 
that  the  model  was  working  well.  Unfortunately,  the  validations  with  cold 
air  Immersions  did  not  extend  Into  the  time  period  when  core  cooling  occurs. 

This  problem  was  not  observed  In  the  Initial  testing  at  Clemson  because  the 
simulations  were  also  of  relatively  short  duration. 

The  Gordon  model  employs  the  most  detailed  model  log  of  the  human  physical 
structure  of  the  three  models  examined.  It  was  anticipated  that  It  would 
prove  to  be  the  most  accur.te  of  the  three  selected  models.  For  this  reason, 
while  It  was  not  an  objective  of  this  effort  to  Identify  end  correct  modeling 
errors,  an  attempt  to  do  so  was  made.  The  eductions  for  the  physical  system 
were  verified  by  reder I  vat  Ion.  The  coding  was  checked  and  several  simplifi¬ 
cations  were  made.  However,  the  problem  could  not  be  Identified.  It  was 

finally  decided  to  terminate  efforts  to  correct  the  Gordon  computer  code.  In 
favor  of  concentrating  on  the  Montgomery  model. 


4.3  The  Montgomery  Model 

As  suggested  earlier,  the  Montgomery  mode!  Is  basically  an  adaptation 
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of  Stolwijk's  model  designed  for  use  In  evaluating  the  thermal  aspects  of 
diving  activity.  The  model  divides  the  body  Into  ter  compartments/segments. 
Each  segment  is  subdivided  Into  eleven  concentric  layers.  Each  compartment 
represents  a  lumped  thermal  capacitance  with  appropriate  modes  of  heat 
production  and  heat  transfer  to  other  compartments.  Each  body  layer  gene¬ 
rates  metabolic  heat  at  a  basic  rate  and  exchanges  convective  heat  with 
the  central  blood  pool.  The  eleven  compartments  of  each  segment  exchange 
heat  via  conductive  trans'fer  with  adjacent  compartments  as  function  of 
layer  geometry  and  tissue  thermal  conductivity.  Each  wet  suit  segment 
exchanges  heat  with  the  environment  as  a  function  of  wet  suit  properties 
and  ambient  water  conditions. 

The  starting  point  In  the  development  of  the  thermal  network  for  a 
given  subject  Is  to  estimate  his  percent  body  fat  from  his  height  and 
weight.  His  total  surface  area  Is  also  estimated  from  his  helgitand 
weight.  The  surface  area  of  each  segment  Is  then  calculated  from  the  total 
surface  area. 

The  relative  weights  of  the  various  segmental  layers  are  calculated 
from  the  subject's  total  body  weight  and  a  percentage  weight  distribution. 

The  various  compartment  weights,  when  multiplied  by  the  corresponding 
specific  heat  value,  yields  the  thermal  capacitance  value  for  each  compart¬ 
ment.  Since  the  core  of  each  segment  consists  of  both  skeletal  and  visceral 
tissue  which  have  different  specific  heat  values,  they  must  be  treated 
separately  and  averaged. 

The  central  blood  compartment,  representing  the  blood  In  the  heart  and 
the  great  vessels,  is  assumed  to  contain  2.5  liters  of  blood.  The  thermal 
capacitance  of  the  central  blood  compartment  Is  subtracted  from  the  total 
thermal  capacitance  of  the  trunk  core.  The  metabolic  heat  generation  in 
each  body  compartment  Is  calculated  using  the  distribution  given  by  Stolwljk 
and  Hardy  (1 966). 

The  convective  heat  exchange  that  takes  place  between  each  body  compart¬ 
ment  and  the  central  blood  pool  as  a  result  of  blood  flow  Is  calculated 
using  the  basal  blood  flow  values  for  each  compartment.  The  thermal  con¬ 
ductances  of  each  compartment  are  assumed  to  be  uniform,  concentrated  at 
the  compartment's  center  of  mass  and  only  dependent  upon  compartment  tempera¬ 
ture.  thermal  conductance  between  layers  Is  a  function  of  the  thermal 
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conductivity  of  the  material  between  compartments,  the  distance  between 
compartments,  the  area  of  the  heat  transfer  surface  located  at  the  mid¬ 
plane  between  compartments  and  the  temperature  of  the  two  compartments. 

The  wet-sul t-to-amblent-water  heat  transfer  coefficients  are 
dependent  upon  the  geometric  shape  of  the  body  segment,  the  ambient 
temperature  and  pressure;  the  viscosity,  thermal  conductivity,  heat  capa¬ 
city  and  density  of  the  surrounding  water;  and  the  water  velocity  relative 
to  the  body  segment.  Wet  suit  compartment  thermal  capacitance  values  are 
calculated  as  a  function  of  wet  suit  specific  heat  and  wet  suit  density 
for  each  segment. 

The  Stolwljk  (1970)  blothermal  model  was  used  to  form  a  basis  for  the 
controlling  subsystem.  The  first  change  that  was  made  to  the  Stolwljk  model 
was  to  provide  more  compartments  to  represent  the  core  and  muscle  tissue 
of  each  body  segment.  Finite  difference  methods  of  solution  using  lumped 
nodes  will  produce  errors  when  new  gradients  develop  In  the  relatively 
thick  muscle  and  core  layers.  This  type  of  error  may  be  decreased  by 
Introducing  additional  compartments  In  the  core  and  muscle  portions  of 
tne  controlled  subsystem.  This  method  was  used  by  WIssler  ( 1 96^1 )  to  Improve 
the  simulated  response  to  cold  exposure. 

The  core  and  muscle  portions  of  each  segment  were  divided  Into  four 
compartments,  each  having  one-fourth  of  the  cor«  of  muscle  mass  of  the  given 
segment.  An  additional  compartment  was  also  provided  to  represent  the  wet 
suit  covering  each  body  segment. 

The  effect  of  evaporative  heat  loss  from  the  skin  compartments  Is 
negligible  under  diving  or  totally  Immersed  conditions.  The  evaporative 
heat  loss  from  the  trunk  core  Is  equal  to  that  amount  of  heat  that  Is 
carried  away  from  the  body  during  expiration  of  the  respiratory  gas.  The 
quantity  of  heat  loss  from  the  respiratory  tract  for  any  gas  mixture  can 
be  calculated  from  the  physical  properties  of  the  gas  mixture  and  the 
thermal  and  dynamic  character  1st les  of  the  respiratory  system.  Respiratory 
heat  loss  Is  proportional  to  the  respiratory  minute  volume,  which  Is  In 
turn  related  to  the  amount  of  oxygen  required  to  provide  energy  for 
metabolic  needs.  The  respiratory  loss  of  heat  Is  somewhat  offset  by  the 
work  of  breathing,  the  net  amount  of  evaporative  heat  loss  from  each  of 
the  trunk  core  compartments  Is  taken  to  be  one-fourth  of  the  difference 
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between  the  total  respiratory  heat  loss  and  the  total  heat  generation 
due  to  respiration. 

Net  heat  flow  Into  or  from  each  compartment  Is  then  calculated.  The 
skin  compartment  In  each  segment  loses,  through  conduction  to  the  wet 
suit,  an  amount  of  heat  equal  to  the  heat  transfer  coefficient  multiplied 
by  the  temperature  difference  between  the  skin  and  wet  suit  compartments. 

Since  a  diver  does  not  receive  solar  heat  input  and  does  not  transfer 
radiant  heat  to  his  surroundings,  the  environmental  heat  transfer  coefficients 
used  by  Stolwljk  (1970)  have  been  replaced  by  convective-conductive  heat 
transfer  coefficients  between  the  wet  suit  and  ambient  water.  The  water- 
neoprene  skin  surface  heat  transfer  coefficient  Is  dependent  upon  the 
geometrical  shape  of  the  body  segment;  the  ambient  temperature  and  pressure; 
the  viscosity,  thermal  conductivity,  heat  capacity  and  density  of  the 
surrounding  water;  and  the  water  velocity  relative  to  the  body  segment. 

Little  modification  was  necessary  to  the  physical  system  of  Montgomery's 
model  since  It  already  had  provisions  for  an  additional  layer  between  the 
skin  and  environment.  Initial  modifications  were  directed  toward  increasing 
the  flexibility  of  device  modeling  to  allow  simulation  of  other  than  full 
body  suits.  In  addition,  modifications  suggested  by  Montgomery  (personal 
communications)  were  implemented  to  simulate  Immersion  of  the  body  with 
the  head  and  neck  (modeled  as  a  single  segment)  exposed  to  air.  These 
changes  were  made  In  the  WETMAN  subroutine  of  the  model  computer  code. 

Typical  Results  with  Montgomery's  Model 

Runs  were  made  to  simulate  several  of  the  devices  Included  In  the  cold 
Immersion  test  described  in  Part  I  of  this  report,  Harnett  et  al .  (1979)- 
The  results  of  one  series  of  these  simulations  Is  shown  In  Figure  I  I  1-6. 

The  article  modeled  (WP3)  was  a  Jacket-type  device  providing  protection 
basically  to  the  trunk  and  arms.  The  figure  presents  the  simulation  results, 
for  a  man  172  cm  tall  weighing  Jk.k  kg,  and  the  experimental  observations 
for  five  volunteers  testing  this  device. 

Considerable  variation  In  response  Is  exhibited  by  the  five  experimental 
observations,  largely  due  to  somatotype  differences  among  the  subjects.  None 
of  the  individual  responses  Is  represented  well  by  the  model  results.  The 
model's  metabolic  control  subsystem  Initially  fails  to  recognize  the  heat 
drain  caused  by  the  cold  immersion.  This  results  In  a  too  slow  Increase 
In  the  rate  of  thermogenesis  and  a  rapid  cooling  shortly  after  Immersion. 
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RESULTS  WITH  MONTGOMERY'S  MODEL  AND 
FIVE  COLD  IMMERSIONS 


. 


FROM  IMMERS 


When  the  metabolic  controller  does  recognize  the  loss  It  elevates  the 
metabolic  rate  by  large  steps  to  compensate  for  the  decline  In  average 
skin  and  head  core  temperature.  The  result  Is  a  complete  stabilization  of 
core  temperature. 

Montgomery's  Metabolic  Rate  Controller 

The  controller  for  metabolic  rate  formulated  by  Stolwijk  and  used  by 
Montgomery  calculates  changes  In  metabolic  rate  based  upon  the  following 
formulat Ion. 


AMR 


where  AMR 


C  •  (AT  +  R  •  T'  )  +  C  •  (  AT  +  R  •  T')  + 

c  c  cc  s  s  ss 

C  •  (AT  +  R  -T'  )  •  (AT  +  R  •  T'  ) 

C  CCS  s  s 

change  In  metabolic  rate 

C  ,  C  ,  C,  R  and  R  are  weighting  factors 
c  s  c  s 

AT  deviation  of  head  core  temperature  from  a 

"  set  point  value 

T'  ■  rate  of  change  of  head  core  temperature 

AT^  "  deviation  of  average  skin  temperature  from  a 
set  point  value 


T'  -  rate  of  change  of  average  skin  temperature 


The  weighting  factors  (constants)  are  defined  as  follows. 

R  -  0 

c 

Rs  -  0.03 

C  -  0 

c 

Cs  -  0 

C  -  21.0 

The  definitions  of  C^,  Cs  and  C  ware  used  by  Stolwijk  based  upon  the  experi¬ 
mental  evidence  of  Benzlnger,  et  al.  (1963) ■  With  these  definitions  the 
metabolic  rate  controller  simplifies  to  the  following. 

AMR  -  21  • AT  • (AT  +  0.03  T '  ) 

c  s  s 

*Average  “sSTTT T  temperature  Is  defined  as  the  average  of  the  outer  surface  tempera 
of  each  compartment,  weighted  by  their  proportional  amounts  of  surface. 
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This  simplified  model  Is  completely  Insensitive  to  the  rate  of 
head  core  cooling  and  Is  not  particularly  sensitive  to  the  rate  of  skin 
cooling.  The  mode!  Is  largely  based  on  the  amounts  of  cooling  In  these 
two  temperatures  and  places  essentially  equal  emphasis  on  the  two. 
Initially,  there  will  be  a  rapid  decrease  In  average  skin  temperature 
but  very  little  change  In  core  temperature.  This  Implies  that  the  product 
In  the  above  equation  wilt  be  small,  resulting  In  small  Increases  In  the 
metabolic  rate,  until  thdre  has  bean  a  significant  decrease  In  head  core 
temperature. 

This  helps  explain  the  behavior  of  the  model  as  depicted  In  Figure 
111-6.  It  is  obvious  then  that  there  is  considerable  room  for  Improvement 
In  modeling  metabolic  control.  While  not  Included  In  the  scope  of  this 
project,  the  development  of  a  new  controller  was  seen  to  be  essentia)  to 
a  positive  finding  that  a  model  can  serve  the  purpose  addressed  In  this 
study. 

The  validity  of  Montgomery's  (Stolwljk's)  metabolic  controller  for 
certain  conditions  (e.g.,  a  nude  man  In  cold  air)  has  been  shown  by  many 
people  Including  Stolwljk.  It  was,  as  stated  above,  based  on  experimental 
data.  It  was  decided,  therefore,  to  formulate  a  new  controller  for  the 
Immersion  environment  based  on  the  data  collected  In  the  human  Immersion 
portion  of  this  study  summarized  In  Part  I  of  this  report. 


Improved  Metabolic  Controllers 


Relevant  data  available  Included  rectal  temperature  at  15  cm,  skin 
temperatures  at  the  toe,  thigh,  forearm,  bicep,  groin  and  subscapular  sites 
and  periodic  measurements  of  metabolic  rate.  The  procedure  used  to  develop 
the  controller  was  to  establ Ish, through  regression  analysis,  linear  models 
relating  changes  In  metabolic  rate  to  changes  In  the  rectal  ar.d  skin  tempera¬ 
tures.  The  basis  of  comparison  for  determining  these  changes  were  measurements 
made  following  a  30-mlnute  rest  period  prior  to  commencing  cold  Immersion. 


The  Initial  attempt  was  based  upon  regression  analysis  applied  to  the 
pooled  sample  of  observations  obtained  In  the  laboratory.  When  the  resulting 
metabolic  controller  was  Implemented  In  Montgomery's  model,  fair  predictions 
of  cooling  rates  resulted  for  most  of  the  wet-mode  suits  but  the  predictions 
for  the  abandon-ship  type  dry  suits  were  much  worse. 


The  experimental  data  was  then  segregated  In  two  subsets  --  one  cbtslned 
during  Immersion  In  ll.8°C  water  with  wet-mode  suits  and  one  obtained  during 
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*5  ,  .»$j\\\!mv7>i.va  wanwy^sy,  A<j-w  wgpy^»:- 


Immer-lon  In  I  - 7* C  water  with  dry-mode  suits.  These  data  sets  warn 
analyzed  separately  to  produce  a  metabolic  controller  for  each  condition, 
for  the  wet-mode  suits  the  following  model  was  obtained. 

AMR  -  7.00**  -  5.822  •  AT  .  -  2.1*07  •  AT..  -  37.382  •  AT 

tn  r  r 

where  AT  ^  -  change  tn  thigh  temperature 

AT ^  .•  change  In  forearm  temperature 

AT^  «  change  In  rectal  temperature 

This  model  had  a  correlation  coefficient  (R-  value)  of  0.77.  The 
accuracy  with  which  regression  relations  conform  to  the  data  Is  often 
expressed  by  "F  test  statistics".  The  significance  of  these  statistics 
may  be  Interpreted  by  the  "leva!  of  significance"  at  which  the  hypothesis 
(that  the  observations  follow  the  model)  me''  be  rejected.  Small  levels  of 
significance  Indicate  that  the  regression  conforms  well  to  the  observations. 
The  "level  of  significance"  for  this  regression  was  0.001. 

The  control  model  obtained  from  regression  analysis  of  the  data  obtained 
with  the  dry-mode  suits  Is  the  following. 

AMR  -  1.987  -  2.65*  -AT  -  1*. 595  •  AT  .  - 

t  th 

5.361  •  AT,  -  8.920  -  AT 
ft  r 

where  ATt  ■  change  in  toe  temperature 

AT^  ■  change  In  bicep  temperature 
other  symbols  as  previously  defined 

This  model  had  a  correlation  coefficient  of  .85  and  a  "level  of  significance" 
of  .0001 . 

These  two  relationships  were  implemented  In  the  Montgomery  model. 

It  was  necessary  to  accept  some  approxlmat Ions  In  marrying  the  list  of 
variables  requlrarl  by  the  controllers  with  those  available  In  Montgomery's 
model .  The  variables  were  matched  as  follows. 


Con t ro 1  1  er  Va r ! afal e 


Mode!  Variable 


Toe  temperature  Average  *ool  temperature 

Thigh  temperature  Average  leg  temperature 

Bleep  temperature  Average  arm  temperature 

Recta!  temperature  Trunk  core  temperature 

These  approximations  are  unavoidable  because  of  the  simplifications  In¬ 
volved  In  modeling  the  physical  structure  of  the  body  In  Montgomery's 
model . 

Model  Val Idat Ion 

The  question  of  model  validity  Is  dependent  upon  the  use  to  which 
tne  model  Is  to  be  put.  The  objective  In  this  case  Is  to  use  the  model 
rather  than  human  experimentation  as  the  basis  to  estimate  the  survival 
time  associated  with  new  developments  In  protection  equipment.  If  one 
accepts  the  survival  time  mode)  and  prediction  procedure  presented  In 
Part  !,  then  all  that  Is  required  of  the  model  Is  a  prediction  of  core 
cooling  rate  which  may  than  be  used  to  estlmato  survival  time.  This  would 
relieve  the  need  to  be  particularly  concerned  with  absolute  temperatures 
predicted  By  th«  model  for  various  body  sites  or  transient  aspects  of 
their  profiles. 

Based  on  this  method  for  estimating  survival  time,  model  validity 
may  be  determined  by  establishing  the  accuracy  of  Its  predictions  of 
the  rate  of  core  cooling.  This  may  be  done  by  performing  statistical 
tests  comparing  the  rates  observed  with  the  volunteer  test  subjects 
(Part  I  of  this  report)  to  corresponding  rates  predicted  for  them  by  the 
model.  This  data  Is  naturally  "paired"  and  so  lends  Itsalf  to  paired 
analysis  as  a  means  of  variance  reduction.  The  "paired  t  test"  described 
by  Steel  and  Forrle  (I960)  was  used  for  this  purpose  The  procedure  Is 
Illustrated  below.  The  test  was  run  at  the  0.05  level  of  significance  with 
a  two-tailed  rejection  region. 

Null  Hypothesis  1H  );  There  Is  no  difference  between  mean 
0  simulated  and  mean  experimentally 
observed  cooling  rater. 


2k 


Subject 

BHHHBB 

mn 

Deviation  (d) 

Simulation 

HTTTSliiHiW  ltfcJI 

Simulation  -Experimental 

MK 

1.037 

2.510 

-1.473 

GE 

1 .1472 

1.800 

-0.328 

TP 

1 . 100 

0.556 

0.544 

BH 

.697 

0.391 

0.30 

PK 

.621 

0.244 

0.376 

$d2  »  2.S08  Id  «  0.575  n  ■  5  (No.  of  pairs) 
d  -  d/n  -  0. 1 15 

S2  *  Id2  -*rT(id)  “  0.137 
d  n(n-T) 

S5  -  0.367 

t-test  statistic  ■  d/S^  ■  0.310 

Critical  value  of  t  at  0.05  level  of  significance  with 
4  degrees  of  freedom  *■  2,776 

Since  the  calculated  statistic  Is  less  than  the  critical 
value  we  cannot  reject  H  . 

A  summary  of  these  tests  for  each  of  the  devices  Included  In  the  cold 
immersion  testing,  except  for  the  Pt-'D  Is  presented  In  Table  MI-2.  From 
the  table  we  observe  that  the  model  performed  well  for  ail  of  the  wet-mode 
suits. 

Those  suits  for  which  we  must  reject  the  hypothesis  of  sameness  between 
model  and  observed  average  cooling  rates  are  largely  the  abandon-ship  suits. 
Reference  to  Figure  l-l  In  Part  I  of  this  report  will  help  explain  the  models 
failure.  The  model  "sees"  all  simulated  Immersions  with  the  subject  com¬ 
pletely  underwater  from  the  neck  down.  As  can  be  easily  observed  from  the 
pictures  of  the  flotation  attitudes  of  these  suits  In  Figure  l-l,  much  of 
the  upper  surface  areas  covering  the  legs,  arms  and  trunk  of  each  of  these 
suits  Is  exposed  to  air.  The  heat  transfer  coefficient  for  air  Is  much  less 
than  that  of  water.  One  would  therefore  expect  the  cooling  rates  predicted 
by  the  model  for  these  devices  to  be  greater  than  that  observed  In  the 
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RESULTS  OF  PAIRED  t  TEST 
WITH  REGRESSION-BASED  CONTROLLERS 


TEST  ARTICLES 

-  -  — --  -  -  -  -  -  - 

|  Calculated 
!  Statistic 

Crl tfcal 
Value  @ 

0.05  Level 

Degrees 

of 

Freedom 

Accept  or 
Reject 

Bay  ley  Exposure  Suit  (PVC  foam) 

5.367 

2.776 

1* 

Rej ect 

Bayley  WeatherMate  Plus 

1.599 

2.776 

4 

Accept 

HeMy-Hansen  Survival  Suit  (D-600-0) 

4.785 

2.571 

5 

Reject 

Henderson  Zlp-On  Exposure  Suit  (208C-4) 

5.237 

2.776 

i* 

Reject 

Henderson  Prototype  Jacket 

1 .  ISO 

2.776 

Reject 

ILC  Industries  Prototype  Survival  Suit 

0.387 

2.776 

1* 

Accept 

Medalist  Ski  Shorty  (7010) 

8.  172 

2.776 

1* 

Reject 

Mustang  U-ViC  Thermofloat  (1661) 

0.009 

2.776 

1* 

Accept 

NADC  Goretex  Experimental  Coverall 

9.  102 

2.571 

5 

Reject 

Dr.  3.  B.  Rentsch's  Prototype  Survival  Suit 
(without  respiratory  heat  reclamation) 

1 . 289 

2.776 

4 

Accept 

S.I.D.E.P.  Seastep  Survival  Suit 

15.931 

2.776 

4 

Rej  ect 

Stearns  Windjammer  Jacket  (FJ-55) 

0.1*16 

2.776 

4 

Accept 

Stearns  Offshore  Survival  Jacket  (FS-500) 

1.823 

mm 

5 

Accept 

Stearns  Heavy-Duty  Offshore  Survival 

Suit  (FS-71) 

6.068 

2.776 

4 

Reject 

U.5.  Air  Force  Modified  Ant  I -Exposure 

Assembly  (CWU-21/AP) 

2.  187 

2.776 

4 

Accept 

experimentation,  as  Is  the  case. 

The  Stearns  Heavy-Duty  Offshore  Survival  Suit  and  the  Henderson  Zlp-On 
Exposure  Suit  failed  to  pass  the  test.  Reference,  to  Figure  l-l  shows  the 
same  type  of  flotation  attitude  for  these  suits  as  for  the  abandon-ship 
suits.  In  addition  to  the  above-stated  reason,  the  inadequacy  of  the 
controller  may  be  contributory  to  their  failure. 

Model  results  for  two  additional  suits  failed  the  paired  t  test.  In 
both  cases  the  model -pred Icted  cooling  rate  was  faster  than  that  observed 
from  our  testing.  Again  reference  to  Figure  l-l  shows  that  both  of  these 
devices  were  tested  with  the  subjects  wearing  a  "water  wings"  type  flotation 
device  which  exposed  a  significant  amount  of  the  trunk  to  the  air.  Addltiona 
the  flotation  device  allowed  the  subjects  to  keep  their  arms  and  hands  out 
of  the  water.  Thus  a  faster  predicted  cooling  rate  from  the  model  seems 
very  reasonable. 

Three  suits  passed  the  test  which,  by  the  above  arguments,  should  have 
failed,  the  ILC  prototype,  Dr.  Rentsch's  prototype  and  the  U.  S.  Air  Force's 
CWU-2I/AP.  All  three  of  these  had  simulation  predicted  cooling  rates  slower 
than  we  expected.  For  the  IlC  and  Dr.  Rentsch's  prototypes  this  Is  most 
probably  due  to  error  In  the  estimation  of  their  thermal  conductivities. 

The  estimated  conductivities  are,  therefore,  probably  smaller  than  the 
reality.  Thus  while  the  model  saw  a  completely  submerged  suit,  it  also 
saw  a  thermal  resistance  probably  larger  than  reality.  The  combination  of 
these  two  "errors"  acted  to  cancel  each  other. 

The  average  of  the  observed  cool Ing  rates  was  higher  than  the  average 
of  the  model  predictions  for  the  CWU-2I/AP.  This  suit,  while  ostensibly 
a  dry  suit,  was  observed  to  leak  during  testing,  us  noted  In  Part  I  and 
In  Chapter  3  of  this  part  of  this  report.  Since  the  only  thermal  protection 
was  the  thin  dry  shell  and  arimld  underwear.  It  is  reasonable  to  assume 
that  the  majo-  portion  of  the  thermal  protection  offered  by  the  underwear 
was  lost  when  It  became  wet.  Thus  the  model  may  have  expected  more  thermal 
resistance  than  probably  existed  In  the  experimentation. 

If  one  Is  concerned  only  with  the  prediction  of  cooling  rates,  the 
model,  when  used  on  suits  that  do  not  expose  a  great  deal  of  body/devlce 
surface  to  air,  appears  acceptable.  Overall,  cooling  rates  predicted  by 
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the  model  are  faster  than  observed.  Thus,  In  general  the  cooling  rate 
obtained  by  simulation  will  lead  to  conservative  estimates  of  survival 
time  when  used  with  the  survival -t ime  prediction  model  presented  In 
Part  I  of  this  report. 

The  failure  cf  the  mode!  lies  In  Its  Inability  to  accurately  predict 
time  traces  of  body  temperatures.  Figure  l!l“7  shows  a  plot  of  rectal 
temperature  as  a  function  of  time  for  two  subjects  wearing  the  U-VIC 
Thermofloat.  The  solid  lines  are  plots  from  experimental  data.  Broken 
lines  show  the  results  obtained  when  the  physical  parameters  of  the  two 
subjects  were  put  Into  the  model.  Two  Important  observations  may  be  made. 
First,  there  Is  a  complete  absence  In  the  simulation  results  of  any  Initial 
rise  In  rectal  temperature.  Second,  the  average  slopes  of  the  observed 
and  simulated  rectal  temperature  traces  (from  maximal  temperature  to  last 
observation)  are  very  nearly  equal.  In  fact,  If  the  model  traces  are 
displaced  to  the  right  to  coincide  with  the  return  of  the  experimental  trace 
to  entering  temperature,  the  simulated  and  observed  traces  correspond  very 
well.  Since  the  slopes  are  approximately  the  same,  the  predicted  cooling 
rates  will  be  very  similar.  Thus  the  model  was  able  to  pass  the  paired 
t  test  even  though  the  predicted  and  observed  rectal  temperature  behaviors 
varied  notably. 


Tuning  the  Metabolic  Controllers 


As  noted  In  Part  I,  for  two  of  the  human  Immersion  tests  we  had  the 
use  of  a  Waters  continuous  reading  oxygen  consumption  meter  computer 
(MRM-l).  The  time  profiles  of  metabolic  rate  observed  using  this  device 
are  seen  In  Figures  1-4  and  l~5  In  Part  I.  It  may  be  seen  that  Immediately 
upon  Immersion  the  metabolic  rate  Jumps  by  about  40  keal/hr  above  a  resting 
rate  of  approximately  70  keal/hr.  Since  even  the  revised  metabolic  control 
models  did  not  show  this  Immediate  Increase,  an  experimental  controller  was 
constructed  which  Included  It  as  a  constant.  The  resulting  experimental 
controller  for  the  wet-mode  devices  Is  the  following. 


AMR  -  40  +  1.237  -  5-339  'ATth  -  4.547-ATf 


The  last  three  terms  of  this  expersslon  were  obtained  by  regression  and  had 
an  R-value  of  .69  and  a  level  of  significance  of  .04.  For  dry-mode  suits 
tested  In  1 . 7° C  water,  the  40  keal/hr  constant  was  added  to  the  regression 
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FIGURE  I  1 1-7 

RESULTS  WITH  MONTGOMERY'S  MODEL 
(WITH  MODIFIED  METABOLIC  CONTROLLER) 
AND  WITH  TWO  COLD  IMMERSIONS 


(3.)  3yniW3dW3i  1VJ.33N  Nl  39NVH3 


TIME  FROM  IMMERSION  (MIN) 


equations  determined  above  for  these  suits. 

Figure  111-8,  shows  the  profile  of  change  In  rectal  temperature  for 
two  subjects  wearing  the  U-VIC  Thermofloat.  The  observations  are  shown 
as  solid  lines  while  the  simulation  results  are  represented  by  the  broken 
lines.  The  agreement  between  observation  and  simulation  Is  much  Improved 
over  that  depicted  In  Figure  111-7  for  the  modified  controller  based  solely 
on  regression. 

Table  1 1 1  - 3  presents  the  results  of  paired  t  test  carried  out  between 
the  simulation  predicted  cooling  rates  using  these  experimental  controllers 
and  those  observed  experimentally.  The  results  are  much  the  same  as  those 
presented  for  the  modified  controllers  In  Table  I  1 1-2.  The  notable  changes 
are  the  ILC  prototype  and  Rentsch's  prototype  which  fall  now,  exhibiting 
slower  predicted  than  observed  cooling  rates  as  would  be  expected  with  a 
higher  average  metabolic  rate.  The  Increase  in  metabolic  rate  also  helped 
the  S.I.D.E.P.  Seastep  and  Stearns  Heavy-Duty  Offshore  Survival  suits  to 
pass  the  test.  Overall,  while  these  controllers  are  not  as  firmly  supported 
by  experimentation,  the  results  are  subjectively  more  satisfying. 

Figure  MI-8  presents  the  profiles  of  three  skin  temperatures  pre¬ 
dicted  by  the  model  for  an  "average"  man,  as  described  In  Part  I  of  this 
report,  versus  the  average  results  obtained  experimentally.  The  experi¬ 
mental  points  are  bracketed  by  one  standard  deviation.  The  profiles  show 
only  the  first  fifty  minutes  of  the  Immersions.  This  Is  necessitated  by  the 
removal  of  one  of  the  subjects  at  that  time. 

The  profile  of  leg  temperature  predicted  by  the  model,  as  shown  In 
Figure  I  I  I - 8 (a ) ,  lies  outside  of  the  bounds  of  the  experimental  data.  However, 
the  temperature  plotted  for  the  model  is  an  average  leg  temperature  while 
that  plotted  from  experimental  data  Is  a  thigh  temperature.  One  would  cer¬ 
tainly  expect  an  average  leg  temperature  to  be  lower  than  a  thigh  temperature. 

The  profile  of  arm  temperature  Is  shown  In  Figure  lll-8(b).  The 
initial  deviation  Is  due  most  probably  to  the  rapid  decrease  In  skin  tempera¬ 
tures  when  the  protective  device  Initially  floods  with  water.  The  model 
does  not  see  this  flooding  and  the  simulated  temperatures  drop  more  slowly. 

In  the  model  heat  must  be  lost  through  the  simulated  device  to  the  water. 

The  frequent  movement  of  the  subjects  (e.g.,  shivering)  during  the  Immersion 
helps  to  maintain  some  continual  flushing.  Thus  the  average  arm  temperature 
from  the  model  would  be  expected  to  be  somewhat  warmer  than  that  observed. 

30 


FIGURE  III -8 

RESULTS  WITH  MONTGOMERY'S  MODEL 
(WITH  EXPERIMENTAL  METABOLIC  CONTROLLER) 
AND  WITH  TWO  COLD  IMMERSIONS 


TIME  FROM  IMMERSION  (MIN) 


TABLE  III -3 

RESULTS  OF  PAIRED  t  TEST 
WITH  EXPERIMENTAL  REGRESSION-BASED  CONTROLLERS 


TEST  ARTICLES 

Calcul ated 
Statistic 

Critical 
Value  § 

0.05  Level 

— 

Degrees 

of 

Freedom 

Accept  or 
Reject 

Bay  lay  Exposure  Suit  (PVC  foam) 

4.088 

2.776 

4 

Reject 

Bayley  WeatherMate  Plus 

2.03 

2.776 

4 

Accept 

Helly-Hansen  Survival  Suit  (D-600-0) 

2.926 

2.57! 

5 

Re  Ject 

Henderson  Zip-On  Exposure  Suit  (2080-4) 

4.552 

2.776 

4 

Reject 

Henderson  Prototype  Jacket 

0.945 

2.776 

4 

Accept 

ILC  Industries  Prototype  Survival  Suit 

4.750 

2.776 

4 

Reject 

Medalist  Ski  Shorty  (7010)  with  Flight  Suit 

3.143 

2.776 

4 

Re j  ec  t 

Mustang  U-VIC  Thermofloat  ( 1 66 1 ) 

0.310 

2.776 

4 

Accept 

NADC  Goretex  Experimental  Coverall 

4.369 

2.571 

5 

Reject 

Dr.  S.  B.  Rentsch's  Prototype  Survival  Suit 

14.609 

2.776 

4 

Reject 

(without  respiratory  heat  reclamation) 
S.I.D.E.P.  Seastep  Survival  Suit 

2.597 

2.776 

4 

Accept 

Stearns  Windjammer  Jacket  (FJ-55) 

1.170 

2.776 

4 

Accept 

Stearns  Offshore  Survival  Jacket  ( FS- 500) 

1.178 

2-571 

5 

Accept 

Stearns  Heavy-Duty  Offshore  Survival 

Suit  (FS-70 

1.706 

2.776 

4 

Accept 

U.S.  Air  Force  Modified  Ant  I -Exposure 

Assembly  (CWU-2I/AP) 

1.323 

2.776 

4 

Accept 

DEVIATION  OF  TRUNK  TEMPERATURE  FROM  BASAL  (°C) 


Figure  lll-8(c)  shows  the  relationship  between  model  predicted  and 
observed  temperatures  for  the  trunk.  Experimentally  the  trunk  temperature 
was  taken  to  be  that  monitored  at  the  subscapular  site.  Again  ar  initial 
deviation  may  be  seen  which  is  easily  attributed  to  initial  flooding  of 
the  suit,  in  the  long  run  the  experimental  profile  Is  at  a  slightly  higher 
temperature  level  than  the  simulated  one.  The  simulated  data  Is  for  the 
average  trunk  temperature  which  might  reasonably  be  expected  to  be  cooler 
than  the  subscapular  site*. 
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5,0  CONCLUSIONS  AND  RECOMMENDATIONS 


Tha  potential  usefulness  of  a  simulation  model  for  evaluation  of  an 
immersion  protective  device  Is  without  question.  In  addition  to  eliminating 
the  risk  associated  with  human  experimentation*  there  ere  economic  advan¬ 
tages  In  conducting  evaluations  by  simulation  rather  chan  hy  human  sxperl- 
mentatlon.  Furthermore,  simulation  of  an  Immersion  may  be  accomplished 
In  a  short  time.  The  computer  resources  required  by  the.  mode!  are  smell, 
most  will  run  In  under  one  minute  of  CPU  time  on  an  IBM  >035  and  use  less 
than  256k  of  core.  The  cost  of  a  typical  run  made  with  Montgomery's  model 
In  this  study  was  approximate! y  $7-80.  Obviously  this  la  a  much  more 
economical  approach  than  direct  human  experimentation. 

An  additional  advantage  Met  In  evaluation  of  new  postulated  designs. 
Currently,  these  designs  must  be  fabricated  end  tested  by  human  experi¬ 
mentation.  Using  a  simulation  model  on«  need  not  physically  construct  th«s 
prototype.  Rather,  one  need  only  describe  It  to  the  simulation  model. 
Obviously,  this  approach  would  greatly  Increase  the  number  of  prototype 
devices  which  might  bn  examined  and  the  cost  would  he  much  reduced. 

The  modified  Montgomery  mod  *  I  ,  discussed  In  Charter  ’i ,  represents  a 
good  start  toward  a  useful  model.  One  may  be  reasonably  confident  of 
survival  times  calculated  from  predicted  cooling  rotes  generated  by  this 
model.  There  Is,  however,  some  room  for  Improvement.  Tho  major  areas 
Include:  Improvement  of  the  controller,  more  detail  In  the  physical  medal 
and  Improved  methodology  for  the  determination  of  protective  device  thermal 
properties. 

The  most  Immediate  need  is  for  an  Improved  controller.  Review  of  the 
literature  has  shown  that  all  control’ers  presently  In  use  were  designed 
from  data  collected  In  low  rate  of  heat  loss  experiment:,  gere^ally  nude 
cold  air  Immersions.  The  validation  of  models  with  those  controllers  has 
been  accomplished  under  the  same  conditions.  There  Is  ro  reason  to  expect 
that  a  control le"  thus  devuloued  will  function  correctly  under  the  condi¬ 
tions  one  encounters  during  cold  water  Immersions  (e.g.,  Mgh  rates  of 
heat  loss). 

As  has  bean  demonstrated,  Improved  model  performance  may  be  obtained 
by  using  a  controller  developed  from  data  collected  during  cold  water 
Immersions.  The  controllers  Investigated  In  this  study  were  all  based  on 
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simple  deviations  of  temperatures  s'  various  body  sites  from  pradescrlbad 
resting  temperatures  for  those  sites.  The  behavior  of  the  model,  using 
a  controller  based  solely  on  a  linear  regression  of  deviations  In  surface 
cemperatures  leaves  much  to  be  triad. 

It  was  seen  In  Figure  I! 1*7  that  a  basic  problem  appears  to  be  a 
too  slow  start  of  the  rwitabol  1c  lurnaca.  That  Is,  the  simulated  metabolic 
rate  did  not  Increase  sufficiently,  compared  to  the  rate  of  neat  loss,  until 
a  large  recti.!  temperature  drop  occurs.  This  leads  one  to  believe  that 
controllers  based  on  rato  of  heat  loss,  rata  of  change  of  various  body  sit* 
temperatures  or  whole  or  partial  body  heat  flux  may  be  required  In  order 
to  obtain  Improved  accuracy  In  temperature  prediction  from  the  model. 

One  of  the  reasons  Gordon's  modal  was  chosen  for  Inclusion  In  this  study 
was  because  Its  controller,  unlike  others,  used  whole  body  heat  flux  In 
determining  deviation  of  metabolic  rate  from  Its  basal  level.  Unfortunately, 
most  of  the  experimental  work  has  been  carried  out  under  low  rate  of  heat 
loss  conditions.  That  work  which  has  been  performed  under  conditions  of 
high  rates  of  heat  loss  has  not  bean  done  with  controller  formulation  In 
mind.  One  finds  that  metabolic  rates  were  not  always  taken  or  If  taken  were 
taken  at  long  Intervals.  Referring  to  Figure  1-4  we  observe  a  great  deal 
of  fluxatlon  In  metabolic  rate  as  a  function  of  time. 

Improved  controllers  can  certainly  be  developed  with  existing  data. 

More  human  Immersion  work  may  need  to  be  done  continuously  recording  meta¬ 
bolic  and  temperature  data  In  order  to  formulate  an  accurate  controller. 

An  attempt  could  then  be  made  to  correlate  data  collected  In  this  fashion 
with  data  from  low  rate  of  heat  loss  work  In  order  to  develop  a  controller 
good  for  all  modeling  dealing  with  conditions  of  haat  loss. 

As  was  discussed  In  Chapter  4,  most  of  the  suits  which  failed  the  paired 
t  test  were  believed  to  fall,  at  least  In  part,  because  of  the.  model's 
Inability  to  simulate  the  flotation  attitude  observed  In  the  experimentation. 
The  Montgomery  model  assumes  Immersion  to  the  neck.  Experimentally  (Figure 
l-l)  It  was  often  seen  that  much  of  the  trunk,  legs  and  arms  was  exposed 
to  air.  Obviously,  the  modal  should  and  did  predict  higher  than  observed 
cooling  rates.  Therefore,  to  be  applicable  to  the  evaluation  of  suits  of 
this  type  (basically  abandon-ship  suits)  the  model  should  be  modified. 

An  additional  problem  encountered  with  Montgomery's  model  was  the 
inability  to  properly  describe  some  of  the  suits.  The  mcdel  represents 
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the  human  body  .-.«  composed  of  a  complete  hand  and  trunk,  complete  arms, 

Hands,  legs  and  feat.  Thar®  Is  no  way  to  distinguish  between  thorax  and 
abdomen,  upper  and  *owwr  arm  or  upper  and  lower  legs.  As  a  consequence , 
protective  devices  which  only  partially  covered  the  extremities  were 
modeled  as  nos.  covering  them  at  all  (devices  such  ss  the  Medalist  Gkl 
Shorty  and  d-VIC  fhermo float) , 

It  Is  certainly  possible  to  nodal  this  structure  In  nor®  detail.  As, 
for  6xumpla,  Cordon'1-  model  does.  This  affords  one  the  capability  to 
moru  accurately  describe  the  physical  structure  of  the  device  to  be 
evaluated.  One  would  expect  this  to  lead  to  much  more  -«l table  predictions 
from  the  model. 

No  model  thus  far  reviewed  has  the  detail  In  the  physical  system  to 
allow  adequate  description  of  flotation  attitude.  The  Inability  Is  Inherent 
In  the  modeling  philosophy  adapted  hy  all  authors.  Specifically,  each 
body  element  has  been  modalad  <it  rtlther  a  cylinder  or  a  sphere.  In  order 
to  derive  the  differential  equations  describing  heat  conduction  within 
cylindrical  body  elements,  a  simplifying  asv.umpt  Ion  has  been  made:  all 
her:  transfer  Is  radial  and  uniform.  To  model  an  element  partially  exposed 
to  air  would  necessitate  consideration  of  iongl tud I n«I  heat  transfer.  The 
derivation  of  model  equations  including  this  consideration  would  ,iot  be 
easy. 

Fhss  changes  necessary  to  give  proper  consideration  to  flotation 
attitude  in  the  model  may  not  be  necessary,  have  seen  that  p^rts  of 
the  legs,  arms  and  trunk  are  exposed  to  air  for  some  s"lts  (Figure  1-1} . 

In  rough  sea  or  othe*‘  open  water  conditions  the  continual  washing  of  water 
over  surfaces  exposed  to  the  air  should  resui'c  In  a  rat:'  of  heat  less  very 
similar  to  that  which  would  occur  l1?  the  surf&cw  wai  continually  covered  by 
wate-.  Under  this  assumption,  tha  model's  prediction  of  cooling  rate  r.<ay 
be  very  much  In  I'ne  with  reality  for  thesa  types  of  protection  devices. 

No  method  has  been  Included  In  the  modlfrad  Montgomery  model  for  simu¬ 
lating  flushing  In  wet-mod*)  suits.  Flushing  Is  believe  to  do  two  things. 
First,  It  pnriodl-al 1y  places  a  :«rge  amount  of  cold  wator  between  suit 
and  skin  thus  causing  a  rapid  hrat  lo3$.  Second,  the  water  acts  tis  additional 
Insulation,  once  warmed,  between  the  subject  and  his  environment.  Little 
Is  known  of  the  dynamics  of  flushing  and  no  attempt  was  evade  tc  Include  Its 
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effacts  In  t h«ns«  simulations.  In  extremely  rough  water,  flushing  may 
play  a  significant  role  In  Increasing  the  rate  of  heat  loas. 

finally,  as  was  pointed  cut  In  Chapter  3,  a  more  accurate,  repeatable 
method  for  determining  device  conductivity  should  be  developed.  The 
method  suggested  In  Chapter  3  has  been  seen  to  work  well.  It  provides  a 
direct  measure  of  conductivity  obtained  by  testing  the  devices  In  the 
wet.  Testing  the  materials  under  dry  conditions  Is  not  recommended  as 
many  devices  contain  In  their  composition  porous  materials  whose  thermal 
properties  are  different  In  water  than  In  air. 

V.tille  may  refinements  may  and  should  be  made  to  Montgomery's  model, 

Its  ability  to  pass  a  significant  numbur  of  the  paired  t  tests  Indicates 
Its  potential  usefulness  for  predicting  cooling  rates  and  thus  survival 
times.  The  model  could,  with  the  present  modified  controller,  perhaps 
be  used  In  screening  devices  prior  to  In  vivo  teste.  In  the  long  term, 
a  modified  model  with  Improved  physical  definition  of  the  uody,  provisions 
for  simulating  flushing  and  a  more  accurate  controller  would  be  very  useful 
In  the  suit  design  and  evaluation  process.  The  physical  system  modification 
can  be  accomplished  easily.  The  modification  of  the  controller  may  require 
additional  human  exper  Imentat  lor.. 
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appendix  a 


THE  WINTON  MODEL 
(MODIFIED) 


UUW 


This  modal  Is  coded  In  CSHP,  a  FORTRAN -based  .general  I  zed  ,cont I nuous  - 

systems  simulation  language.  The  listing  Includes  data  required  to 

simulate  a  "normal"  man  defined  as  180.4  cm  In  height  and  weighing  73. 4  kg 

2 

with  a  surface  area  of  1.831  m  .  The  model  Is  set  up  to  simulate  Immersion 
In  a  full  body  suit  with  thermal  conductivity  of  0.499  kcal/hr  •  m  •  *C  and 
a  thickness  of  O.OO631  m  (1/4  In.).  The  following  table  Identified  the  major 
variables  of  the  model. 


MAJOR  VARIABLES  OF  THE  WINTON  MODEL 


Symbol 

TC 

TW 

TS 

SKDEL 

HE 

HW 

GA 

GW 

TE 


Definition  Units 

Temperature  of  the  core  °C 

Temperature  of  the  water  °C 

Temperature  of  skin  BC 

Deviation  of  skin  temperature  from  set  point  °C 

Evaporative  heat  loss  cal/sec 

Heat  produced  by  the  protective  device 
(If  applicable)  cal/sec 

Conductance  of  environment  cal/sec-*C 

Conductance  of  protective  device  cal/sec-°C 

Environmental  temperature  °C 


A  program  listing  follows  with  Input  data  and  output.  The  model 
Is  set-up  to  simulate  a  protective  device  O.OO63I  m  thick  with  a  conductivity 
of  0.0499. 


^wu,,p.m.,.Uju>m«u  jyn->,m  u i. 


LABEL  THERMOREGULATION  MODEL 

*  this  run  tests  water  immersion  response  of  kodfi 

*  THERMAL  ANALOG  OF  TODY 

TC  =  INTGRL  (  ICTC,-(6V/CC)*(TC“TS)  A  (HC/CC)) 

TS= INTGRL  CICTS,-(GV/CS)*(TS~TC)-<GW/CS)*<TS-TW)-<.HW/CS>) 

TW=  I  NT  r.RL  (  I  CTU,-  (GW/CW  )*  (TW  -TS)  -  <G  A  /CW  )*(T  W~TE)~  (HE  /CW)  ♦  (H  W/C  W)  ) 

*  FEED  PACK  STRUCTURE 
ELI  M=COMF-AR  (33.0,  TS  ) 

SK9F.  L  =  TS-33 .0 

W  =  L E  OL AG ( 14 00.,2  0C. , SKDCL ) 

SKI NFB  =  SGAI  N*W*SL  IM 
CORE  F  B  =  C  G  A I  N  *  TC 
TR  A  =TREF-S<  INFB 
ERROR  =  TRA-C  0REF8 

*  CONTROLLER  FUNCTIONS 
GV  *  A  F GEN (CURVE?, ERROR  ) 

HS  H  =  A  F  6 1'  N  (C  U  R  VE  3  ,  f.K  R  0 R  ) 

HC  =H  SH 

Dill  = DEBUG (3 ,4000.  ) 

DB2  =  DEBUG (5 ,6400  .  ) 

DB5  =  DEBUG  (3  ,  1  7t;00  •  > 

*  THIS  CIRCUIT  PROVIDES  A  STEP  CHANGE  IN  TE  AT  1 S  00  SEC 
Y  =  S  T  F  P  ( 1  -  JO  .  ) 

TE  -  FCN$i/(Y,  3  7.,3?.,1  .h) 

*  0  A  T A  VALUE  S 

*  JNITS  TIME  IN  SEC,  TEMP  IN  DEG  C,  H  IN  CAL/SEC,  C  IN  C AL / D EG  C 

*  JNITS  G , < , R  IN  CAL/SEC-DCG  C,  S  IN  DEG  C/MM  H6,  V  IN  M/SEC 
INC  OH  TCTC=37.0,  I  C  T  $  =  3  3  .  ),I  C  TW*3  5  .0 

PAR  f*  CC  =  58  0C0.0,  CS  *3  0  '.»G .  0,  0  '<*535 .0 
AR  A  h  SGA  IN  =  0 .02  5  ,CG  AIN  =  1  . 1 ,  TRE  F  =  37 . 1 
PAR  A, vi  HE  =0.  C,HW=(,  .0,GA  =  1  3  -U,  GW  =  1  5  .2  3 
AF  S  F  N  C  UR  VE  2  =  ( -1  0  GO .  ,  3  5  .  ,  -  .  5 , 35 .  ,  0.  ,  5.  ,  1  00 0  .  ,  5  .  ) 

AFGEN  CURVE  3=  <-10  OP.  ,1b.,  . 6,18.,1  .4  ,  100  .  ,1  000,  ,  1  0  0. ) 

METHOD  r:  K  S 

TIMER  DF  L  T-  0 . 1  ,F  I  NTI  M  =  9i.'00  .0  ,PRDEL=1  rQ.O,OUTDFL  =  1  00.0 
PR  T  PL  T  TE(32.,42. ,TC ,  T S  , T  W) 

FRII  Lf  Tc  (32  .,42.  ,HC,H£  ) 

END 

STOP 

0  U  T  3  U  T  VARIABLE  SEQUENCE 

SKDEL  7  Z  1 00  6  W  SLIM  SKINFP  TRA  CORE  F  9  ERROR  G  V  HSM 

MC  ZZ1D01  TC  ZZ1  113  TS  Y  TE  ZZ1J05  Tw  7  7  1)0 

DR1  Of?  DP  3 


A-2 


iTLOCK  DATA 

COMMGN/Z Z FO AT/F  (  78 > 

1  /ZZHlST/KEEP,NAlARM,lZODOC,IZaOC1,rt<  IS) 

2  HI  I STO  /  1  (  AO) 

COMMON/Z7rOIN/NP(16) 

INTEGER  NP/  A/  AO,  1  5,  27,  37,  3V,  5A,  5A 

1,  53,  S3,  IS,  3,  78,  0,  1,  A/ 

COMMON/Z  Z  SYMB/S  1 1  (  AO) 

RFAL  *8  SI  1/  *  TIME  ZZ  ’ 

1,  ’  T  I  MTDE  L  T  * ,  *  ZZOELT  *  ,’&  ELMINZZ  ’,  ‘DFINF  INT  •  ,  •  I  »*ZZF  INT  ‘  PPDE  L  ZZ‘ 
1  , ’PR  DEOU  T  D*  ,  •FIZZOUTP"  ,  *D  ELMAXZZ  ’  OfLXTC  *,*  TS  ',*TW  ZZ  ' 

1,,10C7ZZ10,,,01ZZ10O3',,Z  Z10052Z  1 ,  '  1  C 36 1 CTC  1  ,  '  JCTS  8,’1CTW  ZZ* 
l,*10r.VCC  CGA1N  ',!CS  CW8,8  GA  •  ,  9  GW  *,'HE  HW 

1,'  S3  A  1  '  ,  *N  TREF  *  ,  8  C  UP  V  E  2CU  * ,  8  R  VF  3CORE  8  ,  *  F  BDB 1  8,8r>B?  D3  8 

1  ,  8  3  E  R  R  0  '  ,  '  R  GV  8,8RC  HS*,'H  SKDE',,L  SK  I  N  FB  ' ,  '  SL  I M  TP8 

1,  '  T R  A  '  ,  '  W  8 , 8r  * 


SUBROUTINE  UPDATE 
COMMON  TIME 


1,ZZOOCO,DtLT  , Z  ZD  E  L  T 

,DFLMIN,ZZDELN 

,F I NTIM 

,Z?FI NT 

,  P  R  D  E  l. 

,  ZZPR  DE 

1,OUTDEL,7ZOUTD,DELMAX 

,Z7DFLX,TC 

,TS 

,  T  4 

,z  zi  cr 

7,  771  3  01 

1,ZZ1CC3,7Z10GS,ZZ190* 

, ICTC  , I C  T  S 

,  IC  TW 

,  Z  Z  1  0  09  ,  C  C 

, CG  A  I  N 

1  ,  C  S  ,  C  V!  ,  G  A 

,  GU  ,  H  F. 

,  HW 

, S3 AI  N 

,  T  R  L  F 

, C  UR  V  E  2 

1, CURVE 3, C ORE  FB,Dni 

,  DP2  , DD3 

, ERROR 

,  G  V 

,HC 

,  H  S  H 

1  , SKD  EL  , SKI NFB,SLIM 

,  T  F.  ,  TR  A 

#  w 

,  Y 

COMMON/7  ? HI S  T / K t  E  P,  N ALARM #17  OOOO ,  I  Z  OH0 1 
REAL  ICTC 
1,1  CTS  ,1  CTW 
REAL*R  ZZTIME 
EQUI  VALENCE  (ZZT1ME,TIME  ) 

GO  TO  <39  0  95 ,399  9  6,  399  9?,  3  999  8),  I  7  0  DOC 
C  SYSTEM  SEGMENT  Of  MOTEL 
39995  CONTINUE 

50  TO  39  9  99 

C  INITIAL  SEGMENT  Of  MODEL 
39096  CONTINUE 

GO  TO  39999 

C  DYNAMIC  SEGMENT  OF  MODEL 
39997  CONTINUE 

SK  DE  L  =TS  -  33 . 0 
ZZlrC6=<SKDEL-7ZlO0  7)/?0Q  . 

4  =  1  A  CO.*  2  Z1  0C6*ZZlQL'-7 
SLIM  =  COM P AR  <33.C,TS  ) 

SK IN  FB  =  S5AI N*W*  SLIM 


C 

c 


TRA  =  TREF-SKINFB 
CORE  FB  =  CO  AIN*TC 
ERROR =TRA  -COREFf? 

3V--AfGEN(  1  ,CIJRVc2  ,CR  POR) 

HSH  =  AFGEN  (  6  ,C  UP  VE  3  ,  -  RRO  R  ) 


HC=H  SH 

u  in  :>i 

Tf. 

ZZ  1  .03 

TS 

Y  =  S  T  E  E  ( 


=  - ( G V / C C )*<TC- TS)  +  (hC  /  CC  1 
=  INI GRL  (ICTC  »  Z  Z  1  Ct  1  ) 

=  -<G  V/CS  >M  TJ>  TC)  -(OL'/C  S)  *  (  TS-TW)  -<M  W/CS) 
=INTGRL  (I  CTS  ,ZZ1°f'3  ) 

14,1  POO.  > 
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C  TV  =  I N  T  G  R  l  ( 1  CT W  ,  Z  7  1  '  US  ) 

C  Z  Z  1  ‘  C 7  =  INT'-RL  <Z  Zlt  09  ,72V  -:6  > 

DB1  =  DEBUG  (  11  ,3, 4  "'DO.) 

DB  2  =  D  £BU  5  (  12,3,64(  0.) 

DO  3  =  DEF'U  G  (  1  3,3,  1  7000.) 

GO  TO  39  9  99 
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APPENDIX  B 


THE  GORDON  MODEL 
(UNMODIFIED) 


This  model  Is  coded  In  standard  FORTRAN  and  has  been  run  on  an 
IBM  370/168  and  an  IBM  3033-  Run  time  on  the  IBM  3033  Is  about  two 
CPU  minutes  for  a  simulated  experiment  of  three  hours  length.  As  stated 
In  the  body  of  this  report,  the  mode!  never  performed  properly.  The 
only  modification  appearing  In  the  listing  was  made  for  ease  of  Implemen 
tatlon.  The  following  table  lists  the  major  variables  of  the  model. 


MAJOR  VARIABLES  OF  THE  GORDON  MODEL 


Symbol 

Definition 

Dimension 

C  (N) 

Heat  capacitance  of  compartment  N 

kcal  .  °C 

T  (N) 

Temperature  of  N 

°C 

F  (N) 

Rate  of  change  of  temperature  In  N 

°C  .  h-1 

HF  (N) 

Rjte  of  beet  flow  into  or  from  N 

kcal  .  h 

TC  (N) 

Thermal  conductance  between  N  and  N  +  1 

kcal  .  h 

TD  (N) 

Conductive  heat  transfer  between  N  and  N+l 

kcal  .  h 

QB  (N) 

Basal  metabolic  heat  production  In  N 

kcal  .  h 

Q  (N) 

Total  metabolic  heat  production  in  N 

kcal  .  h 

EB  (N) 

Basal  evaporative  heat  loss  frcm  N 

kcal  .  h 

E  (N) 

Total  evaporative  heat  loss  from  N 

kcal  .  h 

BFB  (N) 

Basal  effective  blood  flow  to  N 

1  .  h“' 

BF  (N) 

Total  effective  blood  flow  to  N 

1.  h’1 

BC  (N) 

Convective  heat  transfer  between  central 
blood  and  N 

kcal  .  h 

HC  (1) 

Convective  and  conductive  heat  transfer 
coefficient  for  Segment  1 

kcal  .  m 

s  (1) 

Surface  area  of  Segment  1 

2 

m 

HR  (1) 

Radiant  heat  transfer  coeff.  for  Segment  1 

kcal  .  m 

H  (1) 

Total  environmental  heat  transfer  coeff. 
for  Segment  1 

kcal  .  h 
-1 

m  .  sec 

V 

Air  velocity 

TAIR 

Effective  environmental  temperature 

°C 

RH 

Relative  humidity  In  environment 

TIME 

Elapsed  time 

h 

PAIR 

Vapor  pressure  In  environment 

mm  Hg 

(TIME 

Elapsed  time 

min 

!NT 

Interval  between  outputs 

min 

B*  I 


Symbo  1 

Definition 

Dimers  Ion 

OT 

Integration  step 

h 

p  05 

Van or  pressure  table  from  5  “  50* C 

■nm  Hg 

•I  MAX 

(1) 

Calc,  max.  rata  of  evaporative  heat  Ions 
fiom  Segment  1 

kcal  .  h  * 

WORK 

Total  metabolic  rate  required  by  exercise 

,r  | 

kcal  .  h 

PSKIN 

(1) 

Saturated  water  vapor  pressure  at  skin  temp. 

mm  Hg 

't  SET 

(N) 

"Set  point"  pr  reference  point  for  receptors 

®C 

ERROR 

(N) 

Output  from  thermoreceptors  In  compartment  N 

Jc 

RATE 

(N) 

Dynamic  sensitivity  of  thermoreceptors  In  N 

h 

COLO 

(N) 

Output  from  cold  receptors  In  N 

°C 

WARM 

(N) 

Output  from  warm  receptors  In  N 

®c 

COLDS 

'ntograted  output  from  sk*n  r.old  receptors 

°c 

WARMS 

'ntegrated  output  from  skin  warm  receptors 

°c 

w  1 

SWEAT 

Total  efferent  sweat  command 

kcal  .  h 
kcal .  h  * 

CHILL 

Total  efferent  shivering  command 

D 1  LA  f 

Total  efferent  3kln  vasodilatation  command 

<  .  h"1 

STRiC 

Total  efferent  skin  vasoconstriction  command 

n.  n. 

SKI  NR 

(0 

Fraction  of  all  skin  receptors  In  Segmrnt  1 

N.  D. 

SKINS 

(1) 

Fraction  of  sweating  command  applicable  to 
skin  of  Segment  1 

N.  D. 

SKINV 

(1) 

Fraction  of  vasodilatation  command  applicable 
to  skin  of  Segment  1 

N.  D. 

SKINC 

(1) 

Fraction  of  vasoconstr Ictlon  command  applicable 
to  skin  of  Segment  1 

N.  D. 

MWORK 

(1) 

Fraction  of  total  work  done  by  muscles  In 

Segment  1 

N.  D. 

MCHIL 

(1) 

Fraction  of  total  shivering  occurring  In  muscles 
of  Segment  1 

N.  D. 

-  t 

CSW 

Sweacing  from  head  core 

kcal  .  h  . 

_  > 

ssw 

Sweating  from  skin 

kcal  .  h  ."C 

—  1  -  I 

COIL 

vasodilatation  from  head  core 

i  .  h  .°C 

SOIL 

Vasodilatation  from  skin 

1  .  h  1  .ec  1 

CCON 

Vasoconstriction  from  host!  core 

•c"1 

1 

SCON 

Vasoconstriction  from  skin 

®  r  ™ 

v> 

1 

CCHIL 

Sblvtrlng  from  htad  core 

kcal  .  h  .* 

kcal  .  h  *  .  *C 

SCHIL 

Shivering  from  skin 

PSW 

Sweating  from  skin  and  head  core 

kcal  .  h'1  ,*C 

B-2 


Symbol 

PDIL 

PCON 

PCHIL 

BULL 


Def  I  n 1 1 Ion 

Vasodilatation  from  skin  and  head  coco 

Vasoconstriction  from  skin  and  head  core 

Shivering  from  skin  and  head  core 

Factor  determining  temperature 

sensitivity  of  sweat  gland  response 


Dlmens ion 

•MMMM  111  W«MIILWI'»ll 

I  .  h~  . °C 


-i 


rmmmmmw  K’ra'WT'  ‘^WTwt  ",l "' "  ’fi’-’t  ■wrr^S’V?- : T7n- 


C 


c 


IMPLSC  IT  NEAI.  *8  (A-M*0-J  ) 

pi  m  fns i on  sain  < 3 ■  t ( ) ,  i  /.» ,  io<3'  n,  73),  i»i  n  <in>  ,\  «ax  <id  , 

*AT  <54,  ?  1  1,1  N  D(34>  ,GE0<3Ml?  >  ,V  <  ?14>  ,YY  (34  ,?1  ),FF’S(3  4 

*  )  -PC  (3  4  ,  1 1 )  ,  U«  <  34  *  V)  ,R  (  34,?  1  ),  V  <34  ,  ?';)  ,  VT  <  34  ,9>  ,  AC  (  T, 
*4)  ^RMO  <  3 4 , <-> >  ,C  (  ,V.„9y  ,AK  (34,  0)  ,0M<34  ,?H  )  ,  CME  T  <34, 9),  V 
*EP  (3*«,  9  )  ,  TAR  (74  )  ,  A  (  74,  ?  0,0  <  34,  ?0>  ,  D  (3  4,?U>,  t  (74,cr  ) 
*,T(34,  ?  1  >,1  P  (I  ?  O  ,  TV*N  (  34)  ,U  A'  34},  VtL  N  (34  )  ,*  i1(  1?6>, 

*  CNGR(  3  4,?T  ,  FPT  (34  >,09(  34,9)  ,BFP<34 ,9)  ,TV  (34  >  -  -MRMi  3 
*4  )  /  COL  D  (34)  ,  FLUX  (  H  )  ,t  ff  ROR  (  3  A  )  ,  SK  1  N  H'  ( 3  4  )  *  SK  I  MS  1 34  )  ,  S 
*KI  NV<3  4  )  ,$K  I  TC  l  74) , WORK  M< 34)  ,C  M  M.M  <  34)  ,E  8  (  34  >,F  SE  T(  3 
*4  )  /PS<  INC 74 ) ,TO ( 34/ ?1 ? * OAT* ■  (  74 ) ,T0l  M 34)  , I  f  (  34) ,EMA  X 

*  <  34  )  ,4  N  A  MM  .3  4, 3  )  ,  I  SK  (  34  \  ,  I M  U  (.7  4  )  .  I  W  CR  M(  34  >  ,  r  COL  0  (  34  ) 
COMMON  /  r.NV/BO/l  DR/OR/LR  /R/L  V  /V/LVT  CVT/14C  /  AC  /LRHO/P  H 

*0/LC/C/LAK/AK/Lr-M/HM/(.CMtT/CMET/LVfr/VCF./LTAR/TA»/LB 
*L  0  /  R  ID  U  ,  CH  ,  T  AF  F  /  C  OA  /  A  /C  00/  8  /  CO  0  /D  t  C  OF  /  fc  /  l.  T  /  T  /l  I  D  /  NIJ  Y 

*  ,  I  S  pH f  R  ,NUM  S  ,1  S  S,NOP, t  D/LTVF  Hi  YVCN/  LUA/t'A  /  LV  ELfl/ VFl  P 

*n  I  4/1  A /LtNGR/t  NGR/LFBT  /til  T  /  LOB /G3  /  LB  FB/ B  F  9 /  L T I M F  /CO 
*,  FT?  ME  ,  FWRITF/L  T  V/T  V  /L  TOT/TOTV,  TOT  A/LWAKM  'WA  RM/LCOL  P 
*/{<  LO/L  F/F7L  EkROR/ERROR /I.SK  I  ’1R/SK!  IH7LSICNS/  SKINS/IS 
'K  IMV/SK  INV/  L  SKI  I  *C  /  '.H  0  P  KM / WO R  KM/ L  CHILM/CHHM/IE 

*B/FB/LHR/HR/LHr/HC/LPSKIN/PSKlN/LT(7T0/LPATF/R/T»7LT 
*OLP/TOLP/LF  E  /F  t-  / 1.  E  AiA  X  / E  MAX  / 1  IMU  /  I  MU  /  l  F  LUX  /FLUX/STOlP 
*X  /TA1R  ,  'll  ,KH  #WQKK  ,  f :CMF_  T  cCTOL  0/  L  NAM  t  /  A  Nl  A  M  C  /’II  5K/ISK/L 

*  5T0L/C  S  W,S$ W ,PS V, CD  tl  ,  SO  II  ,  PDI  L  ,  C  C  0  H,  S  C  ON  ,  PC  ON,  C  CHI  L 
*,SCM  IL  ,  PC  HI  L  ,PULI  ,  AOL/P  CO/LO  TEN  /«T  t  0/L  FSf  T  /F  Sf  T/LFW  A 
*RM/F  W4  RM  /LF  COLD  /  F  C OL  f'/L  QTHL  /  C.TPL/LS  TOP  /STOPI  T/LMUS  /  M 

*us 

LOGICAL  CSYFL 

C 


c 

c 

c 

c 

e 

c 

c 


<  ><>' 


<><><><><> 


<  >  <> 

<>  INPUT  AND  OUTPUT  OF  MOOr.L  ST/fT,  MATERIAL,  AND  <> 

<>  NODAL  CONFIGURATION.  <> 

<  >  <  > 


■<><>  <><><><  ><> 


DATA  T I KF,I Z ,VH , IU,VT, I DA/n . PO COM , 2 ,0 . 0, 0 ,0 . n, 3 / 
PFAD(Wft1>TIM$T,I  SPUE  R,  M  US 
61  FORMAT ( F TO. 7, cl f) 

VPI1  £(  3  #  31 )  7  IMST^NUN#  I  S  PHfcK  ,  KIIS 
31  FORMAT  (  ■1,,T  20,  'CXPERIMFNT  DURATION  =’  ,F1'J.Z  ,'  HOURS1 
*//T20,  *  TOTAL  NUPPLR  OF  OODY  FLEMENTS  EQUALS  :*,I5//T?0 


R-l, 


* »  '  T  0  T  4 1  NUMBEh  Of  SPH  F  R 1  CAL  ELEMENTS  EOUALS  ;'xI5//T20 
'r.MREMITJ  f  ;>  STANi  WITH  ROOT  ELFMfNT  : 1  ^  I  5  /  /  > 

PF  AD  <1  x  32)  (  I  A(M  ),N*1,JO 
3  2  FORMATU012) 

READ  <1  x  33)  U  D(M  )  xM*1  rU  ) 

33  FORMAT  (H9) 

DO  2  5  M=TxNU/ 

IFUA(H).NE.C)  HOP*'? 

25  CONTINUE 

RE  AD  <1  x62)  (li  A‘Vl )  xM  =  1  x  H  > 

RE:  A  D  ( 1  /  133)  (  IMU(M)xvi*1  x  3 #*> 

133  FORMAT ( 1615) 

SE  A0<1  x35)  C  (  »NAMf.(MxlO  *  K=1,  3  jxM=d/34> 

35  FORMAT (  1 8 A4  ) 

WR]  TE(  3x34MMxI  A(M)  x  I  »  <M),AN  AKE  (Mx  1  )  xANAME  <M  r2>  xANAVF  < 
*M,3>  xUA  <M)x  IMU(M>xM-1,.U; 

3*4  FixRMhT  (1204.  •  FLEMF.NY  S  YMC'ET  P  IC  A  I.  NO  OP  *  .  1  0  X,  *  AN  AT  P  ’  , 

*  1  f* )  C  i\L  CCHf.  (KCAL/HR/t-E  C  >  M  IT  CL  E  •  /  T2?  x  * NO  .  *x 

*'  ELEMENT  DISTRIBUTION  NAMF  COL K  F  ’  x 

3X,  ’D  I  STRI  PUT!  OH  •  //  3',  (  T  ?2  x  I  2x733/  I  2xT4f  <  1 1  0  ,1  5?<  3A4<  T 

*  76#  F  5 .  2xT96x  14/ )) 


C 

C 

c 

c 

c 

c 

c 

c 


<> 

<  > 

<  > 

<><><><><  ><> 


INPUT  AND  OUTPUT  OF  GEOMETRIC  S  I  Z  F  C  <• 
PARAMc’iJ.RSx  AND  30UNDRY  CONDITIONS. 


><  ><><><><• 


<><><><><> 
<> 

D A  5 A L  BO-5  V  <> 
<> 
<  > 

<><<><><><>^d<  ><><><><>  <><><  ><> 


CALL  3E0MRG (GEOxNUMx-Z) 

WRITE! 3^66) 

G6  FORMAT  (  *1  'xT?Ox  "F<C  Cl  1*  =  CONVECTIVE  H.T„  EOE  FF  .  (CAlr, 

* ' /SO .CM /HR/ DFC) ' /T20x J FPS  =  EM  I  SSI  V 1  T  Y  *SH  A  0c  -F  A  C  TOP  'x 
*'  ffl.D.  /  ,/T2f'x,RC  (3)  =  T  t.KPE  R  ATURE  OF  CONVERT  1VE  EMVl'x 
CENTIGRADE  >'/ T  2C>x 'PC  (  4)  =  TTMPfcRC  TURF  'x 
ENVIRONMENT  (DFG  C  EN  T  IG  R  Ah  E  )  *  /  T?  ''x  *  2C  I’  x 
TEMPERATURE  CHANGE  PERMISSIBLE  AT  EKD  'x 
(DEG  CFN  j  IGRADr)  ' //T7x*  PLEMENT 
B  C  (  3  )  B  C  ( 4  >  B  C  (  7  )  ' ) 


RC  (  1  '  x 


0  f  <  M  x4  )  xB  C  (1x7) 


*  '  RON  ME  N  1  (DFG 
*J0F  RADIATIVE 

*  '  7.)  =  INITIAL 

*  *  0  F  ITERATION 
*'1)  FPS 

DO  2?  M-1/N0F 
I F( ASY  FL  (M) ‘  GO  TO  T 
READ  (1  x67)3C  (Mx  11)  »EPS<  MTxBC  (f?,  T) 

67  FORMAT  (4F1'1. 5x1  >  xF  U.9> 

WRI  TF<  3x64)  A  NAM  r  CM  x  -  )  x  A  N  AM  F  <  M  j  ?  )  x  A  N  AM  E  (M  x  3  )  ,  PC  (  M  x  1 1  )  x  * 
*PS  (M)xBCU-x  3)xPC(Mx4)  x«TC  (Mx  7  ) 

64  FORMAT  <  1Xx3A*-x1X  »5F  10.5) 

9C  (Mx3T  -BC(Mx3)  *27  3.  16 
BC  CM, 4  )  =BC«M.-4)  >273. 16 
2  0  CONTINUE 
URITFC  3x66) 

6  5  FORMAT  <  /  //riCx'RHO  =  DENSITY  (GM /CUBIC 
*'  SPECIFIC  H5AV  <  C  AL /CM  /  r P.G  )  V  T 1  Jo  '*•! 


CM)  ’/  T1  :• »  ‘f 
CCND  or  TI  VE 


*’.T.  COErHTJFf.’T  (CAL  /l‘M/Hk/!>F  G>*/T  10,  >«.MET 
* 'MET  seOL’SM  (CAL/U*- .•  VT  lOx'VIB  =  BASAL  PLOCD 


H  '  x 

*  BASAL 
rLOW  (M‘, 


e-s 


&B5i  . 


• '  L-’KWfcSWS*  Sit 


' '  •■.-Aatyi - 


,^x-j  n  ymr,';™rffl®1^?’- 


l,  :t^,-r,U w ‘ a Ht  ^ ' 


6? 

1  1 

6  3 


1  0 


7  A 


*  *  l, /MR)  •  unt,,'  F  CEMENT  MAIL  NO,  K  HO  •  ,  1  2  X  #  *  C  ’  ,  1  <'Y  , 

*  ’  A  K  *  #1  0  >  *■  *  C  M  f.  T  ’  #  1  G  X  *  *  V  C  B  *  ) 

R HOC  f»\ /  -o.r, 

DO  17  M  « 1  ,  N  0  P 

<)A(M)=UMM)  *1GCC. 

IMASYFLOl))  00  TO  1?. 

CALI.  D  E  F3RC,  (M/LMAT^MODL  ,11  ) 

I  NO (M)  -NODE 
IFK  (M>  =  L VAT 

SIREPsEPS  (M  >  V,  . F756E-9 
BC(M,2 )  -S IGF  F*A  C (H) 

DC  11  J=1,LMAT 

PEAD(1#e?)RH0CM#J)#C(M/.J),AK:<«^J)^CriET<'1/J)/VE3(M/J> 

Rh^CMV  =  RHOC  PV+RMO(M*  J  )  *  C(M,  J  >*VT  <M,J) 

FORMAT  (5F12.5) 

P  F  EH  M.  J  )  -  V  E  B  J  ) 

OB (M  *>  )  =  CMF  T  (Mi»  J  > 

VTsWTf  RHCXM,  J)*  VT<M,.J)  /  1000. 001 

W R  J  T  F  (  3  ^  63)  A  NAM  E  (M 1  )  ,  A  N AME  (  M,  2  )  ,  ANAMF  <m,  3)  -  J,«  HO  (M  .  J  ) 
*,  C(M/J  )  *  AK!  (M  ,J  >  ,  CM  FT  CM#J  ),VEf  CM  ,i  ) 

FORMAT  (  1X^3A4,4  X,I  2,7:X,  5<1X  /•  F 1  ?  .5*  IX)) 

IF  (M.3T  .  ISPHER)  GO  TO  10 
ABR-RC  M  *  NOD  F  )+DR(M#LMAT) 

ARRR=R  (  M,NOP  F>  +  D «  (  M ,  LM A  T  )  t  2  . 

HC  (M,9)  =-DR(M,LMAT)/«K(K,LMAT)  /AC(M)*APH*  APR  /APRR/APRR 
GO  TO  12 

AC  C  = AC  (M  )/R  <M,NODt:> 

DC  (?*?)  ~-D  R  (  M/-LNA  T  )/  <R(  MrNOD  C)  +  DR(MHNAT)/2.  )/ABC/AK(M 
*,LMA  T) 

CONTINUE 
WRITE  (3W4)  WT 

FORMAT  (///T2C/iTME  TOTAL  BO  D  Y  WEIGHT*/ T  ?  0  *  *  E  DUALS  -  =>  * 

*  ,  f  1  0 . 5  ,  ’  <  <  G  )  *  /  /  > 


C 

c 

c 

c 

c 

c 

c 


<  >  <> 
<>  CALCULATION  OF  CO C F F I C  I  TN T S  FOR  ITERATIVE  ROUTINE  <> 
<  >  •  <  > 

CALL  COtFSG 

RE  AD  (1  /-1D2)  C  tfH  IM)  *1*1=1  f  NUM  S  I 
ID?  FORMAT ( 16f5 . D 

READ  (1  y  1  60)  TA  J  P,VL,RH„  WORK  #  D  l  N  T 
15  0  FORMAT  <3  F10. 5,  F  F1P  .fe,?J  7) 

IK-  0 

DO  3  02  N  =  1#N0P 

1F( AST FL<M> >  CO  TO  30? 

I  M  =  I  fit  1 
BUM, 8 )=RH 

fH  (M,5  > -DINT 

BC  <M,1  )  =BC  <W  ,1  1  >  *AC  <M> 

33?  CON  T I S  UE 


B-6 


Tl 


C<>  <> 

C  <>  STEADY  STATE  TEMPERATURE  D  1STKI8UTI0N  AND  BODY  <> 

C  <>  PARAMETERS  ARC.  CALCULATED  AND  PRINTED.  <> 

CO  <> 

C 

CALL  S  TE  PRG  C  YY*  Y  ) 

WRI TE  C  3  *73) 


73  FORMAT ( //IX *  THE  ARTERIAL  PLOOD  TEMPERATURES  ARE  A", 
*  '  S  FOL  LOWS:  ' //  ) 

r 


C 

I  M  “  0 

DO  15  M  =  1/M0P 
IF(ASYFL(M) )  GO  TO  15 
IM=IM+ 1 
L  5  K  =  I  $  K  CM) 

L=IND(M)+1 
FSETCI  M)=FLUX(IM) 

DO  5  M  =  1  *L 
5  T  0  C  M  *  V  )=T  <  M  *  M ) 

TOLD (I M) =T(M*L-  1) 

V  M  s  V  B  ♦  VT  CM*  l  SK) 

AT  M)  =  TAR  (M)  -273. 16 
15  CONTINUE 

C  T  0  L  D  =  T  C  1  *  1  ) 

WRI TEC  3*62)  ' AT  C 1*IM> ,IM  =  1,NUMS ) 

RFADC1  *  10?)  STOPIT 
W R I  TEC 3  *  24  )  STOP  IT 

24  FORMAT  C//T20*' SUCCESSIVE  ITERATIVE  TEST  EQUALS  =  =>*,F5 
*.4) 

10  3  FORMAT  C16F5.4) 

RFADC1  *  103)  C  SKINRCID*IM=1*NUMS  ) 

READ  Cl  *103)  CSKIN5C1M)*1M  =  1*NUMS) 

RE  AD  Cl  *  103)  C  SKINV  C I  D  *  I  M=1*N  UMS  ) 

RFSDC1  *103)  CSKINCCIM)*IM  =  1*NUMS) 

RE  AD  Cl  *  1  03)  C  WORKMC  ID*  I  M  - 1 ,  N  UM  S  ) 

READ  <1  *  103)  C  CHI  LMC  ID*  IM  =  1*N  CHS  ) 

WR!  TEC  3*  104) 

IT  4  FORMAT  C  '1  '//  '  ER=CAL/H^■,/,  F  L  UX  =  f  AL  /  HR  l  SO.  CM  '//WO* 

*' ENVIRONMENT  AL  CO  f  F  F  .  *  *  T  30*  ’  WF  I  GH  T  I NG  COEFFICIENTS*// 
*  6  X  *  'ELEMFNT  SKINR  SKINS  SKINV  SKINC  WOK", 

*'KM  C HILM*  *  10X * ' FLUX ' *  11X*  ' EB  1  ) 

I M  -  0 

DC  1  05  M  =  1  *  N  OP 
I F  C  ASY  FL  (M)  )  CO  TO  1  '5 
IM= I M* 1 

WRITEC  3*1('7)  ANA  ME  (M*1  )  ,  ANAM  T  <M*?)*ANAMECM*3>  *SKINR(IM) 

B-7 


*#  SKINS  <  IM>*  SKIN  V(IM)  *SK  INC  (  IPO  #WORKM(IM>*  CM  LM<1M>*  FLU 

*x  um)«  eb<  in> 

ITS  CONTINUE 

ID  7  FORMAT  <1X*3A4*SX*6F7.4*E?0.?*F7.()> 

PE  AO  <1  *  1 1 0 >  C  5W,SSW*PSW*CDIL*  SD IL*PDIL*CCON*S  COM  *  PCO  N*  C 
*CHIL *SCHIL*PCHil *B UL L *  A 8L* Q T EM* QTR L *  PC  0 
11  0  FORMAT  (<SF0.3  ) 

W  P  I T  E  (  3/11?)  CSW*SSW*PSW*CD1  L  *S  OIL  *PDI  l  *  f  CON#  SCON  *  PC  ON  * 
*CCHIL*  SCH1L  *  PC/(  1 L  *  BUL  L  *  AHL  *3TEN*QTBL*PCO 
11?  FORMAT  (  1  5  (  /  )  *1X*T20**CONTROLLEk  COEFFICIENTS  *///T?5* 
*'CSW  =  ' * F9.4*  '  KCAL/HM/DEG ' / T?5* 1 SSW=* *F9.4,  '  SO.M’/T 
*22*'  pc.W«  VF9.4#  *  LlTFR-SQ.M/KCALVT?5*aCDIL***F9.4 
LI  TFR/hR/DEG '/T25*  '  $DIL=  **F0.4*  '  LITER’-SO.M/KCA'* 
*,L,/T2  5*,P0IL=,*F9.4*,  LI  T  fc  R-SQ  -  M /KC  AL  '  /  T?5  *  *  C  CON=  *  *  F 
*9.4,'  1/DEG  '/T?5,  'SCON  =  ‘*r  «  .4#  '  L  1TE  R--SQ*  R /KC  AL  '  /  T2  5 
**,PC0N5:•*F9.4*,  LITER-SQ.NI/KCAL,/T?5*,CCHIL=,*F9.4*2X 

*  *  'KCAL/NR/PFG'/TPS*  *  S  OH  1 L  =  *  *  F9 . 4*  *  SQ .1 ' /T2  5,  ■  PCHI  ' 

**F9.4*'  KCAL/MR/SQ.OEG  ' /T2 S *  *  BULL = ' * F 9 . 4  * ’  1/0EG'/T2 

*5* 'ARL= ’ *F9. 3* •  LITER/O PG /HR' /  T?  b  * ' QT  F  N  =  '  *F 9.4/’  1*, 

*  *  /DEG*  /  T25*  '  GTBL  ='  *F  '5.4  *  1  1  /  DEG  '  /  T  2  5  *  '  P  C  0  =  *  *  F9.4,  •  tf", 

* 'L/C AL •  /  //) 

C  S  W=  C$  W  *  1  0  TO  . 

SSW=SSW*10PC) 0. 

psw=ps  w*iO(;oo. 

CDIL=CD  IL*109C. 

SOI  L=S  OIL*  10  000. 

P  01 L  =p  0 1  L  *1  0  uOO  . 

SCON=SCON*lOCOO  . 

pcom=pcon*i o;:oo . 

CCHIL  =  C  CHIL  *1000. 
schil=schil*ioo:o. 

PCHIL=PCHIL* 1  TOT. 

A9L  =  ARL  *1000  . 

GO  TO  44 
90  COnTIMUF 

c 

c  <> 
c  <> 
c  <> 
c  <> 
c  <> 
c  <> 
c 

READ (1*150)  TAI P* VL*RH* WORK *DJNT*iBC*IDA 
DO  ? 0?  K  =  1*N(iP 
IF(ASTFL(M))  GO  TO  2r‘2 
IF(IPC.NE.I)  GO  TO  22 
IF(IOA.EO.n)  GO  TO  23 

RF  AD  <1  *  1  3>E  PS(M)*PC(Pl*3  )*RC  ( M,  1 1)*PC(M*8  )  *RC  <M*  10) 

1  3  FORMAT  (F5.4*f5.?*F8.4*F6.4*FA.4) 
BC(M,2)=rPS(M)*AC(M)*4.R756E-9 
BC  < M *3  )  =BC  (M*3>  -»2/3. 16 
GO  TO  23 


INPUT  AND  OUTPUT  OF 
F)OUN  DR  Y  CONDITIONS. 


TR  ANSI  0 NT 


<> 

<> 

<> 

<> 
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22  RC(M,8)=RH 

DC (M,3  >  =TAIR  +2/3. 16 
BC(M,1  f)  )  =VL 

23  BC ( M,7  ) =1.000000000001 

BC ( M,5  )  =  DINT 

FC(IM,4 ) =  PC(M  ,3) 

BC  (M,1  >  =  BC(M,1 1  >*AC  <M)  *  (RC(M  ,1  0)/.  10>*  *0.5 
212  CONIJMUE 


C 

C 

C 

C 

C 

C 

C 

C 


<><><><>< 

<> 


<  > 
<  > 
<  > 


TRANSIENT  DETERMINATION  OF  TEMPERATURE 
DISTRIBUTION  AND  PHYSIOLOGICAL  PARAMETERS 


<  > 
<> 
<> 
<  > 

<><><><><  > 


CALL  S  T  EPRG (  YY, Y  ) 
T  I  ME  =T  IME  +  EJC  (1/5) 


44  I  Q  =  I  Q*  1 


I  F  (  TIM  F.  .GE  .  T  IMS  T)  GO  TO  82 
F-  T I  ME  =  T  IMF  *  6  G.  G  0001 
TPOOL=TARP-273. 16 


C 

C 

C 

C 

C 

c 

c 

c 


<  >  <><><> 
<  > 

<  > 

<> 

<  > 


OUTPUT  OF  TRANSIENT  P  H  Y  S  I  C  LO  S  I  C  A  L 
PARAMETERS  FOR  EACH  T I  NF  STEP 


<><><><  ><> 
<  > 
<> 
<> 
<> 

<><><><  ><> 


WRITE!  3#  40)  ETI ME  ,  T°OOL,ECMFT 
4  0  FORMAT  (//IX/  *T  l  V  t  EQUALS  •  ,  F5 . 1,  *  MINUTES* /IX, 'ARTE 
*  '  R I  AL  POOL  T  TM  P  F.  P  A  T  U  R  E  EQUALS  *,F6.3,'  OTG  '/ 1  X  ,  '  P  L 
♦'SPIRATORY  LOSSES  EQUALS  *,-3PF6.3,*  KCAL/HR*/) 

WRI TEC  3,42) 


FORMAT ( KY, ' E LEM 

TC  ORE 

TSX  IN 

INBAL 

RBEC 

BLHTR 

*  *  MET 

CONV 

RAD 

EVA  P 

SBF 

TBF 

QCC 

*  *  T»R 

FLUX 

FOT 

FOOL  *  ) 

SCUT  (I  Q  ,  1  )  =  F  TIMF 
SOUTd  Q  ,13)  =  C.O 


EWE  T=3 . 0 


SOUTd  Q,12)  =0,0 
SOUT  (I  Q,3)=rCMFT 
S0UT(IQ,4)  =  ECf1ET 
SOUT  (I  Q  , 6 )  =  0  .0 
SOUTd  Q  ,  5)  =  0  .? 
SOUTd  Q, 11)  =C.O 
SOUT  (I  0 , 1  2 )  =  CO 

fo<  ia,  i  )=’). 

EOCIQ,  2 ) =0. 

E0( IQ, 3) =0. 
EOCIQ, 4) =0. 

FC  (10,  5)  =  3. 

FO(IQ,6)=0. 
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EO (IQ# 7 >=C. 

EOU^8  )  =r\ 

E  G  (  I  U#  9  )  =E  C  M  E  7  /  1 .)  It). 

EQ(IG#10)=TAIR 
EO  <10.11  )  =VL 
F  0( IQ#  1 2  )=H  H 
EC(10,1  1)=W0RK/1:’!'JC. 

C  T  0  L  0  =  T  (  1#1  ) 

1M  =  1 

DO  3 ')  N  =  1#N0P 
I  F  (  AS*  Ft  («)  >  GO  TO  VO 
L’-IND(M)  +  1  * 

LSTOP=INn<M) 

I  M  =  I  M  ♦  1 

TOLD  (I  M)=T(M,L-1) 

AT  AR  =  T  A  R (M) -?73  .16 

<H  C  - 1 1 A ( M  )  *  C TARP-TAR(M) ) 

EO  ( 13#  1  )  =  E0  (  IQ#  1 )  ♦EOT  CM  ) /I  DOR.  001 

ro  do#  ?  >=eo  < in #?>  +  £■'< or  (m#/>  /  ioog.ooi 

E0( IQ#  3) =E0< 10,3) *£ MGR  <  M#1 ) / 10  00. 0C1 
EO ( IQ#  4  )  =C0 (  IQ#  4> ♦£  MGR <M#2> /  1000. °01 
EO  ( IQ#  S>  =E0  (  IQ#  S)  +E  NCR  <M#3)  /  13  90.901 
E0<  1(3#  6  )  =  E0  (  lu,  6)4  f .MGR  (M#4>  /  1000.001 
EC  (IQ#  n  =  E0  (  IQ#  7)  ♦ENGR  <  M  #5)  /  10  OC  .  001 
F0<  IQ#  8  )  =  E0  <  I«#8)  +  LNGR  <  «#6>./  10  00.001 
AT  (IM# 1  1  )  =  T  <M#l 1-273. 16 
AT (IM#  2 )=T(M#1) -273. 16 
!  SK  =  ISK  (M) 

S0UT<I  Q#6)  =  S0UT  (IQ#6)  +  (  T(M#l  )-2/3.  16)*VT(M#L  SIO  /  V» 

EWE  T=E  WFT+E  E  (I  M  )  /E  M  A  X  ( I  M  )  *  A  C  (0)  /TOTA/1  .3000  3  IE 04 
SO  1.1  T  ( I  Q  #  4  )  =•  S  OUT  ( I  G  #  4  )  ♦  F  E  <1  M  ) 

TPLF -0 . 0 
SKHF*0  .  0 
,1  =  1 

CALL  OEFKG(M#IMIN#MAX#LKAr#NOnE) 

DO  SOI  N  =•  1  #  L  STOP 
IF  (N.3  T.  IMAX  (J)  )  J  -  J  ♦  1 
AQTEN=F<M(M#N  )**  (OTFN/QTOL) 

VS-VE3(M#J)*V(M#iV)/VT(M#J) 

CM  =  CME  T  (M#J  >  *V<M#N  )  /VT  (M#J  ) 

IF(J.VE.LSK)  GO  TO  A 
SKDF=SKPF*VS  *AGTEfJ 
A  SO!JT  (I  Q  #  3)  =  S  dIT  (I  G  #  3) +C  M*PM  (M#N  ) 

831  TBLF=TPl  F+VS  *AGTEN 

SOUT(IQ,T1)  =  SOUT(JU#11)*SKnF  /6  0000  .01 
FOOL=FCOLI><  I  N>  *  SK  I N  R  <  I  M  > 

FOT=rWARM(IM)*SKINR( IM) 

WRI TEC  3  #43)  A  NAME  (M#1)#ANAMfc(M#2)#ANAME(M#3)#AT(IM#?),A 
*T<  >M#1  1  )#EPT  (tf>  #I  NGR  (M#  1  )#EiMGR  CM# 2  )  #ENGK  (h#5  )#ENGR(  M#4 
*)#FNGR  (M#5)  #FNGR  (M#6)#StcnF#  TPL  F  ,QCC#AT  AR  #  FLU  X<3M)  #F  OT# 
*  F  0  0  L 

A3  FORMAT (  1X#3AA#0P?F  ?.?#-3PlOF  7.2#OP1F7„?#OP3F  7.2) 

Jj  CONTINUE 

IF(IDA.EQ.O)  60  TO  36 
WRITEC  3.3?) 
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c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


5?  FORMAT  <  MX#  1  25(  *X' >/T30# '**  WARNING  ** ' / T?  C,  ’BOUND*# 
**RY  CONDITIONS  HAVE  UtF.N  CM  A  FT.  E  D  '  /  7X  ,  0  EL  t  M  EPS(M)  °  , 

*2X# aBC <M#11)  RC(M,8)  BC  (  M,  1  0)  *  /  ) 

DO  38  M  =  ?,NOf> 

IF(ASYFL(M)>  GO  TO  38 
BCD UM  =  B  C CM# 3  )~?73 .16 

W  PI  T  E  <  3,39)  AN  AMT  (M  ,  1  )  #  A  N  AMF.  (  W,  2  )  ,  A  N  AME  CM  #  3  )  #  EPS  (  M  )  ,  BC  D 
*UM,BC(M,11),BC(M,8},BC  (  M  ,  1  0  ) 

39  F0RMAT(1X,3A4,2X,F5.4,4X,F5.2,5X,F8  .4,3*,F6.4,4X,F8.4> 
3b  CONTINUE 


36  E  C  ( I Q,  5  )  =  F.O  {IQ,5)  +  E0  ( I  Q  ,  9) 

GOUT  (I  Q  ,  4 )  =  S OUT ( 1 0 ,4 ) / T OTA / 1  000  .00 1 
SOUTCI  Q#3)  =  SCUT(I«,5)/TOTA/1C0C.001 
SOUT  (I  0  #  ? ) *: f  C (  T  6  #  3  }  /  TOT  A 

SO UT (I  Q ,13) = ($0 UT ( 10,3) -ECME T/T OTA/1 00 O.Q01- SOUT (IQ »  2  ) 
*)/(TARP-?73. 16-S0UT(10,6)) 

SOUT (IQ ,5)=E 0(1 0,2)* 1000.001 /RHOCPV 

SOUT (I  Q ,10)  =  T(?#2) -2  73 . 16 

$0UT(IG,9)  =  T(2,?>-27  3.16 

SOUT (I  Q ,fc)  =  TARF-?73. 16 

SOUT  (I  Q,  7)  =  T  (1,2)  -  273. 1  6 

SOUT  (I  Q  ,12)  =  SOUTC  I  ft,  1  2)  /  600  00. 01 

SO'JT  (I  CS,14)  =  1C0  .  *  F.  W  F  T 

GO  TO  90 


32  K  =  60 
L=  1 

LX  =  I Q-  1 

77  WRITE! 3,75) 

I F (LX. IT  .AO)  K=K-6C*LX 

75  FORMAT  ('l',*  TIME  HFLOW  HP  tV  TP  TS 

*'  TH  TO  TR  TM  S?1  CO  COND  PWET') 

WR1TFC  3,76)  (  (SOUT  U  I  ,J  J  )  ,J  J  =  1,14),  I  1=1  ,“■ ) 

7  6  FORMAT  (  1  X  ,  F  5  .1#  3F  ?.2,6r  6.2,  7\7  . 3,2?  7.2  ) 
irax.LC.4Q)  GO  TO 
LX  =  L  X- 6  C 
L -L*  63 
K=K*6D 
GO  TO  77 
7  S  K  =60 
L=  1 
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LX  = IQ-  1 

37  WRITE* 3*85) 

TFILX.LT. 60)  K  -  K  -  6  '}♦  L  X 

85  FORMAT  (  *1 ’*•  ALL  UNITS  IN  KC'L/HR  EXCEPT  TIME  =  MIN.*’* 
T  A  f  R  =  DE  G  -C  *  VEL  =  M/SEC/-  AND  RH  IS  NON-DIMENSION’* 

*  ’  A  L  1  /  /  1  X  *  *  T  IMF  INBAL  BF.  C  RBFC  BLHTR  ’* 

*•  MET  C  0  M  V  RAD  EVAP  RfcSP  T  A  1  *  * 

*  *  R  V  EL  RH  WORK*) 

WRI  TF.<  3*86)  <  SOUT(1 1*1  )*  <E0(  J  I*  J  J>*J  J  =  1  *13)  *1  I  =  L*K) 

36  F OR F  AT ( 1 X , F  5  .1*13F9.3) 

IF(LX.LE.60)  G.0  TO  88 
L  X  =  L  X -  60 
L  =L  «  63 
K  ~K  +60 
GO  TO  87 

38  CONTINUE 
STOP 

END 


c 

<> 

<> 

c 

<> 

* 

*** 

kk  k  k  it 

<> 

c 

<> 

•*  * 

* 

* 

k 

it  k  k  k 

<> 

C 

<> 

**  *  * 

* 

k 

Ct 

it  k  kk 

<> 

c 

<> 

*  * 

k 

it 

k 

*  k  k  k 

o 

c 

<> 

**  *  k 

★  it  k  *  v*r 

kk  k 

k  .V  k  k  k 

<> 

c 

<> 

<> 

c 

<> 

k  *  *  * 

it 

k 

kk  k  k  k 

<> 

c 

<> 

k 

1c 

it  k 

k  k  k 

<  > 

c 

<> 

k 

it 

k  k 

k  k  k 

<  > 

c 

<> 

k 

★ 

W  k  k  kk 

k  kk  kkk 

<> 

c 

<> 

*  *  *  * 

*  k 

k  k  k 

<  > 

c 

<> 

$> 

ooooo 


c 

BLOCK  DATA 

IMPLICIT  RE AL*8 (A-H*0-Z  ) 

DIMENS  ION  BC(34*11),0R(34*9)*R<34*21>,\/'34*2n>,VT<34,9 
*)*At  (34  )*KHO  (34  *9)  *C  (34  *9)*  AKt  34*9  )*BM  (34*20  )*CM£  T<  34* 
*9)*VE3 ( 34 *9 >  *TA  P( 34) „ A( 34*2  0 ),B  <3  4  *20)  ,D  <34* 20}  *E (3  4*  2 
*0  )  *  T  (3  4 , 21 )  *  ID(126)*TVEN(34>  *IJ  A  (34  )  *VE  LB  C  3 4 )  *1  A  (  1  26  )  *  E 
*NGR  <  34  *7>*FFiT<34>  ,QB  (34  *9),P  F3  (  34, 9)  *T  V  <34  )  *  WARM  (34  }*  C 
‘OLD  I  34  )  *  FLU  X  (34  )*  ERROR  <  34)  *SKINR  <  3  4  )  *S  KI  NS  (3  4)*SKINV(3 
*4)*SKI NC  <34 ) *WORKM<  3* ) *  FHILM (3  4)*Fr (34 )*F  SET  <34 ) *PSKIN 
*(34)*T0(34*21}*RATE<34)  *T0LD<34)*EF(34)*EMAX  <34)*ANAMF 

*  <34  *3)  *  TSK  <  3  4)  *  IMU(34  }  *  FV.ARM  (3  4  )*F  COLD  (3  4  ) 

COMMON  /  FNV/P  C/L  DR  /  DR/LK  /R/LV/V  /LVT  /  VT/LAC/  AC  /LRHO/R  HO  / 
*LC  /  C  /L  AK/AK / LPM /BM/LCME  T/CMF  T/LVHP/VFB/LTAR/ TAR/LBL  D/R 
‘HOB  *C3  *  (ARP  /  CO  A  /  A  /  COP /B /COD  /  t  /  COF  /  F  /  w  T  /  T  /L  I  D  /NUM*  I  S  PH  E 
‘R  *NUMS  *  ISS*N0P,  ID  /  LTVFN  /  TVE  N/LUA/UA/LVELB  /VE  l  B/LI  A  /  IA/ 
*L  rNGR/ENGR/l  F_BT  /  F  B  T  /  E  QB  /OB/  L  BF  B  /BF  0/ L  T  I  M  F  /  CO  *E  T  I  ME  *  F  W  R 

*  I  T  E  /  LT  V/TV/LTCT/TCTV*TOTA/LVARM/WARM/LCOL  D/C  ^LD/LF/F/L 
‘ERROR/ ERROR /LSK INR /SK IN  R/LSK IN  S /SK INS/LSK INV  /SK INV/ LSK 

*  I NC  /  $<  IN  C/L  WORK  M/ WORK  M/L  CHI  L  M/  CHI  LM/LE  B/EB/L  HR /MR/LHC/ 
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*HC/LPSK  IN/PSKIN/LTO/TO/LRMTE  /RATE/LT0LD/T0L5  /LEE  /EE  /LE 
*MAX/EM  AX  /LIMl'/IM'l/L  FLUX  /  FLU  V  /»  T OL  W  X  /  T A  I  P #  VL  #  RH #  WORK  #rc 
^FT/CTOLD/L  NAME/ANAMF/LISK/ISK  /LSTOl  VC  SH#  SSA  #PSW#CD  IL# 
*SDIL#Pi>  IL#C  C  ON#  S  CON  /  PC  ON#CCHIL  .SCK1L»P'H  IL  #3  ULL  -ABL  #**€ 
*0/LQTEN/9TEN/l  F  SFT/F  F6  7  /LFWARM/FWARM/L  FCO  L  0/  FC3LD/LOTB 
*L / RTBL  /LST9P/ST0PIT/LM  US/M  US 
LOGICAL  ASy  *'L 

DATA  CSW#SSW#P$W.COIL#$DIL#PDIL/CCON#SCON#PCON#CCHIL#S 
*CHIL#PCH  IL  #  F  UL  L  #ABL#PCO  / 12  *  n  .  J  #  ?*  1  OoODQD .  #1  .  0/ #  QTFN  /ID 

*  JOG  .10.  /#QTBL  /1C000GU./#  STOP  I  T/  .  001  /  #  U  71  4  *30  9./ #  I A  /  1 # 2 
*#4#3#o#5#8#  7 fin, 9, 116*G/#BC/374*n.f)/#EN6R /23  f>*0 . 0/#  EB  T 
*/34*0.  0/#NUM#ISPHFR#  ISS  #ID/3  4#  1  ,1  #  10*1  11 1  #11  6*0/#C0  #FT 

*  I  ME# EH  RITE/ 3 (0000 .#2*0. 0/#TA  R/ 34*309 .936 37/ # RHCB#C« /I . 
*08#  0.9  2/  #  TAR  P/M  9 .93637/#T0/  714*309.0/  #TA I R#  VL#RH#HORK 
*/29.45#  0.10#Q.3#0.0/#RA  TF/  3  4*0.0/#  UA/  34*0./#  SKIN®/ 3  4*0 
*.0/#SKIN$/34*0.P/#SKINV/34*O.f)/,SKINC/34*0.0/#W0RKM/34 
**C.O/#  CPILM/34*  0.0/#FLUX/34*0.  C/#  FSET/34*  0.0  /#RFR/1  36* 
*0 . / #Q3 /I  36*0  ./#BM/63Q*1  .0/ 

END 


SUBROUTINE  S CON T 
IMPLICIT  REAL*8 (A-H,o-2  ) 

DIFFUSION  TS  (11  )#P(1  1) 

D I  MENS  I  ONPC  (1)#VT(1)#AC(1)#3M(1>#CNET(1)#VE3(1)#T(1)#1 
*0(1 ) #E  NGR<1 >  #GB (1 >#9FB( 1)#I AU( 1 >#AK (1 >#R ( 1 )#  FSET (1)  #F W 

*  A  R  M '  1)  #  FCOLD  (1)  #W A  RM (  1)  ,COLb  (1  )  #FL  UX(  1  >#F.  RR3  R(  1  >  #SK  INF! 

*  (1  )  #  SK  I  NS  (1  )  #SKINV(1  )#SKINC  (  1)  #  WO  R  KM  Cl  )#CH  IL  M<1  )  #EB  (1  > 
*#PSKIN  (1  )  #T0  (1  >  #RATE(1  )  #TOL  D  (1  )  #  EE  ( 1  >  #  EM  AX  ( 1  )#ANAMF  (1) 
*#  I  SK  (1  ) 

COMMON /FNV/BC/L  VT/VT/LA C/AC /LBM/BM/ LC MET /CMS  T/LVFP/VFP 

*  /IT/T/LAK/AK  /LID/NU'1#  I  S  PHE  R  #  NUM  S#  I  SS#NOP#  10/  LQ9/0B/LPF 
*P/PFB/L  IMU/  If’U/LTIME/CO#ETIME#EWRI  TE/LR/R/tH  ARM/WARM/t 
*COLO/COLO/L  F  /F / LE  RROR /E  RROR / LSK IN  R /SK l NF / LSK  INS /SKI  NS/ 
*LSKINV  /  SKIN  V  /LSK  INC /SKI  NC/LWORKM /WORK  M/L  CHI  L  P/C  HI  LM  /LC 
*B/EB/L  P SKIN  /  PSK  IN/LFLUX  /  FLUX  /LTO /TO /LRATE/RA  TE/LTOL  0/T 

*OLD/LE  E/EE/LEMAX/EMAX/S  TOL WX  /T  A  I  R#  V L#R  n#WO  R<  #F.C ME  T#  C  TO 
*L  P/LNAME/ ANAMF/L1  SK/I  SK  /LSTOL/CSW#  SSW#  PSW#C0  IL  #SD1I.  #PD 
*IL#CCON#SCON#PCON#CCHIL  #SCH  I  L#PCHIL  #I3UI  L  #A3L  #PCO/LO  Tt  N 
*/QTFN/L  FSET/  FSET/LFH  ARM/FWARM/LFCOl  0/FC3L0/-QT9L  /QTPL  / 
*LMIJS/AUS 
LOGICAL  ASYFL 

DATA  TS/0.0#S.,1P.#15.#00.#?5.,30.  #  35  .  #4  0  .  #4  5.  #  50.  / 
DATA  P/ 4. 579  #6. 543  #9. 2  09  #12 . 78  8  #1  7 . 535  #2  3 . 75  6#3 1  .824,4 
*2. 17  5# 5 5.3?4 #71  .88 #92  .51/ 
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DATA  17/?/ 

R  AT  E  S=  0 
AGO  =  1 . 

strn-o .  G 

WARMS*  O’. P 
F WAR MS =0.0 
COLDS*  0.0 
F  COLDS  =0 .0 
CW  ARM  =  0 . 0 
C  C  OLD*  0  .0 

CRATE* (T (1 ) -cfoLP ) /9C (1 37) 

1  FCWORK  .61.7173?.  )  30  TO  105 
WORK  1=  0  .0 
GO  TO  1 r6 

105  WORK  1=  (WORK-7173?.  )  *fi.  7  8 

106  IM-0 

501  DO  502  M  =1 / N  CP 

IF (ASY F LCM) )  GO  TO  50? 

1*1=  IM+  1 

CALL  DEFSRG(M,LMAT,NODF/IZ) 

L  =  M  ♦  NO  M  *  (NODF-1  ) 

LI  F\|  UM*  (NODE-2) 

FLIIX(IM)  =  (T(LL>-T(L))*AK(MFNUW*(LMAT-1)>/CR(LL)-R(L)> 
RATE  Cl  M  )  =  (T  <L)  -TOLD  C  IM  )  )  /3C  C’HNUM*  4  ) 

RATF  S=  RATES*  RAT  T (IM)*SK INR ( I M) 

00  505  N  =  1#NQDE 
JN  =  M  +  X  UM* (  N  -  1 ) 

306  BM(JN)  *2.**  C  (T(Jn)-TO(JN))/QTEN) 

WA  RM (I  K ) =0. 0 
n •'/  RM  C  IM)=0 . 1, 

COLD  (I  M)=t  .0 
FCOLDC  I  M  )  =  0  .  .0 
MSK=M*  NIIM*  (  NODE  -1  ) 

31  C  ERPORC  IM)=T CMS* ) -TGCVSK  ) 

F ERROR =FLUX (IM)-FSkT (IM ) 

IF(ERRORClM)  .GT.O.J)  GO  TO  5(4 
lTCERRORCIM)  .tO.O.0)  iiO  TO  3  07 
COLD  (I  M)  =-E  R  KOR  ( 1  M ) 

COl  DS=  C  OLDS ♦  COL  D  C  I II  )  *  S  K  I  NR  C  1  M) 

SUMT*S  UMTf( TO(MSK > -2  7  3 . 1  6)  *  S  K I  N  R  C  I  M  ) 

GO  TO  307 

30  4  WARM  (I  It  )*EPRORC  IM) 

UAR"S=  WARD'S*  WARM  C  1  M  )  *  S  K  I  NR  C  1  M) 

30  7  I  F  (FERROR.GT  .0.0  GO  TO  3PF 
IFCFERROR.ro  .0.0  GO  TO  302 
rCOLDC IM  )  =  -F  ERROR 

FCOL  DS  =  F  COLD  SFF  C  0 1  D  (  I  H  )  *  SK  I  N  R  <  I  It ) 

GO  TO  3 r  2 

308  FWARMC I M)=F FRROP 

FWARMS  =  F  WARM  SFF  WARM  (  TM)  *  SKI  MR(IM) 

302  CONTINUE 

TCOLD=COLDS”  SUM  T  f  ?  p, .  5 
CE  RROR  =  T  <1 ) -  TOC  1 ) 

I  F  (TERROR. LT  .0.0  GO  TO  303 


CWARMs  C  ERROR 
F  C  COLD  =  “CE  R  R  OR 

go  to  3  os 
333  CCOLD-- CERROP 
FCCOLD - ~C F R  R  OR 

33  5  SWE  AT  =  C  SW*C W  ARM  *SSW*  FWA  RMS 

D  I  L  A  r=CDIL*CUARM*SDlL*F  WARMS 

SIR  I  C=  r  CON*  C  COL  O+SCONoF  CiU.  0  S  ♦  P  C  CM  *  C  PL  D  S 

IF(TCOLO.LT.G.O)  AGO=0. 

CHI  L  L=  C  CHIL  *  FCC  01.  O  +  SC  MI  L  *FC  OLOS  *PC  H  I L  *  T  F  0  LD  *  AGO 
I  F  (SWE  AT  .GE  .  0.0)  GO  TO  312 
SWfc A  T= 0.0 

312  IF  (D1LAT.GE  .  C.O)  GO  TO  5H 
D  I  L  A  T=  0  .0 

314  IF(STR  IC.GE  .0.0)  GO  TO  .316 
S  T  R I  C  =  0  . 0 

316  I  F  (CM!  LL.GE  .  C.O)  GO  TO  3 IP 
C  H  I  L  L=  0  .0 
318  I  M  =  0 

430  PO  401  M-1,N0P 

I  F  (  A  SY  F  L  <M)  )  GO  TO  4  1 

TM= I Mf 1 

L=ISK(M) 

cus=n. o 

IF  <M.S  F  .MU?)  CHS  =  <ABL*  COLDS*  PDl  L*f  COLDS  )*  C  till  M<  I  M) 

CAl  l  MUSCLE  (M#IF;IIM#<A/KO/KC  ) 

I  MINJH.Lr.O)  GO  10  r  3 
no  TO  (510*520/531)' 1 NU  M 
530  NC=M-MUM*KC 
WB  =M  *\J  UM*KB 
M  A  =  M  4  \|  II M  *  K  A 

VM-VT(NC)*VT  <Nr)*VT(OA) 

C  v  F  T  (  M  C  )  -  OB  (  \'C)  ♦(WORKMl  1  M)  *  W  OR  K  I  *  C  H  1  L  M  ( I  M  *  C  MlLL)*VT(N 

*  C  )  / VM 

VCB  (  NC  >  -HF  B  (  MC  )  ♦  (  WORKM(  IM)  *  WORK  1  *C  HI  LM  (  I  M  )  *C  F’l  LI  )  *V  T  (N 
*C)/V**PCO-CIIS 

I  F  ( V  EB  (NO  .L  F.O  .005*PFP  (NC)  )  V  F. »  (N  C  )  =  .  P>'5  *  PF  F  <  NC  ) 

GO  TO  505 
52<  NP  =  M  U**Ki: 

N  A  -M  +\|  UM*K  A 

VM  =  VT( N  V  )♦ VT  (LA  ) 

525  CM  F  T  (S  E  )  =OF  (  UP)  *  ( WORKM  (  1M)*W0RKI*CH1LM  (IM)*CH1LL)*VT(N 
*P  )  /  V  M 

Vt  B (N3  ) =BFB ( Nb>  +  (*G  RKM  <  IM) *  WORK  1 +CHI LM  ( 1  M  )  *C  I  I ll ) *  V  T  <  N 
*B )  /V**  PCC-CUS 

IF  (VfR  (Nb).l  r  .  P  .  0  j  5  *  B  F  8  ( PO  )  1  V  F  B  ( MP  )  =  .  nr.  S  *v  F  H  ( UB  ) 

GO  TO  515 
510  N  A  •=  M  +  N  U  M  *  K  A 
VM  =  VT(  NA  ) 

515  C  v  E  F  (V  A  )  "GB  <  LA  )  ♦  (Wl.TKM  (  I M)  *  WORK  I  +  CHILM(IM)*C  H] L  L ) *  V  T  (  N 

*  A  )  /  V  M 

V  F  B  (  N  A  )  =■  P  F  0  (f.A)  +  (  WoR  K  *'  (  [m)*UORKI+Cm1LM  ( I  v  )  *  f  HI  L  L  )  *  V  T  ( N 

*  A )  /  VM*  P  f  0-C  U  S 

IF  (  V  E  3  (NA).L  F.O.005*BFB <  N  A )  )  VFP(NA)=.P05*FFB(NA) 

B~  15 


5  3  3  N  J  -  M  4  M  U  M  *  (  l  -  1  > 

vrn(Nj)-HrM<M)-s7Kif*sMNr<in)4niLA7*r.KiMv(ir) 

J  F  (  vr?  (  NJ  )  .  1. 1 ■  .0  .0i  b*Rf  U  <MJ  )  )  VHi(N  J  *PF  H(N  J  ) 

5  01  ff!  <  !  I«>  =  <FR<  I  M)  4SK  I  NS  <  J  M  )*SWf  AT  W.O**  < »  KR0R<  1M>  /HUL  l  ) 

CALL  0  F  f  SRG  <  I*  ,  L  T  ,  N  0  0  F  <•  1 7  ) 

T  7  ~ I)  C  (  M  +  N I IM  *  7 )  -  ?  7  3.  1  r> 

CALL  TFRCrtf.  (TS/P^TZ/PAI  B  ,  1 1  ) 

P  A  I  K  sP  A  I  R^PC  ( M  *  N U  xi  *  7  ) 

NU-NOi)  F  4  1 

T  SK  *  T  (  M  ♦  MJM  *  (MJ  -  1  )  >  -  ?  / 3  ..  16 
C  *'>  l  L  T  F  RPRG  (TS#P#T3K»PSIUN(1M)#11) 
f  M  A X  <  I  M  )  ~  <F  S KI  N  <  I  M  >  - P  A  1  R  ) *  7  .  7*  P  C < M  > 
n  (f  e<  im;  .lf  .r.AXdi))  goto  4  3 

f  I  (  1M)  —  F  MAX  «  ]M  ) 

4  0  3  PC  (M4N  IIM*  5)  =  FF  <  IM)  *  T  <M4  Mlif**  (  NU-  1)  ) 

4  0  1  CON  T  IV  Ilf 

F.C'IF:  7*  m  73?  .♦WOKKI  /  ./H  4CHI  I.  L)  *0.0  023*  «44  .P-PAI  R> 

CMF  T  <?  )  :|)BI  3 )-[  CMF  T*  .? 

C  M  Z  T  (  3  h  )  =  cm  (  38 )  -  f  C  M  t  T  *  .  1 
CMF  T  (7  1  )  =0P  (  71  )  -  I  CME7  *  .  7 
V  F  B  <  *i  >  =  C  0  H  F  I  (  S  ) 

Vt  P  (6)  =RFP  <  6  )-PSn*  F  NARMS 
RF  TURN 
F  NO 

SUP  ROUT  INF  CPEFFG 
I  VP  l  1C  I  I  RF4l*<MA-H,P-?) 
niPflJSION  I  M  IM  Hi)  ,IM*X  (  10) 

r>  IMF  MS  ION  PR  11)  rR(1)/VT  (  1  )  ✓  T  AR  (  1 )  ✓  R  HO  (  1  )  «  C  (  I  >,AK.(1)/RK(1), 

*.CMf  T  (1  )  ,  Vl.i'  (  1  >,  A(  1  >,u  <  1  )  ,P  (  1  ),  F  (1  )  ,  10  (  1  )  ,  11  ft  (  1  >,PF  B(  1) 

COMMON  /LOR/D  W/  I  R/H/LVT/  VT  /  L  l<  HO  /  RMO  /  l  C  /  C  /  I  AK  /  AK/l  RM/PM 

COMMON  /irm  I  /  CVFT/LVf  P  /  VFFJ/l  TAR/1  AR /LIU  P  /  R  Mf)  P,  C  R  *  T  A  RP  /  C  OA  /  A  /  CP  P  /  R 

COMMON  /COP/ P  /CO  I  /k  /l  IP  /  NL'KI,  I  SP h  r R , N UM S , i S S , N  OP,  I  P 

COMMON  /l RR/UU/L  P  F  R /I F  3 

COMMON  /  L  0  T  k  N  /R  T  F  M 

COMMON  /  l  UTP  l  /RTPl 

LOGICAL  ASYU 

I  CM  -  1 

!  S  P  -  1 

po  <n  i  - 1  ,npp 
I  F  (]  .3  T  .  ISI’MFR)  I  SP=? 

I  C  t  N  -- 1 

1  f  <  ASY  FLU))  GO  TO  9  ’ 

C  A  l  L  D  E  FNGU  ,  l  M  I  N  ,  I  M  A  X  ,  L  MA  T  /  NO  P  t  ) 

K  =  1 

I  r  (  I  MI  N  (K  >  .  L  L  .1  MAX  (K  )  )  K  =? 

GO  TO  1?W4),K 
?  4  icro=? 

R  H  -  A  K  (  I  4  MUM  *  (K-1)  )  /PR<1  4  MUM  *  <K  -1)  >  /  4. 

PO  HI  N  =  ?  *N  0  OF 
M  P  -  1  ♦  N  U  M*  K 
MMs  I  4N  U**  <K  -  1  ) 

MM  -  I  4  N  U M *  ( N  ~  1) 

1  F  (  iMI  N  <K  )  .  I  0.  I  MAX  <  lO  )  GO  TO  14 
I  F  (N  „E  0  .  IMIN  <K>  )  GO  TO  6 
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If  CN.l.  V  .  IMA  X  (K  )  )  GO  TO  S 
4  ULOWnl 

1MN.EQ  .NODE  )  GO  TO  12 
AA-BB 

R  P  =  D  R  (  H  M  )*D  R  (MP) 

BB  =  ( AK  (MM)*DR(MM)4AK(MP)*DR  Owp)  )/RP/WP/?„ 

I  GO*? 

GO  TO  40 

6  I  F  <N .t Q  .2)  GO  TO  16 
AA=BB 

BB “AK ( (MM) /4. 

GO  TO  4  0 

8  AA=BB 

GO  TO  40 
19  I  F  LO  2 
12  A A  =BR 

GO  TO  <46, 47), ISP 

46  BB=AK<MM>/DR <MM) /2. 

GO  TO  48 

47  BBT=R(MN)/DR  (MM ) 

P8  =  AK<  MM)*< BPT+ .5 ) /<? .*SPT+1  .)  /  OR (MB) 

48  GO  TO  <  40,1 7), I  FLOW 

14  I  F (N.EQ  .NOD  E  )  GO  TO  19 
RP=DR(MM)*DR  (HP) 

BB=(AK (MM)* OR (MM) ♦A<(MP )*DR  <MP) )/RP/RP/2„ 

I  G0  =  2 
GO  TO  1  7 

16  BB=AK<MM)/DR <MM> /4. 

17  RP=f)R(MM)*DR  (I  +  NtiH*(K-2)) 

AA=( AK ( MM)*  0R(MM)4AK< I*  NUM* ( K- 2 > )* DR <1 +9UM * ( K-2 ) ) ) / RP /R P/2 . 

4  AQTf N=QM  (MN ) ** ( 0 T tN/OTOL ) 

10  CM  =  CME  T  (MM)  /  VT <  *M) 

VB-VEB ( MM) / V  T ( MM) 

9  GO  TO  (  43,44  ),  ISP 

43  AP=„75*<R<MN)+R(I+NUM*N))**?/(3.*DE'fMM)*R(MV))**2  +  DR(MM)**3) 

AM  =  .75  *  (RCMN  >4fM  I  4N'JM*  <N-?>  >  )**27  (3.*D  R<MM)*  R(MN)  **  24  DR  <MW)**3  ) 
GO  TO  4  5 

4  4  AP”<1 . ♦ R ( I  +  NUM*  N) /R(MN) ) /2. / OR (MM) 

AM  =  (1.  +R(I+NUK*(N-2)  )  7R  (MN)  )  72  ./PR  (MM) 

4  5  GO  TO  ( 27,28  ),ISS 

27  RC=1. 

GO  TO  51 

28  RC=RHO (MM)* C  (MM) 

51  A (MN)- < -RXOB*CP*VB*AO TEN-AA *  AM -RB*  AP ) / RC 
B(MN)  =  AA*AM/RC 
D (NN)=BB*AP/ RC 

E  (MN  )  =  (  CM*OM  (MN)  ♦  RH  00  *  C  R*  Vi)  *  AO  TEN*  TAR<  I;  )  /  RC 
GO  TO  (80,18), I  GO 

18  K=K  tl 
I  GO-  1 

30  CONTINUE 

DUN=pR (  I  ) *  2 . 

GO  TO  <33, 34), ICE  N 
33  AK  BOE=  A  K  (I  ) 
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GO  TO  <49, 50), ISP 

3  4  A  K 13  0  B  =  <  4KU  >  *DR  <  I  >  ♦  AK  <1  *NUM  )  *D  R  <  I +  NUM  )  >  /  <  fiS  <  I  >  ♦  D  R  (I  +MJM  )  ) 

HO  TO  <  45/SO), ISP 
49  A  $  t-  ~  6 . 

r-0  TO  3  5 
5  0  A  S  p  =  4  ♦ 

55  no  TO  <25/26), ISS 

25  RC  =  1. 

GO  TO  52 

26  P  C  =  R HO <  1  )*f  (  I) 

5  2  A  (  I  >  =  <  -  ASP*  AKOOB/l>UN/[>UK-RHOb*CR*VE£i  (I  )/VT  (1  >  **  (OTFN/QTPL  )  )  / 

■  R  C 

n  <  i )  -3 .  o 

D ( I ) =A  S  P  *  A  KB  08/ PC /DUN  20  UN 

f  (!  )  -  (  C  M  F_  T  (I  )/VT(I  )  *  F  M  (  |  )  ♦  R  H  OP*  CH*V  ER<  I)  /  VT<  I)*TAR<I)*BM<T)**(OTEN 
*/QTBL)  )  /  R  C 
ISP  =  1 

90  CONTINUE 
RFTURN 
END 

SUBROUTINE  f  HA  L  (  <  ,  T  T  ) 

IV P LIC I  T  RE AL*8 < A-H/O-Z  ) 

DIMENSION  1  M  I N  < 10) ,  I M  A  X  (10) 

D  I  MENS  I  ON  BC  (1 )  ,l>'-  <1  ),R  (  1>  ,  V  (1  )  ,VT  (  1),AC  (  1  ),  RMO  (1  ),  C  ( 1  )  ,AK  (1  )  , 

1PM  <  1 ),  CM  FT  <1  ),V  FB<1  )  ,TAR  <1  )  ,  ID  ( 1)  ,  rfJGR  <1  )  ,ER  7(1  )  ,T  (  1)  ,  TT(1  )  ,GR  (  1  )  , 
4HKB  <1) 

COMMON  /LNV/PC/LDR/DR/LR/R/LV/V/LVT  /  VT/L'C  /  AC  /LRHO/RHO/LC/C/LAK/AK 
COMMON  /  LEM  /  BM/LCME  T/  CMF  T/LVE8/  VFB/  LTAR  /TAR/l  BLD/RHOP,  CB  /L  T/T 
COMMON  /  L  ID/NUM,  I  SPHER/NUMS  /  N/N  0F>,  I  D/LFNGR/EN  GR/LEBT  /EBT/LBFB/BF  3 
COMMON  /LQT  EN  /Q  TEN /LOT  BL  /QTBL  /LQR/GP 
17  SUM=0. 

ST0R=3  . 

DEC  =0. 

GEN  =  0. 

I  -1 

CALL  DEER  Ci  (K/IM  IN/I  MAX, l  MAT  /  NODE) 

N  A  -  K  +N  U  M  * N 0  D  F 
DO  30  J=1/NQDF 
I  F  <  J  .3  T  .  IMA  X  (I  )  )  I=IH 
MN=K+N  UM*  (J  -  1) 

M  =  K  +  NU  M  * <1~1  ) 

AQTF.N-PM<MN)  **  <QTEN/Q  FPL) 

10  C*-CME T (M)* V <MN ) /VT <M) 

VB  = VE3 <M)*V  <MN) / VT<M) 

16  $  TuR=S  T0R  +  RH0('1>  *  C  <M)  *  V  <MN)  *  <T  <MN)  -TT(MN)  )  73  C(K+NUM  *4  ) 

PEC  =  PF.  C  *RHO  (M)*  C<M  )  *  V(MN  )*<  T  <MN  ) -2  7 3.1  6) 

E N (> R  < <  +  NUM*6  )  =R  F  C 
ENGR  <<  )  -  STO  R 

15  5HM  =  SUM  +  PH0B  *CB  *  VB*  <  T(MN)-T  A  R<  K  )1*  AOTF.N 
FNCR  (<  +  NUM)  =  SUM 
3  .  G  t  N  =  GF  N  +  CM  *  R  H  ( M  N  ) 

ENGR<< ♦  N  U  M  *  2  )=GEN 

OU  T  A  -  3  C  <  K  )  *  <  T<NA)-t3C  OC  +  NUM*?  )) 

ENGR«<  *NUM*  3  )  =  0  U  T  A 
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0UTB*3C  <K  +  Ni'M>*  <7(NA>**4-BC'  K*  N VM *  5  >  *  *  4 ) 

EN  6R  (.<  +  NUM*  4  >»Ol.mj 
0UTC--3C  (  X  ♦  i4  U  M»  5  )  /T  <  N  A ) 

ENGR  (<  *NUM*5  >"0UTC 

our*ouTA+oiirn+ourc 
18  GO  TO  !1J,H  )fN 

13  EBT  ( K)  =  OUT-G  EN  +  S UM 
RETURN 

14  EBT(K) =OUT-GEN+SUM+STOR 
20  RETURN 

END 

SUBROUTINE  RLTPEG 
IMPLICIT  REAL*  8  (A-H*0-Z  > 

DIME  NS  ION  I  MINI  10)..  I  MAX  (10)  ,DtLTC  126? 

DIMENSION  V(  1)*>VT(1)  »HM(1)*9FB(  1) 

DI  MENS  I  ON  VE  UC1  )  ,  T  AR  M  )  /  T<1  )  ,1  D  (1  3  y  1  VEKi(  1  >  ,U  AID  .VE  LB  (  1  ) 

COMMON /L  V/V /LVT /VT /LPFB/PEE 

COMMON  /LVl.B  i  VEB  /L  TAP  /TA  R /LB  L  D /  R  HOB  , C8  a  T  A  R  P  / 1  T/T/LBM  /BM 

COMMON  /L  ID  /NUM.- 1  SPUE  R  /  M  UMS  /  I  SS  ,NOP  ,  I  D /  l  1  V  EN /  TVF.  IJ /  LU  A/  V A /L  V  ELO  /  V T LB 

COMMON  /  LTIME  /CO^tTIME,r  W»ITE 

COMMON  /l.DTEN  /QT  EN 

COMMON  /LQTRL  /CT 

LOGIC*  t.  AST  EL 

SUM 8=3 . 

SUVC  =D . 

DO  5  I  =  1  <»N0 P 
I  F(  ASY  Ft  (I  )  )  GO  TO  5 

suma~:,  . 

V  E  L  9  (I  )  =0. 

CALL  DEFRGI I#IM JN,IMAX,LMAT,NODE> 

J=1 

DO  1  5  K  =  1  #N 0  OF 
MAT*  I*  NUM*  (  J  -1  ) 

KN=I+NUN*  (K-1> 

AQTEN  -  BM(tCN  )  **<  GTEiN/RTBL  ) 

«  VB-VE3  C  *’AT)  *  V(KN)  /  VT  (MAT) 

10  V  ELO  <  I  )  =VELB  <1  >  +V0*AQTF  N 
25  SUMA=SUMA+T < KN) *VB*AOTEN 
1  F  (K  .E  si  .  IMAX  (J  )  )  J-J+1 

15  CONTINUE 

EVEN  (I  )  "SUM  A  /VE  l.B  <  i) 

5  CONTINUE 

DO  45  1=1, NOP 
I F  <  ASY  F  L  ( I ) )  GO  TO  45 
55  PELT <1 ) =TARP-TAR( 1 ) 

TVEN  (I  )  —  T  VE  M  (I)  +DELT  (I  ) 

40  SUMC=SUMC  +  VELB(I)*rvFN(  I) 

SUMB-S  U^B  +  VE  LB  <  1) 

45  CONTINUE 
CO  =  SSjMB-CvT 
TARP=S  UMC/SUMB 
DO  5  0  I~1,N0P 
I  F ( AST  FL  ( I  )  >  GO  TO  50 
5  5  A L A M  =11 A  (  D/VELU  (I)  /  k  H  /  CB 
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TAR  <I)  =  <TARP  *AL  AM*TVCN  (  I  ))/  (  AL  AM-H  .  ) 

5  1  CONTINUE 
A ?0  Re  TURN 

END 

SUBROUTINE  STEPRGI  TT,DD) 

IMPLICIT  REAL*®  (  A~H,0-Z  ) 

DIMENSION  X  A  7  (  2  1  ) 

D  X  ME  N$  I  ON  OP  (1  )  ,7ARC1)  ,  E:iT(  1  >,  EMGP  <  1),  ANA  ME  (  1) 

DIMENSION  A<1)/P(1)^(1)/E(1)/9C(1IeT(l)/3D(1)/TT(1)/AT(1) 

COMMON  / L  TAR  /  TAR  /  IE  BT  / E  B  T /LF.  N  GR  / ENGR 
COMMON  /  ENV/13  C/COAT  A /  COS  /  P/  C  0  D  /  D/COE  /  F/  L7  /  T 

COMMON  /L  TJMF  /CO, ET  IMF, E  WRITE  /LNAME/aNAM.E/LID  /NUM, I S  PH E  R,N UMS, X  S  S, N 
!iOP*J  D 

COMMON /LSTOP/STOPI T 
1A*T 

LOGICAL  ASYFL 

IT-1 

IL*^ 

CALL  SCONT 
53  12*2 

DO  90  K  =  1,N0P 
YRC-RC  (  1  +NUM  *6 ) 

IF  (ASYFL(K)  )  GO  TC  9" 

DO  (K  )=  TAP(K  ) 

CALL  DEF  SRf>(K,LMAT  /  NODE  ,17) 

ON  -  K  +N  UM  *N0  D  E 
Nl'=NOD  E  *  1 
12  DO  30  1=1, MU 
N-K+NUM*  (1-1  ) 

X  A  T  (  T )  =T(N) 

GO  TO  (31,30), IT 
31  TT(N)=T<N) 

30  CONTINUE 

5  2  IFCISS  .NF.1)  GO  TO  20 

T<K)=-D(K)/A(K)  *T (K*NUM  )-F.  (  K  )/A  (K) 

TT  ( K  )=T  (K) 

DO  A  I =2, MODE 
N  =K  ♦  NU  M  *  (I-  1  ) 

T(N>  =-F>  (N)  /  A  (N  )  *T  (K  +  NUM*  (I  -  ?  ))  ~D(N  )  /  ACN)  *  T  (K  *NUM*I )  -F  (N  )/  A  (N> 

6  r  T  (  N  >  =  T  (  N  ) 

T(NN)  =  8C(K  +  NUM*t)*<;U(K)MT(NN>-UC<K+NUM*2>)+BC(K+NUM)*(TlNN>**4- 
?PC  (K+NUM*3)**4>  *F?C(K  +  MUM*5)  /  T(NN)  )  ♦T(K*NUM*(  NODE-1)  ) 

GO  TO  41 

20  T  (K)*<  -  D  (K  >  *L-C  (K  >NUM*4  )  *T(K  +NUM  >-F*C  <K+NUM*4)  *F(K)-T  T(i()  >  /  (MK)  *3C  ( 

*K+NUM*4 )-1 . ) 

DO  40  1=2, NODE 
N  =  K ♦ NUM  *  (  I-  1  ) 

4  0  T  IN)  =(  -3CCK  ♦Ml|M*4)*(e  (N  )  *T  (K  +NUM*  (  I-?)  )*D  (*  )  *T(K  +  NIIM*I  )  )-BC  IK  +  MUM  * 
$4)*  F(N ) -TT(N )) / ( rC<<  +  NU M *4 ) *  A(N  )-1 . ) 

T  (\N>  =  -PC(K  ♦  Nil  M*y)*(PC(K)*(  B  C(K  ♦NUM*?)  -TT(NN  ))*PCCK+NUM)*(PC(K*NUM 
<*3)**4-TT(NU)**4)-BC(K  +  NUM*5)/TT(NN))  +  T(K+NUM*(NODf-1)) 

41  DO  51  J  =  1  ,N  II 
N  =  X ♦ NU  M  *  (J-1  ) 

TFST=D APS(X AT(J ) - T ( N ) ) 
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X  AT  <  J  )  =  T  (N) 

IMTEST  .GE.XHC)  GO  TO  52 
51  CONTINUE 
9  0  CONTINUE 
I  T  =  2 
I  L  =  I L  + 1 

35  CALL  3LTPRG 
I  AD  =  1 
ENERGY  “0. 

F.  f*  A  T  =  r> . 

DO  73  1=1, NOP  , 

IMASYFL(l))  GO  TO  7? 

ETEST=DABS(DC(1)-TAR<1) ) 

0D(I)=TAR(I ) 

I F  (ETEST.GE .PC(I+NUM*6) )  IAD=2 
CALL  EBAL(I,TT) 

ENERGY  =  FN  E R  G  Y  +  E  B  T (I ) 

?MAT=EM  AT*ENGR<  I  +NUM*?) 

73  CONTINUE 
CALL  S  CONT 
CALL  COEFRG 

GO  TO  (  74,53  >,IAD 

74  I  F  <  3C<  1  +  NUM*6)  .  LE  .  3T0PI  T)  GO  TO  77 
DO  75  1=1, NOP 

75  PC (I +  NUM*6)  =  EC( I +  NUM  *4 ) MG. 

GO  TO  53 

77  EMAT=.05*EMAT 

I  F(DA3S  (ENERGY  >  .LT  .EWAT  )  GO  TO  70 
7P  WR1 TF< 3,79) 

7V  FORMAT  <  /IX,  '  CAUTION  THE  ENERGY  BALANCE  HAS  NOT  NT  T  CRITERIA') 

?•:  IF(ISS.NE.I)  GO  TO  9  5 
WRITE(3,7?)  IL 
DO  IT  K  =  1 , N  0  P 
IF (ASY  FL(K> )  GO  TO  10 
CALL  DE FSRG (K,LMAT,N0DE ,IZ  ) 

NU=NOD  E  +1 
DO  11  1=1,  Nl) 

N=K  +  NUMMI-1  ) 

11  TT <N)*T (M)-2 73. 16 

WRITF.(3,71)ANAMf<K),ANANF<K+NUN),ANANE(K  +  NUN*2>,(TT(K  +  NUM*(l-1)),I 
$=1 , NU) 

1  CONTINUE 

71  FORMAT ( IX, 3 A  4,1 V,2(T15, 1 1F  10.5/)) 

72  FORMAT  </1X,  '  STFADY  STATE  CONDITIONS  HAVE  BIEN  REACHED, . THE 

SNUM&ER  OF  ITERATIONS  FXFCUTFD  =  =>  1  ,  1 1  0  /  /  1  X  ,  '  TH  E  STEADY  STATF  TF"PE 
JRATURE  DISTRIBUTION  F OL LOWS ,  ' / 1 X , ' S T AR T I N G  FROM  THE  CFNTER  NODE  TO 
5  THE  SURFACF  NODE  ==>'/) 

I  SS  =  ? 

95  RETURN 
END 

SUP  SOU  TINE  G  LOM  RG ( GE  0 ,N  UT, N  >P) 

IMPLICIT  REAL*8 l A-H,0-7  ) 

DI PENS  ION  GF  D(NIIT,  N39)  ,  1MIN  (  n  )  #  IMAX  ( 1l'),D  (1  0) 

DIMENSION  DR  <1  )  ,R  (  1  )  ,V<  1  ),VT  <1  )  ,AC  (  1),  ID(1  ),  TVU  ),ANAME  (1  ) 
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COMMON /LDR/DK/LR/R/LV/V  /LVT /  VT /LAC /AC/LID /NUM,1SPHE R,NUMS, 1 SS,NOP, 
SID 

COMMON / L  TV/ T  V/L  TOT/ TOTV, TOT  A /L NAME / ANA  Mr 
LOGICAL  ASYFL 
DMA  PI /3. 14159?/ 

TCTV=D . 

TO  T  A  =3  .0 
W  R I T  E  ( 3,6 8) 

6  8  FORMAT (' 1 ', T ?6, ' SKIN  SUB  C  UT  DIA  DIA  DIA  DIA 

JDI.A  DIA  DIA  SPMF R E" /T7, ‘ELEMENT  C1RCUM  FOLD  FOLD 

$  RATIO  RATIO*  RATIO  RATIO  RATIO  RATIO  RATIO  ANGCRA 

tD)  SEGMENT  ‘/T18,  •  (CM)  THICK  THICK  MATL  MATL  MATL 

1  MATL  MATL  MATL  MATL  LENGTH  ANGLE '/ T27,  *<CM) 

S(CM)  3- A  A -5  5-6  6-7  7-8  8-9  9-10 

$  <  CM  )  (RAP)’) 

NUMS=NOP 
DO  31  1=1, NOP 
I SEG=2 

I F  <  A  SY  F  L  (I)  )  GO  TO  3P 
TV(I )«0 .0 

RE  A  D  (1  ,61)  <  GFQ(I,J),J=1  ,12) 

61  FORMAT(F10.6,1X,F6.5,1X,F6.4,7F5.3,F8.4,F10.7) 

WPI  T  b  <  3,62)  A  NAME  <  I  ), AN AMS (I  +NUM  ),ANAME  <1  +  NUM  *2)  ,  (GE  0<  I  ,  J)  ,J=1  ,  1?) 

6  ?  FORMAT  <  1X,3A  A,1X,F  1  ").6,  1X,F  6.5,1X,F  6.5,1X  ,FF  8.4„1X,  F10.  7) 

CALL  OEFRGCI , T M I N , I M * X , LMAT , NO DF ) 

NU*NOD E  +  1 
t)A=6E0  <  I  ,1  )  /PI 
D ( 3) =D  A -GEO (1,3) 

K=LMAT  +1 

DO  15  J  = A  ,  1 0 

D  (J  )=D  (  J-1)  *GEO  (I,J) 

N  =K - J ♦  1 

IF  (J  .EQ  .K)  GO  TO  17 

15  DR(I+NUM*(N-1))=(D(J-1)-D(J))/A./DFL0AT<IMAX(N)-IMIN(N)*1) 

17  IE=?*(IMAX<N)-IMIN(N)+1)-1 

DR(I+NUM*(N-1))  =  (D(J-1)-P(J  )  )  /  2  .  /  D  FLO  A  T  ( I  f) 

D°(I  ♦MUM*  (LMAT -2  )  )  =  <GF.O  (1,3  >  -GFOU  ,2> )  /A  .  /  DF  I  OFT  (  IM  AX  (  LMAT-1  )-I  MIN 
I  (LMAT-  1  )+1) 

DR(T4NUM*(LMAT-1))=GEO(I,2)/A./DFLOAT(TMAX(LMAT)-IMIMLMAT)+1) 

SFG  =  GE  0(1,12)/?. /PI 
K  c  1 

VT  (I  +NUM*  (K-1)  )  =0. 

R(I)=3. 

R  ( I  +  NJ  M  *  <NU  -  1 )  )  =D  A  /  2  . 

DO  10  N  =  2,N  0  TE 

If  (N.GE  .IMIN(K)  .AND.N.LE.IMAX(K))  DD  =  2.*DR(I  +NUM*  (K-1  )  ) 

I  f  (N.ST  .  IMAX  (K)  )  GO  TO  11 
12  R(  I  +  NUN*  (N-1  ))  =R(  I+Nt'M*  (N-2)  )  +  DD 
I  F  (I  .ST  .  ISPHFR)  GO  TO  3  3 
I  S  F_G*1 

SEG  =  GEO  (1,1  2  )*P  1  *  ( 1  .-DCOSCGE  0(  1,11  )*PI  )) 

V  (  I  4NUMMN-1  ))=  (  (  R  (  1  +  NUM*(N  -1)  )  +  DR  (  I+NUM*  (K-  1)  )  )  **3~<R  (  J  +  NUM*  (  N-1  ) 
$) -D  R ( l ♦ NUM* (K-1)))**3)* SEG/3. 

GO  TO  3  A 
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3  3  V  <  I  *NUM  *  (N-1  ))=pi*((jf(I  ♦  NUM  *  (M -1))+ f>R<  I  ♦MUM*(K-1))>**2~<R<I  +  NUM*  (N 
$>-1  )  )  -6  R  ( I  +  NltK*  (K-1 ))  )  **  ?>*GC  0<  1  ,1 1  )  *SE(> 

3  6  VT  (I  *NUM*  <K-  1)  )  =  VT  (  I  +  NUNMK  -  1)  )  *V  (  I  +N  JM*  <  N  ~  1  )) 

TV<  1  )  =  T  V<  l)  +  V<  1+NUM*  (M-  D) 

(  TO  1  n 
11  K  =  K  +  1 

D  0  =  r>  R  <  I  *NUN*  CK-1)  )+DR(I  ♦  NUM*  (K-2)  ) 

VT  (  I  +  NUN*  (K-1)  )  =0  . 

GO  TO  1  2 

1"  CONTI  VUE  * 

GO  TO  (  ?5,36),I  G  r:  G 

3  6  V  (  1 ) I *DR ( I  )*DR ( I)*GT  0  ( 1/1 1  >*SFG 
AC  (  I  )  =  GEO<I  /  1)*GE0  (I  /  I  1  )  *S  E  G 
GO  TO  3? 

35  VCL)=DR(I)*D!HI>*DR<I)*SE5/3  . 

AC(I)  =  R(I+NUf*NOnE)**?*SEG 
3  7  V  T ( I )  =  V  T ( I )  +  V  <  I  ) 

TV< 1 )=  T V< I )  ♦  V( I  ) 

GO  TO  31 
3  ■  M,'MS=Nl»MS-1 
31  roNTiiu  r 

DO  ? a  N=1,M0P 
IF(ASVFL(N))  GO  TO  70 
CALL  D E FSRG ( N,LMAT,  MODE  ,?.) 

M 1 1  =  N  0  D E  +  1 

WR  I  T  H  <  3  ,6  3)  0  /  AM  AM  F  (  V  )  ,A  UAMF  (  N+NUM)  ,  AN  A  NF.  (  N  +  N  UM*  2  )  ,  A  C  (  N  )  ,T  V  ( N  > 

WRIT£(  3 ,64)  (  VT  (FJ  +  NUM*  (K-1)  )  ,K=1  ,LM  AT) 

WRI TF(  3,65)  <  DR( N+NUM* (K-1) ) , K=1 ,LM AT) 

W  R  I  T  F  (  5  ,  *  6 )  (  P  (  N  +  N  U  M  *  (K  -  1 ) )  ,  K  =1  ,  NU  ) 

VeI TT<  3 ,67) ( V (N+NUM* (K-  1 )) , K  =1 ,*00 F  ) 

T 0 T  V  =  T  o  T  V  ♦  T  V  (N) 

TO  T  A  =T  0  T  A  ♦  A  C  (N  ) 

7  U  CONTINUE 

63  FrRf1AT  (  /  /  T  1  0  ,  ‘DATA  FOR  BODY  ELEMENT’, 16, ?X, 3  A4/T10,  *  S II 1  F  A  C  F  AREA  ( 
iSU . CM)  '  ,2X,  F1Q.  3,T(>5,  •  TOTAL  VOL  I’M  E  (CU 1  I  C  -CM  )  1 , 4X  ,  F  1 0  .  *  /  ) 

64  FORMAT  <3X,* VT  ( CU 0  I C - C M ) * , ? X , < T 20 ,  1  OF  1 (  .  3  )  > 

65  FORMAT  (  3*,*  DR  (  CM)  *  ,  ?  X  ,  ( T  ?  u  ,  10  FIT*. 3)) 

66  FORM  AT  < 1 X, * R ADI  US  (CM)  *  ,2X, ( T?0, 10  FI  0.4)) 

6  7  F  OR  MAT  (  1  X,'  V  01.  U'F  ( C  UR  I  C  -CM )  ’  ,  2  X,  (  T  ?f  ,  1  OF  1  G  .  3)  ) 

TOT  A=T  OTA/10  ICO  .001 
W  P I TF (  3,60) T  OT  V , TO  T  A 

69  FORMAT  (///I  X,,TBTAL  BODY  VOL  l)M  E  =  =  >  *,F1'T.3,V  (CUPIC  CM)  '  ,T0!?,’ TOTAL 
S  BODY  AREA  =  =>  *  ,  F  1  0 . 3  ,  *  (SOU  'RE  METFRS)  '/  ) 

RT TURN 
F  ND 

SUB  ROU  TINE  (>  F  F  R  0  (  L  SF.  C  ,  I  M  IN  ,  I  MA  X  ,LM  A  T,  V  OD  E  ) 

IMPLICIT  REAL* A.  (A-H,0-7  ) 

DIMENSION  iriMUi),  I  MAX  (1  ))  ,1  0  (  1  ?  6  ) 

COMMON  7L  ID /NUM,  1  S  I'll  F  R  ,  N  II,';S  ,  I  SS,NOP,  ID 
17  =  1 

ENTRY  DEFSFG ( L S F C , L M A T , NODE ,  17) 

NO-ID(LSEC) 

NC  DF  =  J 

DO  15  I =1  ,1 j 


B- 23 


NE -NO/  1  P 

I C  C  =  NO “NE  *  1 0 

NO--NF 

I  F  (I  CC  .  FQ  .  •)>  »,o  TO  18 
If*  IN  (I  >  =JCC 
15  NOOE=\!OOE  +  l  C  C 
1  =  1  +  1 
1  $  J -1-1 
i  *?  A  T  =J 

so  to  <2(;,3n>#j? 

3'  RETURN 
20  I STOR* I  MINT  1 ) 
in n (t  >  =i 
I^ftXd  )  -- 1  $  t  o  R 
00  16  K=-2,J 

ISTOR=IMIN(K  ) 

lWINf<  )  =  IM  A  X  (K-1  )  +  1 
16  JNA  X  (.<  )  =  IMl  n  (K  )  +  I  $  TOP-1 
RETURN 
F  NO 


5 

6 


LOFICAl  FUNCTION  ASYf  L  (  1  ) 
COMMON / L  IA/ I  A ( 1 ) 
IF(I.NE.IA(I)>  go  to  5 
ASYF1=  .  FALSE  . 

RF  TURN 

IE(2.3T.IA(I>)  GO  TO  6 

asyfl= .false  . 

RE  TURN 

ASYFL= . TRUE  . 

R  f  TURN 
END 


SUBROUTINE  MUSCLFd 
0 IMENS I  ON  1 M UC 1  ) 
f  C vMON / L  TMU  /  1MU 


f  J  *L*M  ) 


J  =  I  MU  (  I  )-(lMU(l  )/11)*1Q 
K=(  I  MU  (  I  )-(  Imj(  1 )  /nr  )  *  1^0)  /  1 J-7 
L“<IMJ<  I  )-<  mu<  I)  /13P0)  *1C0Q}/1(  Q 
M  =  IMU(  I  )  /  1 0 0 C- 1 


RE  TURN 


1 


END 


OIMENSION  TARLF  (INDEX), COT43 (INDEX) 
IT(6INE.LT.T  ABL  F  ( 1  )  )  (,o  TO  15 
IF  (GINE  .GT. TABLE (INDEX) >  GO  TO  35 
IF  (IMOE  X.GT  .  1G)  GO  TO  AU 
L  =  1 


WANT=COTAB(L  ) 

PROD=1  . 

IND  =  INDFX-1 
GO  TO  Tfi 

4  j  po  1J  j  LA  =  1 ,  INDEX 

LA  =  ILA 

IT(TA0LE(LA)  .GF.GIvE)  GO  TO  4  5 
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13  CONTINUE 

45  I  F  <  TABLE  <LA  >  .EQ  .GIVE  >  GO  TO  60 
I  F  <  l  A  .  t  F  .  1  0  )  I N  D  =  1  0 
IF  (LA.  ST. 10)  I  N  D  =  L  A 
L  =IN  0-9 
WANT  =C  01  AR  (  L  ) 

PR0D=1  . 

1  b  DO',  10  l=L#IND 

PR  OD=PROD*(  GIVE-TAB  LEU  )) 

[)  I '» D 1  F  =  P  . 

I  A- l 41 

DU  30  K  =L  / 1  A 
D  •'  R  i)  D  =  1  . 

DO  25  J-L,IA 

I  F  (J  .EG  .K)  GO  TO  25 

DPRO 0=0  PROD*  (TAEL  E  (K  )-T  ABLE  (  J>  > 

2  5  CONTINUE 

3  DI VD IF  =DI VDI F+C  OT  At (K) / DPR  OP 
1  '  WANT  =WANT+PROD*DIVDIF 

R  E  TURN 

15  WANT=fO T AP( 1 ) 

WR  ITF(  3  #50  GIVE 
RFTCRN 

35  VANT=COTAB(INDE'r) 

V  R  1  T  E  (  3  »  5  '  )  G  IVF 

Vi  format  (1 X, *  THE  VALUE  BEING  INTERPOLATED  (',El0.P,*)  EXCEEDS  TABLE 
T  RANGE/  STANDARD  FIXUP  TAKFN') 

RF  TURN 

6'  WANT  =  C  C  T  AB  (L  A) 

RETURN 

END 
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HEAD 
APIS 
FEET 
NOT  USED 
NOT  USED 
NOT  USED 
65. 676620  .  , 

65. 6  76020  .  t 
49.  1,08470  .  i 
35.  6  '4  71  7  .  I 
60.  8  4  724  9  .< 
79.  1  934,)?  .  ? 
?6. 2  693  3  7  .3 
14.  2071  63  .  4 
35.968834  .6 

72. 6  04  526  .4 
. 70601 

. 265 oi 
.31301 
. 32701 
. 11 "Cl 
. 15801 
.32001 
. 14301  l 
.23901  ! 

.18301  : 

1.05 

1.50 
0.55 
1.9,1 
1 .05 
1  .  5  J 
0.8  5 
1.00 
1.05 
1.50 
1.05 
0.85 
1 . 00 
1.33 
1.05 
0.8  5 
1.09 
0.58 
1.30 
1.05 
0 . 8  5 
1  .f  U 
1.05 
1.30 
1.05 
0.85 
1  .00 
1.70 


F  0  RE  H  E  A  D 
A  P  n  S 
?  rET 
NOT  USED 
NOT  USED 
NOT  USED 


F  A  cn 

HAH  DS 
NOT  USED 
W0T  USED 
NOT  (IS  FD 


4  0  001 

.8160 

.856  .0  .? 

4  0001 

.816'.; 

•  8  56  .0  .0 

4  0  00 1 
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0  .934  6  2 ' 

.883 
.5256 
.640 
.9269 
.PS  3 
.525  6 
.64  ft 
.9269 
.92  69 
.5756 
.9269 
.64  U 
.9769 
.5256 
.9769 
.64  0 
.9769 
.9769 

.5786 

.9769 
.64  0 
.9769 
.88  3 
.5756 
.9769 
.64  T 
.9769 
.52  56 


7  9  .4  5 
79  .45 
79.45 
2  9  .4  5C 
79  .450 
29  .45 
29  .45 
29  .45 
29  .45 


4.536 

lo.nou 

1  .376 
1.79  7 
4  .536 
1 O.OOn 
1  .3  ?  :'■> 
1  .797 
3.61  7 
1  0.000 
3.61  7 
1  .3  7  5 
1  .79  7 
5.0 J  4 
3.61  7 
1  .3  7  6 
1  .79  7 
7.42  r 
5  .00  4 
3.61  ? 

1  .376 
1  .797 
4. 700 
5  .004 
3.61? 
1.376 
1  .797 

1 


13095. 

0. 

0.34 

100. 

719°. 

0. 

0.0  ft 

17. 

80. 

0. 

1  20. 

0.6 

50. 

0. 

34?. 

O. 75 
30. 

1  7  00. 

n 

7631).’ 

P. 59 
44  0. 

36670. 

ft. 

437C.' 
77.7 
6  ?". 

ft. 


NECK 
H  A  V  D  S 
NOT  USED 
NOT  USED 
NOT  UJED 
NOT  US E  o 
.0  .0 
.  0  .0 
-  0  .  c 

.  0  . n 


.0  .0  .0 

.0  .0  .0 

■P  .0  .0 

0  .  ()  .  fj 

.09999999  ' 

.09999999 

.09999999 

.09999999 

..09°99999 

.09909999 

.09999999 

.09999999 

.  "9  999  999 

.09909909 
5.  4  1610. 


THO KA  K 
L  r  6  S 
T  USED 
T  USED 
T  USED 

T  U'ED 

o 


0. 

1  .55 
7960. 
6950. 

0. 

0.26 
494  . 

75  0. 

0. 

600. 

2.71 

2370. 

0. 

10O0. 

1  .16 
3  55  . 

476  0. 

0. 

8560. 

39.01 

1970. 

15 5000. 

0. 

1 407". 
103.  1 
3774  . 

0. 


A3D0HEN 
LFG  S 
NOT  USFD 
NOT  USFD 
NOT  USED 
NOT  USFD 
■  OC  .8568969 
.00  .  14MU3" 

.8475  3.6041049 
.4258  6.23  3185" 
.  5788  6.28  M&5 
.1564  6.28  3185  • 
.4443  6.2831850 
.0224  6.2831850 
■3424  6.2831850 
-3591  6.283185' 
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1.05 

.9269 

3.612  2860. 

9  iGC* . 

0.85 

.64  0 

1  .376  2.58 

11.7 

1.00 

.9269 

1.70  7  4  8  0. 

26C8. 

1 .70 

.5256 

19.60  8  0. 

0. 

1.05 

.9269 

3.617  160. 

53  0. 

0  .85 

,643 

1.376  1 . G  3 
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1.00 

.9269 

1.797  160. 

2216. 

i  .70 
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1  9.60  8  0. 

0. 

1  .05 

,9269 

3  .61  7  4740. 

1  3830 . 

0.8  5 

.64  0 

1.376  11.4 

51  .8 

1 .00 

.9  2  69 

1.79  7  8  50. 

4203. 

1.70 

.52  56 

19.608  0. 

0. 

1 .05 

.9269 

3.61  2  250. 

810. 

0.85 

.64  0 

1  .375  7.43 

11  .03 

1.00 

.9269 

1.797  220. 

• 

1830. 

393  .  66. 

1  84.  200.1  609. 2843. 1637.  603.  3816. 

8  3,9. 

29 . 4501 

G . 1  0.3 

p  n 

9.9 

01 

0324  .0054 

.31  52.01  65 .1  71  1  . 

3724.0686.1845 .1505. 

0  3  34 

03  77. 006  3 

.3  1  77.01  93.166  3. 

314  2.1  54  0.G310.2189. 

0  3  50 

n6 1  5.0103 

.32  89.3  314.  HI  7. 

2 1  0  3 . 0 9  5  0.  1  21 C .  23 0 n . 

i  rno 

1131.0152 

.1  4  14. 9  0  82.0  99  8. 

0008. 27no.  1911 .2104. 

2961 

003  0 .0000 

.3 000. 30  CO. 1 135. 

1865.0890.  ,1100 .6000. 

0100 

000  0  .ooon 

.3  3  39.9436.4414  . 

4  61  1  .00  76.0004  .011  3. 

0007 

0.0 

0.1  0.0 

9.9  0.5 

4 

.0 

0 . 0 

0.141  .050  300.3  1.35 

75 

.0 

9999 . 99 

33.0  10. 

10.  1.5 

23 . 50 Cl 

0.1  0.3 

0.9 

.16666  466 

> 

0 

23 . 5001 

0 . 1  0 . 3 

J  .  I* 

.1  6666  666 

n 

0 

23 . 5  CL  1 

0.1  0.5 

0  .9 

.  1  6666  666 

n 

n 

4 . 7001 

0.1  0.3 

0.9 

.003  33,23  3 

r 

4  .  70°  1 

0.1  9.3 

9  .  o 

.  00,333  33  3 

r-> 

n 

4 . 7001 

0.1  0.3 

8  "i 

.0C 3 S3  333 

n 

4 . 7001 

0.1  0.3 

o‘.r 

.9  03  33  S3 3 

0 

4 . 70C  1 

0.1  0.3 

0  .9 

.90333333 

0 

4 . 7001 

0.1  0.5 

9.9 

.00833  533 

f 

n 

4  .  7'"  91 

9 . 1  9 . 3 

0  . n 

.0  08  33  533 

0 

4 . 7001 

0.1  0.3 

0  .*  •' 

.00833  333 

0 

4 . 7001 

0.1  0.3 

0.9 

.90833333 

0 

4 . 7001 

0.1  0.3 

0.9 

.00833333 

r 

0 

4  .  79  n  1 

9.1  0.3 

r\  m 

.90833  533 

r\ 

4 . 7001 

0.1  0.3 

o.r 

.90833333 

4 . 7GU1 

0.1  0.3 

» 

.00833333 

0 

4 . 70 C 1 

0.1  0.3 

0  .9 

.01666666 

0 

c 

4  .  700 1 

0.1  G  .  3 

0 . 9 

.01666666 

L 

0 

4 . 7301 

n.1  0.3 

9.9 

.  9  1 666  666 

r 

0 

4  .  7 001 

0.1  0.5 

O.o 

.91 666  *66 

0 

4 . 70 LI 

9.1  9.5 

o .  ■: 

.  91666  66  6 

n 

0 

1  .  7001 

0.1  0.5 

0.9 

.  01665  f>66 

7 

0 

4  .  7901 

0.1  0.5 

0." 

.  0 1  666  f>66 

' 

c 

4  .  7^01 

1  n  .  3 

n  _  -> 
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0 

4 . 79  01 

9.1  0.3 

o.n 

.01666666 

f 

n 

4  .  7001 

0.1  0.3 

o .: 

.01666666 

n 

0 

4  .  7001 

0.1  0.3 

0.9 

.04166646 

p 

0 
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4 . 7") Cl 
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0.0 
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4 . 7001 
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P.n 
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0.3 
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C .  3 

0.0 
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0.0 
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0.0 
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0.0 
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n.i 

0.5 

n  ^  p 

.50000000  0  0 
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0.3 

0.0 

.50000 GOO  0  C 

4 .7001 

0.1 

0.3 

0.0 

.50000000  0  0 

4. 7001 

0.1 

0.3 

o.: 

.  5  0000  nPf;  0 

4  .  7">n  1 

n.i 

0.3 

n  r> 

•  ■ 

.50000000  n  n 
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0.1 

0.3 

P.f 

.50000000  0  C 

4 . 7001 
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0.3 
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.50000000  0  C 

4.7001 
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0.0 

.50000000  0  0 

4  .  7  0  0 1 

0.1 

0.  5 

0.0 
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0.1 

0.3 

O.n 
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0.3 
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4 . 7001 
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0.0 
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0.0 
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4 . 7001 

0.1 

C .  3 

0.0 
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APPENDIX  C 

THE  MONTGOMERY  MODEL 
(MODIFIED) 


The  Montgomery  model  was  Implemented  as  a  FORTRAN  program.  It.  has 
run  on  both  an  IBM  370/168  and  an  IBM  3033-  On  the  IBM  3033  a  nl n<  ‘.y-mlnute 
simulation  requires  approximately  twenty-eight  CPU  seconds.  The  t'ol  lowing  table 
lists  the  major  variables  of  this  model ,  Following  the  listing  of  the  rode  I , 

Its  Input  data  and  a  typical  output,  there  are  Instructions  for  using  It. 

MAJOR  VARIABLES  OF  THE  MONTGOMERY  MODEL 
Symbol  Definition  Units 


A(N) 

Area  of  the  mid-pl&ne  between  compart¬ 
ment  N  and  N+l 

M2 

AL(I) 

Wet  suit  half -thickness  covering 
segment  I 

M 

AT 

Adipose  tissue  weight 

KG 

BFB(N) 

Basal  blood  flow  to  compartment  N 

l/hr 

BFWT 

Body  fat  weight 

GM 

RM 

Diver -s  basal  metabolic  rate 

KCAL 

M^-HH 

C(N) 

Thermal  capacitance  of  compartment  N 

KCAL/' 

CM(N) 

Radius  to  the  center  of  mass  of  compart¬ 
ment  N 

M 

DO(I) 

Outside  diameter  of  Segment  X 

M 

HAR 

Height  of  arms  segment 

M 

rfF 

Height  of  feet  segment 

M 

HHA 

Haight  of  hands  segment 

M 

HL 

Height  of  legs  segment 

M 

HSS(l) 

Thermal  conductance  between  skin  and 
wet  suit  of  segment  I 

kcal/hr-‘ 

HT 

Diver's  height 

CM 

HTF 

Height  of  trunk  segment 

M 

K(L) 

Thermal  conductivity  values  of  different 
materials 

KCAL 

LBWT 

Lean  body  weight 

GM 

NAT 


Non  adipose  tissue  weight 

C-l 


KG 


Symbol 

Definition 

Un  J  ts 

PBF 

Diver’s  percent  body  fat 

N.D. 

qb(n) 

Basal  metabolism  of  compartment  N 

kcal/hr 

QC 

Total  basal  metabolism  of  body  core 
excluding  head  and  trunk  core 

kcal/hr 

QM 

Total  basal  metabolism  of  muscle 

kcal/hr 

QS.F 

1 

Total  basal  metabolism  of  skin  plus  fat 

kcal/hr 

r(n) 

Radius  of  compartment  N 

M 

rmp(n) 

Radius  to  the  raid -plane  between  com¬ 
partment  K  and  N+l 

M 

s(i) 

Outside  surface  area  of  segraant  I 

M2 

SA 

Total  body  surface  area 

CM2 

SG 

Specific  gravity  of  diver 

NcD, 

swt(n) 

Body  compartment  (N)  weight 

Kb 

tc(n) 

Thermal  conductance  between  compartment 

N  and  N+l 

kcal/hr-°c 

THWS(l) 

Wet  suit  thickness  covering  segment  I 

M 

VAR 

Volume  of  arms  segment 

M3 

VF 

Volume  +  feet  segment 

M3 

VHA 

Volume  of  hands  segment 

M3 

VL 

Volume  of  legs  segment 

M3 

VTR 

Volume  of  trunk  segment 

M3 

wc 

Weight  of  body  core 

KCJ 

UM 

Total  weight  of  nLiscle 

KG 

WSCP 

Wet  suit  specific  heat 

KCAL 

KG-°C 

USD 

Wet  suit  density 

kd/m3 

WSF 

Total  weight  of  skin  plus  fat 

KG 

WT 

Diver’s  weight 

GM 

x(n) 

Heat  transfer  length  between  compart¬ 
ment  N  and  N+l 

M 

C-2 


Symbol 

De  f 1 n 1 1 1  on 

Un  1 1.^ 

C0ND3 

Total  body  heat  loss  to  ambient  water 

KCAL 

M2 -HR 

DEN 

Density  of  ambient  water 

ko/m-3 

DILAT 

Total  efferent  skin  vasodilation  command 

L 

HR 

UP 

Integration  step 

HR 

E(N) 

Total  evaporative  heat  loss  from  N 

KCAL 

UR 

EB(N) 

Basal  evaporative  heat  loss  from  N 

kcal/hr 

EH 

Total  respiratory  heat  loss 

KCAL 

HR 

errok(n) 

Output  from  thermoreceptors  in  compart¬ 
ment  N 

°C 

EV 

Total  evaporative  heat  loss 

KCAL 

M2  -HR 

F(N) 

Rate  of  change  of  temperature  in  N 

°c/hr 

H(I) 

Total  environmental  heat  transfer 
coefx'icient  for  Segment  I 

KCAL 

M2 -HR-°C 

HF(N) 

Rate  of  heat  flow  into  or  from  N 

KCAL 

HR 

HFLOW 

Total  rate  of  heat  flow  to  or  from 
the  body 

KCAL 

M2 -HR 

HO(I) 

Water-wet  suit  surface  heat  transfer 
coefficient  for  Segment  I 

KCAL 

M2 -HR-°C 

UP 

Total  metabolic  heat  production 

KCAL 

M2 -HR 

utbal 


Vot&l  body  heat  balance 


KCAL 
M2 -HR 


Un!  ts 


Symbol 

■  ■MWIW*  <•# 

H200U? 

INT 

rrirfE 

JTIH^S 

M 

MU 

P(I) 

PAIR 

PIN 

POUT 

PR 

Q(N) 

RATE(n) 

RE(I) 

RHEAT 

RHL 


S(I) 


0®  f  II  n  1 1  ton 

Total,  water  lou  rate  in  expiratory  gas 

Interval  between  outputs 
Elapsed  tin® 

Elapsed  time 

Metabolic  heat  product  lor  =*  HP 
(Print  out  label) 

E&maiaie  viscosity  of  water 

Water  vapor  pressure  table  from  0-50°C 

Vapor  pressure  of  environmental  air 

Saturated  water  vapor  pleasure  of 
inspired  air 

Saturated  water  vapor  pressure  of 
expired  air 

Pmndtl  Number  of  environmental  water 
Total  metabolic  heat  production  in  N 

Dynamic  sensitivity  of  thermo¬ 
receptors  in  N 

Reynolds'  Number  of  body  Segment  I 
Respiratory  heat  generation  rate 

Respiratory  heat  loss  rate 


Rate  of  body  heat  storage  =  HFLOW 
(Print  out  label ) 


Surface  area  of  Segment  I 


GM 

HR 

MIN. 

MIN. 

MIN. 

KGAL 

m 

KG 

HR-M 

mm  Jig 

mm  Hg 

mra  Hg 

turn  Hg 

N.D. 

KCAL 

HR 

HR 

N.D. 

CAT, 

MIN 

KCAL 

HR 

KCAL 

M2-hR 

M2 


C-*» 


Symbol 

Definition 

Uni  ts 

SA 

Total,  skin  surface  area 

M2 

SBF 

Total  skin  blood  flow 

L 

MIN 

SKINC(l) 

Fraction  of  vasoconstriction  command 
applicable  to  skin  of  Segment  I 

N.D. 

SKINR(I) 

J 

Fraction  of  all  skin  receptors  in 

Segment  I 

N.D. 

SKINS(l) 

Fraction  of  sweating  command  applic¬ 
able  to  skin  of  Segment  I 

N.D. 

SKJNV(l) 

Fraction  of  vasodilation  command 
applicable  to  skin  of  Segment  I 

N.D. 

STRIC 

Total  efferent  skin  vasoconstriction 
command 

N.D. 

SV 

Specific  volume  of  water 

m3/kg 

SWEAT 

Total  efferent  sweat  command 

KCAL 

HR 

T(N) 

Temperature  of  N 

°C 

TE 

Mean  weighted  body  temperature 

°C 

TD(N) 

Conductive  heat  transfer  between 

N  and  N+l 

KCAL 

HR 

T2 

Absolute  temperatuie  of  expiratory  gas 

°K 

TH 

Temperature  of  the  head  core,  repre¬ 
senting  the  hypothalamic  temperature  = 
T(l) 

(Print  out  label) 

°C 

TI 

Absolute  temperature  of  inspiratory  gas 

°K 

TIME 

Elapsed  time 

MIN. 

TM 

Temperature  of  muscle  compartment  in 
the  leg  =  T(J+7) 

(Print  out  label) 

°C 

C-5 


Symbol 

TO 


TR 

TS 

TSET(N) 

U(I) 

VE 

VEMXN 

WARM(N) 

WARMS 

WORKM(l) 

BULL 

CCHIL 

CCON 

CDIL 

CSW 

SCHIL 

SCON 

SDIL 

SSW 


Def I  nit  Ion 

Temperature  of  the  central  blood  compart¬ 
ment,  representing  esophageal  temperature 
T(6l) 

(Print  out  label) 

Trunk  core  temperature  representing 
rectal  temperature  ■  T(ll) 

(Print  out*. label) 

Mean  weighted  skin  temperature 

"Set  point"  or  reference  point  for 
receptors  in  compartment  N 

Overall  body  Segment  I  wet  suit- 
water  heat  transfer  coefficient 


Respiratory  gas  expiration  rate 


Units 


Fraction  of  total  work  done  by  muscles 
in  Segment  I 


Factor  determining  temperature 
sensitivity  of  sweat  gland  response 

Shivering  from  head  core 

Vasoconstriction  from  head  core 

Vasodilation  from  head  core 

Sweating  from  head  core 

Shivering  from  skin 

Vasoconstriction  from  skin 

Vasodilation  from  skin 

Sweating  from  skin 


KCAL 

-HR-C 

L 

HR 


Respiratory  gas  expiratory  minute  volume  L 

MIN 

Output  from  warn  receptors  in  N  °C 

Integrated  output  from  skin  warm  receptors  o_ 


N.D. 

iA 

KCAl/HR-°C 

lA 

l/kr~°c 

kcal/hr-0c 

kcal/hr-°c 

lA 

l/hr-°c 

kcal/hr-°c 


C-6 


Symbol 


Definition 


Units 


PCHIL 

PCON 

PDIL 

PSW 


Shivering  from  skin  and  head  core 
Vasoconstriction  from  skin  and  head  core 

Vasodilation  frcm  skin  and  head  core 

► 

Sweating  from  akin  and  head  core 


l/V 

I/V 

l/hb-^c2 

kcal/hr-°c 


C-7 


COMMON  /PASS/AL(6),D0f6),HS$(6),S(6),TC(60),QP(6Q),BH3(60) 

*,C<6  7>,CS(1),C8(1),RK(6) 

CALL  SIZE 
CALL  WE  f PAN 
STOP 
END 

SUBROUTINE  SIZE 

C  HIS  PROGRAM  IS  I  OR  A  FULL  WET  SUIT  HEAD  OUT 
REAL  *(  A  T  ,K 
D J*FNS  ION  K  <  3) 

DIMENSION  X  (  60 )  ,»MP(60)  , A(6'U 
D  I  ME  NS  I  ON  T  HVIS(o)  ,R(  86)  ,CM(  6  6)  ,N(6> 

DIMENS  I ON  S  W  T ( 6 1  ) 

CCMMOV  /  PASS/AL  (  6>,D0<6>  ,HSS<  6)  ,  S(6>  ,TC  (60)  ,QP(60),«BFB<6n> 

*,C  (6  7)  ,CS(1)/CB(1),R^o) 

DATA  WSCP/(J.?5/,WSO/24l  .C/,S6/0  .0/ 

DATA  THWS/0 . 706  3 12, 3.  (  0  6312,1  .006  31  2,0  .  r  0  6  3 1  2/ 0 „ 006  31  2  ,  i  .  006.31  2/ 
DATA  < /0. 31 93,0.1419,0. 046/ 

DATA  RK  /  0.0.3  71  1 ,0. 0371  1 ,0.0371  1,0.  *'.371  1,0.03  71  1,0.0371  »  / 

R  F  A  D  (  1  , 3  )  H  T  ,  KT 

3  FORMAT  (  2  F 1 U . b) 

R  E  A  T  ( 1  ,  7 ) PM 

7  FORMAT ( FI  0. 5) 

I  F  (SC.  GT  .  0. 0  )GO  TO  5 

SC,  =  0.8*  (  (HT**0.  242)7  (WT**0.  1  ))  ♦  0.1  6? 

8  f-BF=(5  .S48/Su)-8.‘j44 

BFWT=PB F  *WT 

LBWT=V  T -8FWT 

SA  =  71.  8  4*(<  WT/1  OO'')  **0. 425)  *  <HT**H  „  7?5  ) 

S (1  )  =0 . 0  7*S  A  *C . OCl  1 
S  (2)  =  J  .  3602 *  SA *  0  . on0 1 
S ( 3 )  =3  .  1341*  ram:. 0031 
S  C  4  )  =0  .  1 5*  S  A  *0  .  C'OOI 
S<5)=0. 3174* SA*  0 . DUO  1 
S  (6>=0  .  0686*  SA*C.  Dull 
A T  =  B  f  W T  *0  .OU  1 
NA  T  =L  3  W  T  *  U.  0 1)1 
DO  2  I =1,4 

2  SWT (I)  =  (J.J1 64/ (0.85*4. 0  )  )  *  N  AT* (0.024/ (U. 65*  4.0) )  *  \  AT 

DC  4  J =5,8 

4  SWT(I) =(u. 005/(0. 85*4. 0))*N ’  T 
SWT  <  V)  =(0,005/0. 1S)*AT 

SWT  (10)  =  <0.0036/  1.85)  *N  AT 
DO  6  1=11,14 

6  SWT  (I)  =(0.038/<i  .8  8*4.0  >>*N  AT*  (0.1  59/ I  0.8  5*  4  .')>*NAT 

DO  8  1=18,18 

P  SWT(  I)  =  (0.24  9/  (  ).8S*4.0))*NAT 

SWT  ( 19 )  =  <U. 095/U. 1 8>*  AT 
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16 

180 
2  OP 


22 

24 


2  6 


2  8 


300 
3  20 


340 

36 


38 
4  0 


4  2 
44 


SWT (23) =(0.0181 /0.85)*NAT 
00  10  I =21,24 

SWT(I)  = (0,0202/ ( 3.85*4 .Q))*NIAT*  (0.01/(0.85*4  .G ) ) *N A  T 
00  12  1  =25,28 

SWT(I)  =  (0.0453/ (0.85*4. C))*0 AT 
SWT (29 ) =(0.013/0. 1 5  )  *  A  T 
SW T < 33 >=( 0.0065/0. 85 )*N A1 
DO  14  1=31,34 

SUT ( 1 ) ”  (0.0031/ (0.85*4. 0))*NAT+( 0.0 004/(0.85 *4.0))* NAT 
DO  16  I =35, 3d 

SWT  (I)  =<0.001/  (0.85*4  .  0))*N  AT 
SWT<39 ) =(0.002/0. 15)*AT 
SWT (43 ) =(0. 0025/0. 85)* NAT 
DO  180  1=41,44 

SWT  (I)  =  (0.06  73  7  (0.85*4.0))*  WAT*  <0.025  8/ <0.85  *4.0))*  NAT 
DO  203  1=45,48 

SWTC  I)  =(0.1  3 68/ (0.85 *4.0))* VAT 
SWT (49 ) =<0.032/0.1 5)* AT 
SWT<  50  )  =  (0. 0161/0. 85)*NAT 
DO  22  1=51,54 

SWT ( I)  =  (0.00 5/ ( 0.85*4 .0  )  >*N A  T+ (0.0008/ (0 .85*  4.0) ) * V  AT 
DO  24  1  =55,  58 

SWTU)  =(0.001  /  (0.8 5*4.0  )>*N  AT 
SWT(59)*<0.0C3/0.15)*AT 
SWT (63 ) =(0.0032 /O .85) *N AT 
SW T  (  61  )  =2.50 
DO  26  1=1,4 

C  <  I)  *<  0 . 5*  <  C  .01  64  /  0 . 85  )  *NAT  ♦  P.  0  *  <  0 . 024  /C  .  8  5  )  *N  AT  >  /  4  .0 
DO  28  1=5,8 
C  (I  >-~3  -9*SWT  U) 

C(9)=r.l  .6*  SWT  (9) 
r.  (1C)  -  n  .  9*  s  w  t  ( 1  p  ) 

DO  3  JO  1  =  11,14 

C  (I  >*<  0 .5*  (0  .038/0.85)*  NAT*  ?  .9*  CO.  159/0.8  5)*  NAT)  /  4 . 0 
DO  323  1=15,18 
C  < I ) =0 . 9*SWT  (I ) 

C  1 19)  =  0 .6*$W  T(1 9) 

C(20)=U.5*SWT(2P) 

DO  343  1=21 ,24 

C  <  1)  =  <  0 . 5*  (  0  .02  C*2  /  0 . 8 5 )  * NA T  ♦  0.  9*< P  .  01  /  O.  8 5 )  *  NAT  >  /  4 . 0 
DO  36  I  =  25,28 
C U )=0 -9*SWT  (!) 

C <29>=0  .6*SWT(29> 

C  (30)=0  .9*SWT(3u) 

DO  38  1=31,34 

C  ( I  )  =  <  0 . 5  *  (  3 . 3 !  C  3 1  /  0 . 8  5  )  *  8 A  T  ♦  0 .  9  *  (  0 . 0  0  0  4  /  0 . 8  5  )  *  N  A  T )  /  4 . 0 
DO  4  1  =35,  3  8 

C (1) =3 . 9*SWT  (1 > 

C  (  3  V  >  =  0  .6*SWT<39) 

C  <  4  0  )  =  0 . 9*  S  W  T  ( 4  0  > 

DO  42  I =41 .46 

C  (1  )  *<  0 .5*<0  .06  73/0.  85  )  *.«AT*  0.9*  <0 .0258/0.85  )*N  AT)  /  4.0 
DO  44  I =45, 48 
C (I>=3 .9*SWT (1 ) 
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C  (  49  )  =  0 . 6*  S  W  T<  4  9) 

C ( SO)  =  C .9*SW  T(5'D 
DO  4  6  1=51,54 

4  6  C  (  I )  =(  0  „  5*  (0  .00  5  /O  „  8  5  )  *  \  AT*  r,  .9  *  <0  .  v'008  /  85) *NAV> /4  .0 

oo  48  1=55, sa 

48  C ( 1) =3. 9*SWT (I ) 

C(59>=0.6*SWT<59) 

C(60)=C  .9aSWT<60) 

C  (61  >=•'  .  9  *  S  W  T(6  1  ) 

00  53  1=11*14 

so  c  u  )  =c  ( i  )-c  <  6i  >  /n.  o 

CS  ( 1  )  =  C  (10)  ♦C(ZO)  +  C  (  3-7)  ;  0(40)4  c  (50  )  +  C(  61) 

CB  <1 )=0 .0 
00  5  2  1  =1  /61 
CB (1 )  =  C  8  <  1 ) *C<I  ) 

52  CONTINUE 

«  S  f  =  J.  3  *  (0. 0588  2  *  N  A  T  +  1  „  0  *  A  T  > 

<1F'  =  (0.  1  8* (6tt*SA  /  1  u  JO,;  }  )  OS  F 
«C  =  D  .1  0  *  (BM*  SA  /  1  0(  03  ) 

WS F  =  n. 0  5  8  8  2  *  N  A  T  ♦ 1 . 0*  A i 
Wv  =  0 .5  0  5  9*N!  AT 
WC=0.1  7fc5*NAT+0.?588*9 A1 
00  54  1=1,4 

5  4  OB  < I )  =  SWT (I  ) *QC /  WC  +  0  .04*  (BM  *  SA /10000) 

DO  56  1=5,  rt 
56  GO  ( I  )  =  SWT(.t  ) 

Gf  (  9)=  S  WT<9)  * G S  F  / n S  F 
GP<10)=SUT(1r')*tt5F/J5F 
DO  58  1=11,14 

58  OB  <  I)  =  S  WT(  I  )  *0C  /  WC*3. 1  4*<Bfl*  SA  /  10930) 

DO  67  1=15,1? 

60  CP <1 )= SWT<I ) •QM/WM 

GB(  19)  =  S  WT (  1  9)-MJSF/WSF 
Q  B  (  2  ))  =  S W T  (  2  •")  *  Q  S  F  /  JSF 
DO  62  1=21,24 
6?  C.HU  )=  SWT  (1  )  *GC  /WC 

DO  64  1=25,28 
64  00  (  I  )  =  S  WT  (I  )  *<m/WM 

DO  66  1  =29,  3': 

66  wB ( I )= S WT ( 1 ) *GS F/WSF 

DO  68  I =31, 34 
68  QFi(  1)=SWT(X)  *«C/WC 

DO  70  I=35,3/i 
70  QB  (  I  ,=  S  WT(  I  ) 

DO  7?  1=30,4  ’ 

72  QB < I )= SWT  (I  ) *QS f / WS f 
DO  74  1=41,44 
74  OH  II )=SWT  II  )  *0C  /VC 
DO  76  1  =45,4  " 

76  0  B (  I ) =  SWT (I )  *  Q I / W  M 

DO  78  1 =49, 5  j 

7  8  9B(I)=SWT(I)*US  F/WSF 

DO  8  7  1=51,54 

8  0  QB ( I )= SWT <1 ) *oc  /WC 
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DO  82  1=55/58 

82  QB  ( I  )  =  SWT  ( J  )  *OM/WM 
DO  84  1=59,6 

84  QB(l)=SWTll)  *USr/k'SF 

DO  86  1=1,60 
8  6  BFB(I)  =1  .?*«0<I  ) 

DO  88  1=1,4 

83  BFB< I) =1 1  .2  5 
DO  91  1=11,14 

RQ  R  FB ( I )  =  62.5 

BFB(13  )  =5.34  *SVT»(1  2) 

B  F  P  (  23  )=1  .56*SWT(2Q) 

BFB<33)=1.35*SWT (30) 

BFB (40) =5.58 *SWT (40) 

BFB ( 50 ) =0.84 *SW  T (50) 

P  FB  (  60  )  =2.34  *SW  T  (60) 

R ( 1 )=< ( (0.003*SWT(1) )/1 2.56 ) **0.333) 

DO  92  1=2,10 

92  R  (  I  )  =  (  (  R  <  1  -1  )**3.0)+(<0  .003  *  SW  T  (I)  )  / 1  2  .56))*  *0.333 

R<  61  )  =  R ( 1 0)  +  1HW  S (  1 ) 

V  T  R  =  0 .  0 

DO  9  4  I  =  1 1 , 2  C 
94  VTR=VTR+SWT(  I)*O.Ot!l 

8TR  =  (S  (  2  )  *  *  2  .(')  /  (12.56*VTR) 

C  **  *  *  *C  H  A  NGi  A  TO  AA  .  A  IS  PREVIOUSLY  DIMENSIONED.***** 
AA=n.301 

R  <?(■)=  (  2*VTR  )/S  (2) 

DO  96  1=1,9 
J  =  ?1-l 

R(J-1)=((R(J)**2)-((AA*SwT(J)>/(3.14*HTP)))**0.5 
96  CONTINUE 

R ( 6  2 )  =  R (20) +THU$(?) 

V  A  R  =  d  .  0 

DO  93  1=21,3. 

98  VAR=VA R +SWT < 1 )* 0 .001 

HAK=(S<  3>**2  .  ')  /  <  1  2  .  ^  6*  V  'F  ) 

R  <3'0)=  <2  .0*  VAR  )  /S<3) 

do  iro  i-i-9 

J  =  3 1  -I 

R  (J-1)  =  <  <R(  J )•*  ?>- (( A A*  SWT ( J  ))  /  (3. 1 4*M AP)  )  )*  *P.  5 
1J0  CONTINUE 

R  ( 6  3  )  -  R (30) + THWS (3) 

VHA=  0.  0 

DO  1  02  I=?1,  40 

102  VHA  =  VHA*SWT  (  I  )*•''.  .501 

HMA= (S  <  4 )**?  .0) / (1 2. 56*  VPA) 

R  U0>=  <  2  .0*  VpA  )  /S<4> 

D  0  1  !4  1=1/9 

J  =4  1  -I 

«(j-1)  =  (  (R  C  J  )**?>-((  ^  A  *  E  wT  (J  )  )  /  t  3 .  1  4  *HHA ) ) ) *  *0.5 
1  34  C  C  NT  l  M  U  F 

R  (  64  )  =  R  (  40)  +  IHW  S  (  4  ) 

VL  =  n  .3 

DO  1-6  1  =41  ,  53 
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106  VL  =  VL*  SWT  Cl  )  *0.  001 

HL=(S( 5  )**2 . L ) / (12 .56* VL) 

R  <  50)  =  (  ?.0*VL)/$(5) 

00  1  >8  1=1,9 
J  =  5 1  -I 

R(J-1)  =  ((R(J)**2)-((AA*SWT(J))/<3.14*HL)))**0.5 

108  CONTIS  UF 

R(  65)  =  R  (  50)  +  THWS(  5) 

VF=0 .3 

DO  113  1=51, 6U  , 

110  VF=VF* SWT(I) *0. 001 

HF  =  (  S(  6)**2. 0)/  <12.5  6*V  F) 

R(60)=  <  2.0*VF)/ S(6) 

DO  112  1=1,9 

J  =61-1 

R  (J  -1)  =  (  (R(  J  )*  *  2 )  -  (<  A  A  *  S*T  (J  ))  /  (3  . 1 4* H  F  )  )  )**  0  5 
112  CONTINUE 

R (66 ) =  R  ( 60) ♦ THW  S ( 6 ) 

CM(1  )=0.7937*R(  1) 

oo  m  1=2, iu 

114  CM(1)s0.793  7*((R(I-1)**.y)+(R(1)**3))**0.333 

C  M  (  6 1 )  =  '_• .  79  3  7*  (  <  R  <  1  0)  *  *  3  >♦  (  R  <6  1  )*  *  ■*  )  )  *  *0. 3  33 

Cv(  H)  =0.7071*R  (1  1  ) 

DO  116  1=12,20 

116  CM(I)=0.7071*((R(1-1)**?H(P(I)**2>)**3.5 

CM(<>2)  =C.70 71*  (  (R(  20)**  2)+<R  (62  )**2>>**0  5 
CM( 21) =0.7071*K  <  21 ) 

00  118  T  =22 , 3n 

118  CM(l)=n.7Q71*((R(l-1)**2)*<R(I)**2))**).5 

CM(63)=0.7G71*((R(30>**2)+<R(63)**2))**r).5 

CM  (31)  =  u.70  71*P  <  31  ) 

DO  173  1=32,40 

120  C'1(l'  =  0.7r)71*((P(l-1)**2)+(R(i)**?))**o.4 

Ci i  (64)  =0.707  1  *  (  ( r:(  4  ))**?)♦(  *  (64  )**?))*  *o  5 
C*(41) =(  .  70  7  1  *  R (41 ) 

DO  1  22  1  =42,50 

122  C  M  (I  )  =  0.?0  71*(  (R(  i-l)**  2)M  Mi  )**2))**J. 5 

CM  '6  5)  =0.707  1*<  (P(5))**2)+(R(65)**2))**^  5 
CM ( 5 1 ) =0.7071*R (51 ) 

DO  1  24  1=52,60 

124  CM<I)=0.7O71M(?(1-1)**?)*(R<i  )**?))**  r.  5 

CM  ( 6  6)  =  0 .7C  7  1*  (  (*<(*> 3)**2)+('M66)**2))**r’  5 
DO  126  1=1,60 
126  X  ( 1  >  =C  K  (  J*  1  )  -CM  (  1  ) 

DO  123  1=1,6 

128  X(10M  )  =  CM  <  ]  ♦  6  3  )  -  c  M  (  1  'u  l  ) 

DC  130  J  =1,6  1 

130  RMP ( I )  =  CM ( I  )  +X ( I  )  /  7.  f 

DO  132  1=1,11. 

132  A  ( I  )  =1  2 . 56*  (  RMP  ( I  )  *  *2) 

DO  134  1=11,23 

134  A( I ) =6 . 28*RMr ( 1 ) *H TR 
DO  136  1=21,3' 

1  36_  A  (  I )  =5 . 2S*RMP(  I  )  *  H  A« 
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DO  138  1=31/4(1 

13K  A  (  I  )  =6  .  28*RWP<  I  )  *HHA 
00  140  1=41/50 

140  A(I)=6  .28*RMr><I)*HL 

DO  142  1=51/60 
142  A(I)*6.28*RHP(I)*HF 

DO  144  1  =1  /o  : 

144  TC<I)= «  K  C 1 >*A(I >  > / X  < I) 

RMPA=CM (H)+ (  (R ( S) -C 4(8)  )/2. " > 

AA=12.56*<RMPA**?> 

XA  =  R (8 ) -CM( 8  )  * 

Y  A  =  (  K  ( 1  )  *  A  A )  /XA 

RMPB=R  (8)+((CM(9)-R(8))/2.0) 

AB  =  1  2  .  56*  (RMHF-*  *2) 

XF=CM<  9 ) -R ( 8  ) 

YB  =  ( K( 2 ) *AB >  / XB 
iC  (  8  )=  ( YA*YB )  /  (  YA  ♦  Yf)  1 
00  145  1=1/6 
146  T  C  (  1  0*  1  )  =0 . 0 
DO  148  1=1/6 

I  48  T  C  (  1  >  I  -  1 )  =  (  X  (  ?  >  *  A  (  1  I  - 1 )  )  /  x  (  1 '  *  T  -  1  ) 

H  <1  >=0 .0 
H (?)-H  TR 
H  (  3 ) =4  A  R 
H  (4 ) =1 H A 
H  (  5  )  =H  I. 
h  (6)=H F 
DO  150  1=2/6 

RP»PA1=  CM1T*  I-?>  ♦  (  <R  <10*1-2  )  -CM  (10*  W  >  )/?  .0  ) 
A  A  1  =  6.  2S«RMPA1*M  (  I  ) 

X A  1  =  R (  IP *I-2)-Cw(1  0*1-2) 

YA 1 =(< (1>*AA1)/>A1 

RKPR  1=  R  ( 10*  1  -2)  ♦  <  (C'i  <1u*I-1  )  -»  <  10*  1  -2)  >  / 2 .0) 
AH  1  =  6 . 2  8*  RM  P  Hi  *  H  (1  ) 

VB  1  =C8  (  1  u*  I  -  1)  -  8  (  10*1-2) 

VBl =(K  ( 2)*AB 1) / XP1 
T  f  (  1 O  1  -  2  >  =  (  YA  1  *  Y  » 1  )  /  <  Y  A  1*  Y  P  1 ) 

150  CO"ThUf 

SMPR=CM  (ID*  <(R  <  10)  -  CM  (  10) )  /  ?) 

SWF  A  =  1  2 . 56*  (  SMPR*  »  2) 

SX=«<1  n>-Cl"(  in) 

SV=(X(2)*SMPA)/SX 

WSHPR=  R  ( 10)  ♦  <  <  C  (  1 1  )  -R  (  1 '!)  )  /  2) 

WSMPA=1  2.56  *  (WS 8  PR  **2  ) 

WSX=C'1  <61)-R(10) 

W  5  Y  = <K  ( 3>*WSMPA ) /WSX 
8SSU)  =  (SY*WSY)/(SY*WSY) 

OO  152  1=2/6 

S«PPIsC8(K»«l)H(Mn»l  )-CX(1j«I))/2) 

SP»PA  1=  6 .28*  SBPR  !  *8  (  I) 

SX1=R(  1  0*1  )  *  CP  (  K  *  1  ) 

SY1  =  <<  (  2  )  *S  iPAl  )  /SX 

WSMPR1  =  R(1U*I)*(tCH(I  +  6f  )-«  <  1'»«  1  )  )  /  2 . 0  ) 

WSPIPA1  =6. 2fc  *  WSKPR1  *H  (  I  ) 


WSX  1=CM(  1  +  60  >-R  (10*1) 

W$Yl=!K!3)*WSMPA1)/WSXl 
HSS(I) =(SY1 *  WS Y  1 ) /!SY1*WSY1 > 

152  CONTINUE 

DC  154  1=1/6 

154  AL  (I  )=  R  (I+6C  )-CM  (I  +6'  ) 

DO  156  1=1/6 
156  D0(1 )  = 2 .0*R ( 1460) 

S(1>=1 2 . 56* !R!61)**2) 

DO  158  1=2/6 

158  S  <  I  )=6.28*R!l*6*r)*H!I  ) 

C(62  )  =  WSCP*WSD*1  .333*3.  14*!  !R!61)**3)-(R  (  10)  **.*)) 

DO  IfcD  1=2/6 

160  C  ( 1*61  )  =WSCP*WSD*  (  (  R  (  I  ♦  60)  **2>~<RC10*1  )  *  *  2  )  )  *H1I)*3.14 

s  a  a  =  s  a  /innon  .0 

VH  F  A  D=  0 . 0 
DO  162  1=1,10 

162  VHEAO=VHEAD+SWT!l > *0.001 
WRITE! 3,361) S6,PBF,SAA 

361  FORMAT  (IX/'BODY  S  .  &  .  =  •  /  F  7 . 5  /  ?X  /  *?.  BODY  IAT  =  '/F<’.4/ 

12X/ ' SJ RF ACE  AREA= * , F7.4  ,  *SQ  M '  ) 

WRI TE! 3,37) AT/NAT 

37  FORMAT (/,  IX,  ‘AD 1P0SE  TISSUE  WT .  =  * / f  7 . 4 / 1 X / ' < G ' / 2 X / 

& ' VON -A  D I POS  E  TISSUE  W  T  .  =  ' ,  F  7. 4/1 X/* KG*) 

WRI  TF(3/30)R  (UD/VHEAD 

3  Ci  FOP  MAT  (  /  ,  1 X  /  "HEAD  R4  D  I U  S  =  *  ,  f  7.  4 ,1  X  ,  *  M  ’  ,  T  3  0  ,  *  HE  A  D  VOL.=*, 

8F7.4/1  X,  *M  CUBE  *) 

WRITE! 3/31) H IF/ VTR 

31  FORMAT  (  IX/' TRUNK  L  t  N  C- T H  =  •/ F  7  .4  /  1 X  ,  •  M  •  ,  1  30  /  •  T  PUNK  VOL.=*, 

5  F  7. 4,1  X /  'M  C UPF  *  ) 

WRITE! 3 / 32 ) H  AR , V  A  R 

32  FORMAT  !  IX/’ ARM  L  £  N  GT  M=  *  ,  f  7 . 4  ,1  X  /  *M  *  /  T  3  P  ,  *  AR  M  VOL.=  *, 

!v  F7.4,  1  X,  *M  CUBE  ' ) 

WR I TE!  3 / 33) H  HA , VH  A 

33  FORMAT  (  1  X,  '  H  AN  n  LfNjTH=',F7.  4,  IX,  *  M  * ,T  30 /  'HAND  VOL  .  =  *  , 

R  F7.4/1X/'M  CUBE*) 

WRITE! 3/34)HL/VL 

34  FORMAT  !  1  X, ' l  EG  L E N G T H = * , F 7 . 4 , 1 X / * 1  * / T 3  ' / * L E 3  VOL.=*, 

6  F7.4,  1X,’M  CUBE*) 

WRITE! 3/35)hf/VF 

35  FORMAT  !  IX /'  F  OOT  L ENS T H=  1  ,F 7 . 4, 1 X , • M • , T 30 / * F 3 OT  VOL-=*/ 

%  F  7.4,  1 X,  *M  CUHE  ' ) 

WRITE! 3/181) 

181  FORMAT!//*  FOLLOWING  AR  I  G  <  I  )  VALUES  IN  SU  M  * ) 

WRITE! 3/51) ( S( I )/I=1/U) 

51  FO»MAT  !  1!6F8  .4,/)  ) 

WRI Tf  !  3, 1 1 3  ) 

183  FORMAT!//*  FOILO-ING  ARE  AL  (  I)  VALUFS  IN  M*) 

WHITE! 3,19)  (AL(I)/I  =  i,6) 

19  FORMAT  (  1  (6f  8  .4  /  /  )  ) 

WRITE!  3  /  165  ) 

185  FORMAT!//’  FOLLOWING  ARE  DO  (  I )  VALUES  IN  1*) 

WRIT  E<  3/20) ! DO ( 1 ), 1  =  1,6  ) 

20  f ORMAT  ( 1 <6F8  .4, / )  ) 
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WRITE! 3,187) 

187  FORMAT!/,*  FOLLOWING  ARE  HS S  !  J  )  VALUES  IN  KCAL/MR/DFC  C’) 
WRITE(3,18)<HSS(I),I=1,6) 

18  FORM  AT  <1!6F8.4,/)) 

WRITE! 3 ,189) 

1S9  FORMAT  (  /  ,  '  FOLLOWING  ART  TC(I)  VALUES  JN  KCAL/hR  DEG  C') 
WRI Tfc( 3,17) ( TC( I) ,1-1,60) 

17  FORMAT  <  6  <  1  G  F  7.?,/)  ) 

WRI  TE ( 3 , 191 ) 

191  FORMAT  (/,*  FOLLOWING  ARE  A(l)  VALUES  IN  SQ  M  ' ) 

WRI  Tt(  3,163)  !A<  I  >,1  =  1,6  0) 

163  FORMAT  (  6  ( 1  OF  7. 2  , /  >  > 

WRITE! 3 ,193) 

193  FORMAT!/, ’  FOLLOWING  ARE  RMP(I)  VALIES  IN  M  *  ) 

WRITE! 3,164)  (RMPCT), 1=1  ,60) 

164  FORM  AT  (  7(1  OF  7.4,/  )  ) 

WRITE!  3 ,195) 

195  FORMAT!/, *  FOLLOWING  ARE  X !  I  )  VALUES  IN  M') 

WRITE! 3,165) (X( I), 1=1, 60) 

165  FORMAT ( 7(10? 7.4,/) ) 

WRITE(3,197) 

197  FORMAT!/,*  FOLLOWING  ARE  THE  VALUES  OF  CSd)  AND  C3(1>*) 
WRITE! 3 ,47)C  S!1 ) ,CB (1 ) 

47  FORMAT  (1X,'CS(1  )  =  *  ,  F  7 . 3  ,  T30  ,  *C  8  ( 1  )  =  *,  F  7.  3  ) 

WRIT  E( 3,199) 

199  FORMAT!/, *  FOLLOWING  ARE  C(I)  VALUES  IN  KCAL/DEG-C’) 
WRITE! 3,167) (C( 11,1=1,66) 

167  FORMAT ! 7!1GF7. 3,/) > 

WRITE! 3,701 ) 

201  FORMAT!/, *  FOLLOWING  ARE  Q8!I>  VALUES  IN  KCAL/HR*) 

WR  I  TEC  3  ,169)  (QB  ( I  )  ,1  =  1,60) 

169  FORMAT ( 7(1  OF  7.3 ,/) ) 

WRITE! 3,203) 

203  FORMAT!/, *  FOLLOWING  ARE  OFB(I)  VALUES  IN  L  /  HR  *  ) 

WRITE (3,1 71) <Bf P(I),I=1 ,60) 

171  FORMAT  <7(10F7. 3,/)  ) 

WRITE!  3 ,205  ) 

205  FORMAT!/, *  FOLLOWING  ARE  R(I)  VALUES  IN  M') 

WRI TE! 3, 173) (R! I >,1=1,66) 

173  FORMAT  (7(10F7. 3, /)> 

WRITE! 3  ,207) 

207  FORMAT!/, '  FOLLOWING  ARE  CM(I)  VALUES  IN  M ’ ) 

WRITE(3,175) <CM!I),I=1,66> 

175  FORMAT  (7<10F7. 4,/)) 

WRITE! 3 ,209) 

209  FORMAT!/, *  FOLLOWING  ARE  SWT  (I)  VALUES  TN  KG*) 

WRITE!  3,  177) (SWT! I), 1=1 ,61) 

177  FORMAT  <  7!10F7. 3,/)  ) 

RETURN 

END 
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SUBROUTINE  WFTMAN 
REAL  <Ut 

DIMENSION  T  <67),TSET  <61  ),RATC(61),F  (67) 

D  i  'IE  NS  I  ON  WARM  <  61  )  ,COLP  (6l  )  ,HF  (67)  ,  FRROR  <  61  > 

I)  I  MENS  I  ON  E  P  (60), 0(60), E  <6U  )  ,B  F  (oO) 

DIMENSION  TD 1 60), DC  I  60) 

OT IF NS  I  ON  PI  11 ) *  SKINS  <6  ),SX 1 NV (6) , SKIN  C (6) 

DIMENSION  R  F  (6)  ,HO  '6 ) , U (6) 

DIME  NS  I  ON  SK  I  NR (6) 

DIME  NS  I  ON  W0KKM<6),CH!LM<6),H<6>,PSKIN<6),E'IAX<6> 

D  I  MF  NS  I  ON  A  T IHE  (9)  ,  ATA1  R<9)  , AV<  9), ARM! 9) , AW3  RK  <9) 

DIMENSION  A ( ?08 ) , 0 (28 R) 

COMMON  /  PASS  /  Al  <  6)  ,  DO  <  6)  ,HSS  <  6)  ,  S<6>  ,TC  (60)  ,B  0(60  )  ,K  FP  (  60) 

*,C(6  7),CS<1)*CP(1),RK<6) 

DATA  A  T  J  n:/0  .,3'J.,oO.,9  0.,1  ?().,  150  .  ,1SC.,  213  .,250./ 

DATA  ATAIR/1. 7,1. 7,1.7, 1.7, 1.7, 1.7, 1.7, 1.7,1  .7/ 

DATA  A  T  A  IK  /  1  1  .8 -1  1  .8, 1  1  .8,1 1  .0 , 1 1 . 8,1 1  .8 , 1 1 . P,1 1  .8/ 

DATA  A  V  71  ??.,12?.,1?2.,122.,12?.,12?.,1??.,12?.,12?./ 

DATA  A  Rh/P. 3 ,0. 3*0.3 ,0.3,0. 3*0 .3/0 .3*0 .3,0.3  / 

DATA  AWORK/(.,0.,0.,U.,0.,0.,0.,Q.,0./ 

DATA  P/4. 5 7 5, 6. 543,9.209, 12. 788, 17. 5 35, 23. 75 6, 31. 82 4, 4?. 175, 
55.  124,71.88,92.51/ 

DATA  TS  F T/16  .96,36 .96,36 .96,36 .96,35.07,35 .0  7,35.07 ,35.  7 , 
r.  34. HI  ,34. 58  , 36  .89,36  .89,36.69,36.89,36.28,36.28,36  .28,36.28, 

K  34. S3  ,  3  3.6  2,35 .53,35.53,35 .S3 ,35. 53,34. 1 2,34.12,34 .12,34 .1?, 
k  33.  59, 3 1.25,35  .41 ,35.4  1  ,35  .41  ,35. 4U35.  3  b,  3  5.38,35  .38,35.38, 
c  35.  3D  ,35.  2  2  ^3  5 .81 ,35 .8  1,35 .81  ,35.81,35.3,35  .3/3  5.3,3  5  .  3,3b.  31 
3  4  .1]  ,35.14,35  .14,35.14,35  .  14,35.  03,35.03,3  5. 03,33  .03,35  .11* 
i  3  5.04 , 36.71  / 

DATA  F  B  /  60*0  ./,  RATE/ 51*  0./ 

DATA  SK INR/O .0695,0. 4«3  5,0. 0686,0. 1843,0. 1  SO  5,0 .033  4/ 

Data  SKINS  / 0.P8  1,0. 481, 0.154  , 0.031, 0.218,0.3  35/ 

DATA  SKINV/0. 13  2, 0.322, 0.09  5, 0.121, 0.233, 0.1  0/ 

DATA  SKI  N  C/0. 05, 0. 15, 0.  0  5,  0. 35, C. 05,0. 35/ 

DATA  Al  ORKM/  c  .0,  O.G75,  0.  02,0 . 00?  5,0.  15,  0.  0^2  5  / 

DATA  CH ILM/O .00 5,0.2125 ,0.01 25 , 0.0 ,0. 0 1 75 ,0 . 0/ 

DATA  C  S  U/32C  .0/  ,f>SW/?V.  0/,PSU'/0 .0/ 

DATA  C  D IL/1 1 7.0/, SOIL / 7  .5/ , P  D I L  /  0  0/ 

DATA  CCHI  1/0.0  /  ,SCHlL/0.f)/,PCHIL/21  .0/ 

DATA  CCON/5.  V  ,  S  C  ON  /  5 . 0  /  ,PC  0  N  /  Q  .  '>/ 

DATA  MULL/  10.0/ 

DATA  T  /  3  /  .7 , 37 . 7 , 3  7 . 7 , 3  7 .7,  3  7  .  ,  37 .  ,  ^7  .  ,  37  .  ,  3  6  .  ,  3  5 . 3  , 
k  37. 7,  37. 7, 37. 7, 37. 7, 36  ., 3  A., 3  6. ,36. ,36. ,3*'.  3, 

£  -5. ,35. ,15. ,35. ,35. ,35. ,35. ,34. ,3  4  . ,34., 

8  7 5. ,35. ,35. ,3 5. ,35. ,  3 5., 35  . ,35. ,35. ,34. 5, 

K  3  5.81  ,  35. 81  ,35.81  ,35.0  1,35.3,  35.3,35.3,35.3  ,35.31,  34.  1, 

?,  3 5. ,35  .,35. ,35  .,35. ,35  .,35  .,35. ,35  .,32.8,35  .71, 

?,  34 .55 , 3 3. ,33. ,3 3. ,34  .,  3?./ 

DATA  3M/39.4/ 

KK  =  ■' 

J  J  =  4  n 

DO  99  1=1,6 

RATE  (1  0*  I  )  =  0  .03 

READ  (1  ,  5  )  PR  ,  MO,  SV,COND 
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5  FORMAT  (4F1.1.5) 

R  F  A  D  ( 1  *A)TAM*PRA*MU4*SVA*CONDA 

6  FORMAT  (  5f  IT.  S) 

S  A  =  0 . 3 

r>o  1  1  r>  K  =  1  .•  6 

110  5  A  ■-  S  A  ♦  S  ( K ) 

T  T  M  F  -'O  .  0 
IMT  =  S 
I  T  I  M  E  -  > 

J  T  I M  E  =  0 

00  10?  N  =  1  *  6  1  ' 

r  <  M  )  =0 . 0 
1  j2  CON  T  I  V  U  E 
1=1 

199  I  f  <  I  TI  ME  .LT  .  ATIME  (L*1  )  >  00  TO  189 
L=L  +  1 
GO  TO  199 

189  TAIR=A TAIR(l > 

V  =  A  V  (L  ) 

RH  =  A  RH  ( L  ) 

WORK -k WORK (L ) 

I  F  (WORK -OiM*SA)  )  1  U4 , 1  0  4  A  0  5 

104  WORK  1=  Cl.  0 
GO  TO  1  CO 

105  WORK  1=  (  WORK  -  CBM  *  S',  )  )*0.  '8 

1 06  CONTINUE 

0 FN A  =  1  .  P/SVA 

RC (1 )= ( (PENA-V*  00(1 ) ) /MU A) 

MO  (1  )=  (  CONDA/DO  (  1  )  )*  (  0. 9  7  +  0 . 68*  (RE  ( 1 )  *  *h.  5) )  MPR  A**  0.31 
DC  ?  I -2*6 
DFN=1. O/SV 

RE  ( I  )  =  (  ( DFN*V*DO(I) )/MU) 

HO  (I  )=  (  COND  /  DO  ( I  )  )  *  (O.ftOMRF  (I  )**0 .50>  *  (PR**  0.31  > ) 

?  CONTINUE 

HO  ( 1  >  =  (COND  /D0(  1  )  )*  (0 .9  7*0. 6  8*  (PE  ( 1  >**  0.  5)  >*  <PR**0.  3) 

U( 1 )=1 . 0/( ( 1 *0/h0(1) )♦ ( AL(1 )  /RK  (1 > 

U(?)=1  .  CJ/(<  1  .0/  F.0(?)  )  +  <  AL(?)  /RK  (?) 

U(3>=1  .()/(<  1  .0/P0(3> )♦ (AL<3>  /RK<3) 

IK4)  =1 . 0  /  (  (  1  .0/00(4)  )♦  (  AL(  4  )  /RK  (4) 

0(5)  =1  .(}/((!  .0/10(5)  )*(AL(5)  /RK  (  k  ) 

U<6)  =1  .  (;/((  1  ,0/HO  (6)  )♦  (  Al.(  6)  /RK  (O 
DO  4  I = 1 *6 
(HI  )=|J(  I  )*S  (I) 

4  CONTINUE 

I  =  (  T  A  M  /  5  HI  .0 

PIN=P<  I  )  HP  (  1  +  1  )  -PCI  )  )*  (TAM -5*  (1-1  )  )/5 .0 

PAIR=RH  *PIN 

K«(T(11)/5M1.f) 

POUT  =P  (K)*((‘{K*  j)-P(K)  )  *  (T  f  1  1)  -5MK-1)  )/5 
301  CONTINUE 

DC  3D?  N*1,fc1 

W  A  RM ( V )=C.O 

COLD (V ) =0.0 

IF  (  F  (  N  )  )  31^311, 311 
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311  r<N)~n.o 

3  1  (1  C  0  N  T  I  M  ti  f 

F  RROR<  N  )  -  T  <N  >-TSf  T  <M  >-*R  ATC  <N  W  <N> 

IF  (  £R  k  0  R  ( N  )  )  7.)3,3'12,3U4 
703  COLfX'O  »-f  RROR(N) 

GO  TO  SOP 

304  WARMIM) =ERROR(N) 

302  CON  TIM  lit 

WA  RMS-0 .0 
COLDS=0.i* 

DO  303  1-1,6 

K  -  10  *  I 

WARM?- WARMS  t  WARM  (K  )  *$<  I  N  K  <  I  ) 

COLDS-C  OLDS  ♦  COLD  (K  )  *  S  M  ’  V  R<  T  ) 

3 05  CO M  T  I V  U  f 

SVf  AT=CSW*tRR0R(1)tSSW*  fVARMS  -  COL  DS  )tp  yj*  WAR  A  <  1  )  *  WARMS 
OIL  A  T  =  C  OIL  *  F  RDOR ( 1 )tS  DII  •< w A  RMS  -COL  DS) *PD I L*  V  4RM( 1 )  *WARMS 
STRIC-  -CCON*  E  k  R  0  R  ( 1) ~  S  C  0 N*  (  UARMS-COLDS  )  t  p  C  0  M  »  r  0  L  D  (  1  )  *  C  0  L  D  S 
C  THIS  CONTROLLER  RASED  0‘1  f  II  LI  SUIT  DATA 
C  CHILL*  .  9  «  70  3?  5 II- R.  V?.9S  1  66  *  ( 7  <  1  1 )  -  TS  f  T  ( 1  D) 

C  W  .6S7S  669?  *  (T  <6.i)-TSE  T  (60)  )  -4  .  Sf>4*>341  .!*  C  T  (5  C)-TStT  (50  >  ) 

C  t  -  S .  361  ?  6  fc  1 9  *  <T  (  3  0  >  -  T  S  E  T  (  30)  > 

C  THIS  CONTROLLER  HA  STD  ON  PRELIMINARY  DATA 

CM  II  L=  ? .001 -S.R??*  (T  (SCD-TSE  1<  Sfl)  )  -  2.4  0/*  <T<  70)  -  TSI  T<  70  )) 
f-77. 33  H  *  (1  ( 1  D-TSI.T  (  11  )  ) 

I  M  SWEAT)  3  0  9,3  12,312 
7  09  SM  A  1=0.0 

31?  IF (DTLA T)  713,314,314 
317  OILAT=n.:i 
314  IF  (STP  IC  )  31 S, 7  16,316 
3  1 S  S  T  K I C 1  0.0 

316  IT  (CHI  L  )  31 7,71 p, 31 H 

31/  CM  III  =  9 . 0 

3  IK  CONTIMUf 

4  (0  CONTIMUF 

DO  4  01  1=1,6 

N  =  1  0  *  I  -9 

o(N)=an(N) 

H  F ( N ) = P  rF(N  ) 

F  (  fv  )  =  D  .  0 

0  (  N  ♦  1 )  -  OH  (N  ♦  1 ) 

BE  (N  +  1  )  ~ D F B  (  Nt  1  ) 

F  (  N  ♦  1 )  =  G  .  0 
Q  ( N  ♦  2)  -  Q  B  (  N  ♦?) 

F)  F  (  N  ♦?  )  =BFU  (  N  +  ?  ) 

E  (Nt?)  =0.0 
0  (Nt  3)  =0P  (N  ♦  3) 

BF (N+3 ) =RFH ( Nt3 ) 

F  (N+3)  =0.0 

()(!(♦  4)  =  <JP(Nt4)t  WO RK M  (  I  )  *WORK  1 ♦ CHI L  M  (I  )  *CHILL 
B  F  (  Nt4  )  =BFFT  <  Nt4  )  tG  (Mt  4)  -  U8(  N  t4  ) 

F (Nt4) =0.0 

O(NtS)  =0F<(Nt  5)  tpORKM  (  I  )  *  WO  R  K  It  CHl  l  M  ( I  )  ♦CHILL 

BE  (Nt$  )  =HFH  (  f.tS  )  to  (M  *5  )  -OD  (  M  tS  ) 
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F  <  N  ♦  5 )  =  0 . 0 

Q(N*#»>  =  «B(N+6>  +  W0K’K"1(  I  )  *WORK  l*  CHILMU)  *CH  ILL 
PF(N*6)-PrR(N*6)+0(9*6)-on(N*6) 

F (N*6)  -0 .0 

G  (  N  ♦  7  )  =OH(N+  /)*WORK*1<  1  )  *  W  0  R  K  !  +  CHILH(I)  *  C  H  ILL 
PF(N*7)=BFB(N*/)*O(N*/)-GB(0*7) 

F  (W+7)  »P.O 
Q  (  N ♦  8)  =•«  B  CN  ♦  **  ) 

BF  (N*B >  =  PFB ( N*8  ) 

e (n+p) =o  .o 

Q  <  N  ♦  9)  *OP(N*  V) 

0r<N*9)  =  ((BFP(N+9)*SKlNV(I)*f>ILAT)/(1.n«SKISC(])* 

&  STR1I  C  )  )  *c  .  0  **  <  ERROR  (N*  9)/6 . 0) 

E (N  +  9) - f '  . 0 

4  n1  CONTINUE 

IF  (AWORK  <L)  -  <SA*HM)  )  14,14,15 

14  (3VJ0RK=S  A*BM 
FO  TO  16 

15  BWORK=  AWORK  (  L) 

16  CONTIMIIfr 

VF  =  <  22  .  9*BW0RK  )  /  4  .  g  3 
VFMT M= V [  / 6 H 
BP  =  763 . 0 
f  1  =TAI  K *  ?73  .  b 
T  E  =  T  ( 1  1  )  +  ?  7  3  .  P 

AIRDCP- (n.?9''S*  (HP/76P)  *(?9  7/TE  ))/mf:0 
H?  '-OUT  =VF*  (  <  BP *9 9. 96?)  /  (6?.  4*TE  >>*  (POUT/  (  BP-  POUT)  >* 
ft  (18.31534/99.06?) 

RMrAT-(6.7F. -7*<VEMM**?))-M6.?r-5*(VF'11V**3)> 

RHL=(VE*AIRDCP*(  Tf-TI)  )♦(().  5o*H?00UT) 

DO  6  ■'?  1=11,14 

F  ( I  )  =(  RMl.-U  .  116*  RUE  AT)  /  4  .  T 
6  0?  C0NT19UF 

DO  503  K-1,6(! 

RC ( K )  =  B  F  ( tO *  (T  <  K  ) ~  T ( 61 )  ) 

T  D (K )  =  T  COO  *  (T  (K)-T  (  K  ♦  1  )) 

5  00  CONTIMtIE 

00  5M  1  =  1/6 
K  =  1  "  *  I  -9 

H  F  (K)  =  Q(K)-F.  (  K)  -HC  (K)-TD(K) 

H  f  (  K  ■*  1  )  =  0(K  ♦  1)-E  OC  +  1  )-B  C(K*  1  )-TP(K*1)*TI>(K) 

HF(K4?)=0(K  +  2)-F(K*?)-BC(K4?)-TI)(K+?)+TD(K  +  1) 
HF(K*3)=«(K*3)-F(K  +  5)-HC(K*3)-Tr>(K  +  3)*-TD(K+?) 

H  F  (  K  ♦  4  )  =  G(K*-A)-t-(K*4)-BC(K*4)-TD(K*A)*T0(K+3  ) 

HF(K  +  5)=0(K  +  5)-r  (K*5)-BC(K*5)-7D(K*r04TD(K*4  ) 

MF  ( K  *6  )  =  G(K  4  6  )  -  F  (K  *6 ) - P  C  (K  ♦  (  )-TP(K*6)*Tr'(K*t'  ) 
HF(K*7)=G(k*?)-E(K.*7)-nC(K*7)-Tf>(K*7)*1D(K*6) 

HF  (ft  *3  )  =«<K  ♦  i  )  -F_  (K  ♦»  )~F.  C  («♦  h  )-  TC  (K  +  6>*  TD  (  K  •*  /  ) 

HF(K49)=0(K4c/)-E(K  +  ))-BC(K+9)+TD(K4vS)“MSS(I)*(T(K>9)-T(I*61)) 
HF  (1*6  1  )=HSS  (1  )  *  (  T  (K  +  9  )  —  T  <  I  *61  )  )-H  (  I)*  (T  (  1+6  D-TA1P) 

501  CONTIVUr 

H  F  (  6 1  >  =  0  -  0 

DO  5n?  K  =  1  #60 

MFC  61)  =MF (6 1 )*HC (K) 

5(1?  FONT  F  9  II  F 
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I '•(  M,WU'.  mt — *U 


6  j1 

6  30 


7jG 


70 1 


fc  'JO 


S"? 


8  1 


9tn 

999 

91? 


911 

915 


3C  1  ,  701 , 7Q1 


i>  T~0 .01  66^ 66 .‘'6ft (S6 6  7 
00  6  )3  K  =  1  #  6  / 
f  <«0=MF  (K)/C  <K  ) 

AU-ABS ( F  (K)  > 

i f  (  au*  ot-c.  i )  <sir„<so% f.ui 

PT  =  0.1  /AC 
fONTnUF 
0  0  7  0  K  =  1 »  6  7 
T f<)  -r  ( K )  +  r  ( K)*  01 
CONTIVUF  ' 

TI>1F  =  T  I  MMDT 

liimf  =  ap.*timf 

IF  <L  11  VF~J  NT -  IT  JM£) 

comtiiuf. 

I T I M  F  a  I  TIKE*  IN T 

CO  =  0  -3 
H  P  =  Q  .  3 
E  V  =  0 . 3 
T5  =0 .3 
TF5  =  0  .3 
H  FLOW*  0  .r* 

EOF-  0.  0 

DO  F  '<0  N  -  1  ,  f,  '' 

ro  =  co*pf  <u) /ec.  o 

liP'HPf  0  (  N  ) 

F.V-FV+E  (N> 

CONTIVUr. 

Do  8-2  1=1,  ft 

SPF  =  s«ir+Pr(KM  )/6  ).n 
Ts-rs+r  (in*i  )*r  (ic*t)/cs(i) 

CONTI'JUF 

DC  8  !1  N  =  1 , 6  1 

7P=T«*  T  (M)*  C  (N)  /C')(1) 

H  FLOWS  H  T  L  0  W  >  H  F  (  M  ) 

roNTisiir 

FV  =  f  V/SA 

ht^h0/  sa 

H  FLOIV-  H  FLOW  /  SA 

IF  =  E(11)  +  f(1P)4f<i3)4P(Uj 

C()\f)?=  (MP-FF  /S A -HF LOW)  /  IT (61  J-TS) 

I  F  ( I  TI  MF -iNT  )  910,0^0,911 
URITFC 7,' 99) 

E<RMAT(M*,1rX,‘HAYLFr  FVC  C  R  •) 

WRITE(3,91?) 

8'TLFG., -CO*, 2x,.crND.;?;/S 

NN  =  ■ 

DO  ?  ;5G  1=1  ,t 


C-20 


qgyr^gqqp  gyygyg^vwt^ 


t}  F3p»:«|wvs,«S’’ 


2  05  0 


C  1 1  3  0 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


1101 

cm? 

c 

C3001 

c 

c 

C  3  id  ? 

in? 


180  . 
39  . 

9  . 
21  . 


ro.MD3=  C0ND3  ♦«(  I  )*(  I  (  1  *6  1  )-T  A  IR  ) 

CONTI  \i  Uf 

COND  3-CONb  3  /  SA 

H  THAI  =  H  P  -  uf  L  nw-f  6  /  S  A-C  ON  0.3 

\  '1  -  N  N  ♦  1 

KK=XK*  1 

A  (  K  K  )  =  1  T  I M  fc 

A  (  J  J  ♦<  K  )  =T  A  I  R 

A<?»JJ*KK)  =  T(11> 

A (y * JJ ♦ KK)=  T  S 
A  <  4  *  J  J  ♦  K  K  )  =  T  (1) 

B  <KK  )=  I  1  IMt 
B ( J  J ♦ <  K >  =  T  A  I k 
P(?*JJ+MO  =  H»’ 
n<3*JJ ♦KK)=SLF 
B  (  4  *  J  J  )~CO 

JT  IMt  =  J  1  IMF  ♦  1NT 
T  f ( J  TT  ME -30 )  301 ,1 1 01,1  m 
J  T  I  K  e  =  0 

1  F  (  ITIMF-VU  )  10?  *  1  Ki?  r11  C  ? 

CALL  PLOTLfc  (  1,A,4  6,5#4H,0) 

MR  1  T  E  (3/^n1) 

FORM  AT  (  *  -*  ,  1  >'x,  1  1  :  T 4 T  H  '  /  C>X /  *  ?=■  T  R  V  5  X/  *  3=  T  S  *  »  S>  /■  ' 
CALL  PLOTEf  <  ,*#U,4K,5#4£m> 

WRITE  (M  f0  2> 

FORMAT  (  '-'/iry,  •  1=  T  JT  R  *  ,  m  *  ?  =  M  5X  ,  '  3  ‘SB  F  ’  ,  ‘'X,  • 

0  F  T()R*J 

ENO 


4000069  3  00. 0  0 HU 

4  or,  no 

5  5  0  U  U  4.70800  n.mou  n  .4  8  400 

OOOtO  0.71  0  k  C  0.06580  0.81620  0.35120 


4  r.  T  H  ’  ) 

4  -  C  0  *  ) 


l>rJBT»»VT?»’»Trff7,f  JfKVFV^OF?7 1f-f??*,‘rJ*ri''f  3* 


H 10  Y  S. 6.  =  1.0846?  X  BTDY  FAT=  0*0711  SURFACE  A  RE  A  =  1.66I5SO  IN 
ARIPQSt  TISSUF  W1.=  4.9797  K  C.  N0N-AD1P0SF  TISSUE  WT.=64.3700  KG 


HEAD 

RAD  IUS  = 

0 

.09  9  2 

M 

HEAD  VOL. 

=  C 

.  0  03° 

M 

CU  OF 

TRUNK  LENGTH 

- 

0.97  ('9 

THUNK  VOL 

m  - 

0.03  77 

M  C  URE 

a  sr 

L  t G  T  H  = 

0 

76  3  5  M 

ARK!  V  9t.  = 

0. 

0066  3 

CURE 

HAND 

L  ENG TH  = 

1 

.18  2  3 

M 

H  AND  VOL  . 

=  0 

.  000  A 

M 

CU  PE 

l  EG 

LENGTH= 

1 

„ 

4  42  4  M 

L  F (-  VOL.  = 

c. 

0197  5 

CU?  f 

FOOT 

LENGTHS 

1 

.  54  9  3 

M 

FOOT  VOL. 

=  0 

.0009 

M 

C  J  PE 

FOLLOWING  ARE  SCI)  VALUES  IN  SO  M 
0.1  398  0.7  1  61  0  .  ?  S  2  5  1.1409  9  .6  5  43  0.1905 


FOLLOWING  ATE  AL  (  I  >  VALUES  IN  M 
0.3  32  8  0.30  31  0  .0031  9.0028  0  .0C>71  C'.DfJ?9 


FOLLOWING  ATE  DO  FI  )  VALUES  IN  M 
0.2  110  0 .  ?  ?  4  9  0.11  79  0.C38O  0.1  4  45  0.039  2 


fOLLOWING  ARE  HSS(I)  VALUES  I N  K C A L/HR / D FG  C 
1.5  430  8. 9372  3.33'9  1  .2735  7  .8793  1.7441 


following  atf 

TC  (  I  ) 

values 

IN  KCAL 

/HP  DEG 

C 

0.61 

1  .6  8 

2 . 92 

6.99 

3  0.33 

31.63 

32.93 

13.73 

8  .  6 

o.n 

3.63 

7  .64 

11  .60 

14.14 

16.57 

29. 4P 

24.38 

22.29 

27.61 

c.c 

2.0  6 

6.03 

9 . 12 

10.05 

11.17 

14.25 

17.32 

1  ° .  0  8 

18 . 2  3 

0.0 

4.4? 

9.3  3 

14  .  12 

?  a .  1 3 

71,10 

75.84 

8  0.59 

20.1  5 

7.  19 

9.0 

5.  4(1 

11.39 

1  7  .  23 

19.18 

21.43 

27.24 

3  3  .  n  4 

4  0.11 

41  .88 

n.o 

5 . 8  0 

1  2  .2  5 

13.51 

4U.14 

1  50.35 

15  6.56 

162.78 

27.71 

1 C  .  4  8 

9. 0 

FOLLOWING  ARE 

A  (I  )  VALUES  IN  SO  M 

0.04 

0.0  6 

0  .  'ft 

J .  1  0 

0.  1  1 

0. 1  1 

0.  1  1 

0.  1  1 

0.1? 

(.13 

0.  21 

0.3  0 

0.  37 

0.43 

0,49 

0.54 

9.59 

0.63 

0.66 

".68 

G .  C  7 

0.1  fj 

0.  13 

0.15 

o,  i  y 

0.20 

0.2  2 

0.2  5 

C  .24 

C.  26 

0.0  3 

0 .0  4 

0 . 05 

0.16 

0.0  6 

0.07 

9.07 

0.07 

0.08 

0.10 

0.  1  7 

0.2  5 

0.  31 

0.36 

C.42 

0.4  7 

0.52 

0.56 

0 . 5  8 

0.61 

0.0  4 

p  .06 

0 .  98 

0.09 

0. (JO 

0.1  n 

9.19 

9.10 

0.11 

0.14 

FOLLOWING  ATF 

RMP (  I) 

VAL  IES 

1  N  W 

0.055? 

0.071  5 

C.C82? 

0.0891 

0  .0920 

0.0931 

0.0942 

0 

.0954 

n 

.0971 

0.1004 

0.03  3  8 

0.0491 

0.06  04 

•0.0703 

o  .■  79d 

0 .0887 

0.0^68 

J 

.1037 

n 

.  1  0  8  4 

0.1122 

n.01 4  9 

n .  9  21  7 

".0267 

"  .0313 

n  .  "3  6  1 

9.9498 

0.  "4  5  0 

0 

.0484 

.0507 

9.053? 

0 .004  1 

0.005  9 

r>.90  73 

c-.rC8  3 

l . n08 ? 

0 . °090 

C. '093 

0 

.0097 

r, 

.  0  1  0  8 

0.0139 

0  .01  9  1 

0.027  7 

r  .03  40 

r  .0  399 

0  .0460 

0.C51V 

0.0571 

0 

.0613 

n 

.0  6  3° 

0.'  670 

0.0046 

0.006  7 

0.00  *2 

f  .009? 

0 .909  7 

0  . n099 

0.0101 

■J 

.0105 

L. 

.0115 

C  .  i  1 4  5 
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FOLLO  WING 

AR  E 

XU)  VALUES  IN 

M 

1.02v  2 

n. 

01?? 

0.0093 

0.0  04  6 

0 

.001  1 

0.001  1 

0.  nr-i  1 

0  .0014 

0  .0021 

0.004  6 

0.01  81 

0 . 

0 1  ?  5 

0.01  Cl 

0  .  ;j  G  O  7 

'j 

.0094 

0.0084 

0. 0077 

0  .0061 

C  .  C  0  3  4 

0.004? 

n .  0  3  '■  0 

c . 

0  05  5 

0.0  0  AS 

0.0048 

0 

.0050 

0  .0044 

0  .  C  (  4  0 

3  .0028 

C  .0019 

0,004? 

0.002? 

~  ft 

j  11  5 

0.0012 

0. -0  0  7 

0 

.COO'S 

0.*'001 

0.  •■001 

a  .0006 

0.00  1  6 

0 . GO 4 6 

0.01 .1? 

0 . 

•1 37  0 

0.0057 

0.OG6- 

n 

.  006  2 

0.005  5 

O/'^n 

0  .0033 

n  .  0  0  ?  C 

0.904  3 

0.03  2  5 

J. 

001  7 

G .  0  U  U 

0.0  0)7 

U 

.000? 

0.000? 

0.  00? 

3  . HOC 6 

0  .001  5 

0.004 4 

FOLLOWING  ARE 

THE  VALUES  OF 

C  S(1 ) 

AND  CB  <  1 ) 

csn  )  = 

3. 43  a 

CH  <1  )  = 

56.34  7 

FOLLOWING  ARE 

C  (I  )  VALUES  IN 

K CAL  /  DEG-C 

0 .564 

G  .56  4 

0.5  64 

0.564 

0.085 

0  .085 

0.385 

r' .  3  8  5 

0.099 

0.24  5 

2.506 

2.506 

2  .  5  C6 

2.5  06 

4.105 

4  . 1  C  5 

4.105 

4.105 

1.87? 

1  .234 

0.3  6? 

0.  36  2 

0,  362 

0.36? 

0.77  2 

0.77? 

O.7?? 

r .  7  7  2 

0.256 

0.443 

0.036 

0 . 05  6 

0.036 

0.036 

".01  7 

9  .01  7 

0.017 

0.01  7 

0.959 

0.170 

1.077 

1  .  07  7 

1 . 0  77 

1  .  J  7  7 

2.33  1 

2.331 

2.331 

2.  331 

r .  6  3 1 

1.09  7 

0.36  1 

0.36  1 

0.0  61 

0.  G61 

0.01  7 

0  .01  7 

0.017 

C.01  7 

C.059 

0.218 

2.250 

0.05  0 

0.265 

3.102 

0.04  5 

0.238 

FOLLOWING  ARE 

09(1)  1 

VALUES  I 

*1  KCAL/HP 

3.167 

3.167 

3. 1  67 

3.167 

C  .03  1 

0.031 

0.031 

n.031 

0.04  9 

0.082 

11 .364 

11  .36  4 

1 1  .  3  64 

1 1  .364 

1  .502 

1  .502 

1.5  02 

1.502 

0.937 

0.41  1 

0.151 

n.15  1 

0.1  51 

0.151 

0.28  3 

0.28  3 

0  .2  8  3 

0.283 

0.128 

0 . 14  8 

0.31  8 

3.01  8 

0.018 

0 .016 

0.03  6 

0.006 

0 .306 

C.006 

0.020 

C  .05  7 

0.4  66 

0.46  6 

n .  4  66 

0.466 

3.855 

9  .85  3 

0  .  *•  5  3 

0.853 

9.315 

0.366 

0.3  29 

C.02  9 

0 . 0  29 

0.02  9 

0.00  6 

0.006 

0.006 

0.006 

n.  0  30 

0.073 

following  are 
11.25°  11  .259 

BFb (  1 ) 
11.2  50 

VALUES 

1  1  .25-. 

IN  L  /  H  K 
0.0  3  7 

• ;  .057 

0.037 

0.03  7 

",  0 5 9 

1.454 

52.5 V 3 

5?  .  50  0 

52.  5  9n 

c 

r*- 

w’”"* 

ft 

1  .89  3 

1  .803 

1  .  K  9  3 

1.80? 

1.124 

2.13** 

3.181 

fj .  1 8  1 

0 .161 

0.161 

0.33  9 

0.33  9 

0.739 

0.339 

0.154 

9.66  5 

0.321 

3 . 02  1 

0.0  21 

G.G21 

3.037 

r.  .007 

0 . 0  0  7 

1.30  7 

C  .  0  ?  4 

1  .05  6 

3.5  60 

0. 56  0 

0 . 5  60 

0.560 

1  .324 

1  .024 

1  .024 

1.C24 

0.379 

1.074 

0.335 

'■>  .03  5 

0.035 

3.035 

9.00  7 

■'  .00  7 

0.037 

0.007 

0.035 

0.56  7 

FOLLOWING  ARE 
0.35  7  0.07  ? 

R  ( I )  VALUES  IN 
0.0 h?  0.091 

C 

0.09? 

0.09? 

0.994 

0.  J95 

0.097 

0 . 09  9 

0.  33  5 

0 . 34  9 

0 .  0  61 

C.  070 

0  .08  0 

0  .06  9 

3.097 

0.104 

!",  1  9P 

0.111 

0.31  5 

0.0?  2 

0.927 

0.031 

0.9  3  4 

0.04  1 

0.94  5 

0.04  9 
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9.053 

0.3  04 

0 . 00  6 

0.0  07 

0.0)  P. 

0.90  9 

9.009 
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0 . 0  1  n 
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0.320 
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3.01  0 
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C.1  05 

0.11? 
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0.020 

FOLLOWING  ARE 

CN  (I  ) 

VALUES 

IN  M 

0.04  5  2 

0.9  65  4 
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0.04? 8 
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0.064 5 
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1.  Specific  gravity  of  the  simulate  subject  may  be  entered  by  data 
statement.  If  SG-0  the  program  wl II  calculate  SG  based  on  subject's 
height  and  weight. 

2.  Protective  device  thickness  (THWS)  In  meters  is  entered  by  data 
statement . 

3.  Protective  device  thermal  conductivity  for  each  body  segment  (RK (N) ) 
Is  entered  by  data  statement. 

k.  Water  temperature  (ATAIR)  Is  entered  by  data  statement. 

5.  Subject  height  and  weight  are  read  from  the  first  data  card. 

6.  The  output  may  be  labeled  by  simulated  experiment  by  altering  the 
hollerith  stream  In  format  statement  999. 


